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SUMMARY 

The object of this thesis research- has been to 
develop an improved broad band impedance matching technique. 
The technique is capable of resolving points in the wave- 
guide which generate reflected energy. A version of the 
comparison ref lectometer has been developed and fabricated 
to determine the mean amplitude of the reflection coeffi- 
cient excited at points in the guide as a function of dis- 
tance, and the complex reflection coefficient of a specific 
discontinuity in the guide as a function of frequency. An 
impedance matching computer program has been developed which 
is capable of impedance matching the characteristics of each 
disturbance independent of other reflections in the guide. 

The characteristics of four standard matching elements 
have been compiled, and their associated curves of reflection 
coefficient and shunt susceptance as a function of frequency 
are presented in Appendix A. The characteristics of the 
four standard matching elements (symmetrical capacitive iris, 
asymmetrical inductive iris, inductive metal post, and di- 
electric post) are also programmed into the impedance match- 
ing program listed in Appendix B. 

A sample of each of the four standard matching ele- 
ments have been fabricated, and results measured by both the 
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comparison ref lectometer and a slotted line setup have been 
compared with theoretical predictions. The measurements 
made by the comparison ref lectometer on the sample standard 
matching elements have shown agreement within the accuracy 
of the comparison ref lectometer . 

A shunt slot radiator was fabricated, and its charac- 
teristics were measured and impedance matched to demonstrate 
the computer aided broad band impedance matching technique. 
The shunt slot radiator was impedance matched over a 2 GHz 
bandwidth to a VSWR of 1.16 to 1.0. A theoretical dimen- 
sional perturbation has been made on the matching element 
to determine the extent of the sensitivity of the impedance 
match to machining tolerance. 

As a result of this research, an economical, fast, 
and reliable impedance matching technique has been estab- 
lished which can provide broad band impedance matches. 
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CHAPTER I 
INTRODUCTION 
Motivation for Problem 

For many years impedance matching at spot frequencies 
has been practiced, and the Smith chart is a well known tool 
for narrow band impedance matching. In most cases the 
matching element is placed some distance from the disturb- 
ance to be matched. In such a configuration the phase 
angle of the generalized reflection coefficient resulting 
from the matching element and the discontinuity to be 
matched varies as 2 8L, where 6 is the propagation constant 
and L is the distance between the matching element and the 
mismatch. The propagation constant can be expressed in 
the form, 

8 = (o) 2 y e - K 2 ) 1 * = 2i\/Xg (1) 

where: Ag is the wavelength in the waveguide, K is defined 

for dominant mode propagation in the rectangular waveguide 
by K = it /a , a being the width of the guide, uj is the 
angular frequency in radians per second, y is the permea- 
bility, e is the permittivity. 

Equation (1) shows that 8 is frequency dependent and 
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as a result complete cancellation of waves generally 
extends only over a narrow-band. The greater the distance 
L the narrower the bandwidth of the impedance match. If 
the distance L were reduced to zero then there would be no 
variation in phase angle due to the electrical length sep- 
arating the discontinuities, and two discontinuities having 
reflection coefficients of equal magnitude and opposite sign, 
Appendix C, would present a broad-band match. 

In recent years there has been a need to use many 
microwave devices over a large portion of a waveguide band- 
width. As an example, some antenna systems which transmit 
swept signals require matching over a large percentage of a 
waveguide band. All components which make up this trans- 
mitting system and handle the microwave energy, such as 
ferrite circulators and phase shifters, couplers, modula- 
tors, and many other components, must also have a large 
band pass. Without good impedance matching over the band 
of interest, large quantities of energy are reflected back 

towards the generator and must be absorbed in isolators or 

c 

attenuators. As a result, excessive amounts of energy are 
demanded from the generator in order to supply the load with 
the desired energy. Many times these excessive energy 
demands are economically prohibitive and never desirable. 

Usually broad band impedance matching is done with- 
out much insight into the unmatched component. The loca- 
tion of reflections and the number of these reflections have 
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in most cases been left unknown, exceptions being only in 
the simplest cases. Therefore, impedance matching of com- 
plex components has been largely trial and error. The end 
result is often neither very satisfying nor very economical. 

The location of the discontinuity in the waveguide 
and the complex reflection coefficient of the disturbance 
as a function of frequency need to be known in order to 
gain insight into the problem of impedance matching. In 
this work a version of the comparison ref lectometer ^ is 
used to determine these parameters within limits. In some 
cases there may be some difficulty determining a disconti- 
nuity which will identically cancel out the effect of the 
original mismatch. However, the impedance matching tech- 
nique described provides the engineer with a powerful tool 
to obtain broad band impedance matches that are practically 
acceptable . 


Definition of the Problem 

Discontinuities in an otherwise uniform transmission 
line cause the excitation of a reflected wave in order to 
satisfy the boundary conditions on the electric and magnetic 
fields at the discontinuity. The reflected wave generally 
represents a reduction in transmitted power and is described 
as arising from a mismatch in the characteristic impedance 
of the line. An additional discontinuity placed in the 
transmission line in such a way as to cancel the original 
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reflected wave is said to impedance match the transmission 
line. This general technique provides a good match for 
lossless mismatches if the reflection coefficients of the 
mismatch and the matching discontinuity represent waves 
which cancel identically and if the discontinuities are 
located equal distances from the generator. If, however, 
the mismatch and matching element are not located the same 
distance from the generator, the match will be frequency 
dependent, and therefore of a narrow band nature. 

Most mismatches that occur in practice are in fact 
reflections arising from multiple discontinuities and gen- 
erally cannot be broad band impedance matched by a single 
matching element. Successful matching of multiple reflec- 
tions requires that the location of each reflection should 
be known and that the contribution of each reflecting point 
to the total reflection coefficient be known. 

The objective of this research has been to develop 
an improved broad band impedance matching technique. The 
technique is capable' of resolving points in the waveguide 
which generate reflected energy. A version of the compar- 
ison ref lectometer has been developed and fabricated to 
determine the mean amplitude of the reflected wave excited 
at a point in the guide as a function of distance, and the 
complex reflection coefficient of a specific discontinuity 
in the guide as a function of frequency.^ An impedance 
matching computer program has been developed which is capa- 
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ble of impedance matching each disturbance independent of 
other reflections in the guide. The justification for the 
impedance matching procedure is given in Appendix C. 

The characteristics of four standard matching ele- 
ments have been compiled, and their associated curves of 
reflection coefficient and shunt susceptance as a function 
of frequency are presented in Appendix A. The curves take 
the form of plots of shunt susceptance, and plots of mag- 
nitude and phase angle of the current reflection coeffi- 
cient. 
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CHAPTER II 

BACKGROUND OF IMPEDANCE MATCHING 
AND OF REFLECTOMETRY 

Impedance Matching 

Reflected waves are generally undesirable and are to 
be avoided in waveguides. One method of minimizing reflec- 
tions in a waveguide is to design the system such that the 
load impedance will completely absorb the incident fields 
without reflections. This load corresponds to a character- 
istic impedance termination in a transmission line. A 
second approach to the problem is to create a reflected wave 
near the load that is equal in magnitude but opposite in 
phase from the wave reflected by the load. The two reflec- 
ted waves therefore cancel each other as shown by Equation 
(C-14). Both methods of impedance matching are usually 
used simultaneously. The system is initially fabricated so 
that the load provides as good an impedance match as is 
possible to obtain with a reasonable effort. The reflected 
wave that still remains is eliminated by the use of an 
impedance matching system that creates a compensating re- 
flection. Many waveguide arrangements have been devised 
for generating a controllable reflection. Some of these 
are analogous to the impedance matching arrangements em- 
ployed in transmission lines, while others are unique to 
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waveguides . 

The shorted stub so commonly used in transmmission 

lines has a waveguide analogue, the stub guide or E-H tuner. 

In Figure 1 there are two types of tuners, the series tuner 

and shunt tuner. The series and shunt tuner can be used to 

impedance match a discontinuity or mismatched load over a 

narrow band of frequencies. To achieve this match the 

length of the stub as well as the distance from the load are 

critical dimensions. An analysis of this matching procedure 

2 

is given in many electromagnetics texts. 

A second important procedure used frequently to im- 
pedance match over a narrow band of frequencies is the 
quarter wave transformer. The impedance transformer con- 
sists of a length of line of impedance Z t that begins a 

quarter wavelength from the load, as indicated in Figure 2. 
The proper impedance of the quarter wave transformer Z t is 

3 

given by Equation (1) as : 


When the impedance Z t satisfied Equation (1) the 
impedance looking to the load at point P will be equal to 
the characteristic impedance of the line, Z Q . It should be 

noted that the impedance matching just described is exact 
only for the frequency corresponding to the waveguide wave- 
length A g 0 . Once this match is achieved at a particular 


shorting plunger 



_ - of side arm 

□— O 


(a) Series T 



(b) Shunt T 


Figure L. Series and Shunt Tuners 
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-o 

P 



Figure 2. Quarter Wave Transformer 


Ago, the match does not hold for a general A g ^ Xg Q . This 
results in a narrow band impedance match. There are many 
such techniques giving restricted matches. 

Ref lectometry 

The basic function of a ref lectometer consists of sanv 
pling an incident wave and a reflected wave and determining 
the ratio of the two. The ratio of the reflected to the 
incident electric field is defined as the voltage reflection 

coefficient r , and has both real and imaginary parts. 

4 

The general procedure uses two directional couplers 
located back to back as in Figure 3. The incident energy 
is coupled into the directional coupler D1 and the re- 
flected energy is coupled into the reversed directional 
coupler D2. Each directional coupler is terminated in a 
0 square law detector. These detectors provide a D.C. 

voltage output which is proportional to the square of the 
electric field intensity. If the incident detector is used 
> to provide a feedback signal which enables a high frequency 

oscillator to maintain a constant power level, then the 
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magnitude squared of the reflection coefficient is pro- 
portional to the output of the detector D2, Equation (2) : 

2 

| r I = T E r , (2) 

where T is a proportionality constant, and E r is the 

voltage output of the square law detector. 

Alternatively, if the output of both detectors are 

used as inputs to a ratio meter than the reading of this 

device is a function of the magnitude of the reflection co- 
. . 5 

efficient squared. The ratio meter can be calibrated to 
provide magnitude of reflection coefficient or voltage 
standing wave ratio (VSWR) . The use of a swept oscillator 
in either of the above two cases will provide a display of 
the magnitude of reflection coefficient as a function of 
frequency. 

It is well to note that these techniques do not 
provide phase information. Since the phase information is 
unknown, there is no assurance that two discontinuities will 
in fact provide a cancellation of waves even though the 
curves of the magnitude of the reflection coefficient |r| 
are equivalent. Therefore, impedance matching using these 0 
ref lectometers is largely trial and error. 

Reflecto "meter" 


A technique which takes a different approach to 
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ref lec tome try was introduced by F. C. deRonde. This 
system called the reflecto "meter" utilizes three detectors 
in the main transmission line, Figure 4, instead of the two 
directional couplers as in Figure 3. The advantage of this 
device is improved accuracy of the modulus of reflection 
coefficient, but the setup is more difficult to fabricate 
than the ref lectometer composed of directional couplers. 

The common disadvantage is that it gives no phase informa- 
tion and cannot resolve individual discontinuities in a 
waveguide . 

It is desirable when impedance matching to have both 

the phase information and the magnitude information of the 

reflection coefficient, as indicated in Appendix C. This 

information can be obtained using slotted line tech- 
7 8 

niques, ' but this procedure is extremely time consuming 
when dealing with the many data points needed to cover a 
large frequency range. In addition, the slotted line 
technique has no provision for the separation of reflectioh 
coefficients generated by different discontinuities. 

Time Domain Ref lectometer 

Time Domain Ref lectometer (TDR) is an approach which 
can determine the location of discontinuities in a trans- 
mission line. TDR is an application of a pulse reflection 
9 10 

technique. ' A pulse of energy is transmitted down a 
transmission line. If there exists an impedance discontin- 
uity in the transmission line, energy will be reflected 


unknown 

reflection 



measuring 

detectors 


Figure 4. Reflecto "meter" 
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back toward the generator. The location of the discontin- 
uity can be determined by measuring the round trip time de- 
lay for reception of reflected pulses. The narrowness of 
the pulse and the quality of the measuring equipment de- 
termine the resolution capabilities of the TDR. The narrow- 
band TDR is used when measuring waveguide structures The 
amplitude of the reflected pulse is related to impedance, so 
any slight deviation from the characteristic impedance level 
of the output of the TDR can easily be recognized and 
measured. 

Time Domain Ref lectometry has been applied to open 
transmission lines, coaxial cable and recently to waveguide 
structures. For the waveguide structures a step modulated 
carrier is used^ to remain within the confines of the wave- 
guide passband. Generally, the true value of the narrowband 
TDR is in applications where long runs of waveguide are 
used. The locating accuracy for this device is of the order 
of ±3 percent of the distance measured, depending on the 
quality of the leading edge of the transmitted pulse. 

Network Analyzer 

Recently, the network analyzer has been introduced.’*’^' 
This device characterizes microwave components in terms of 
their scattering matrix. This characterization allows the 
microwave device to be modeled and more efficiently used 
in complex networks. The network analyzer is an important 
impedance matching tool, as a result of its modeling ability. 
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However, the network analyzer lacks one important quality, 
it cannot measure the effect of individual discontinuities 
independent of adjacent disturbances. Therefore, the im- 
pedance matching of the composite discontinuities of the 
microwave device remains a difficult problem requiring a 
composite of matching elements. If however, the composite 
reflection coefficients of microwave devices could be broken 
down into an ensemble of simpler reflections, then it is 
reasonable to assume that the simpler reflections would be 
easier to impedance match individually. The.net result is 
an overall impedance matched device. 

Comparison Ref lectometer 

The comparison ref lectometer^ is another technique 
which provides the location of a disturbance in the guide. 

If there is more than one disturbance in the guide these 
disturbances are shown as separate disturbances and their 
locations are displayed. The comparison ref lectometer also 
can provide a curve of magnitude and phase of the reflection 
coefficient as a function of frequency for each discontin- 
uity in the guide. 

In using the comparison ref lectometer inaccuracies 
in the waveguide components are largely cancelled out due 
to the technique of taking measurements, as described in 
Chapter III. This cancelling effect enables a background 
reflection coefficient noise level typically less than 
0.00005 to be measured, as shown in Figure 31. Because 
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of these attributes a comparison ref lectometer has been 
chosen as the basis of the computer aided impedance match- 
ing technique, and has been fabricated to take the measure- 
ments used in this research. 
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CHAPTER III 

APPROACH TO IMPEDANCE MATCHING 
Introduction 

The computer aided broad band impedance matching 
technique, Figure 5, consists of first collecting data to 
determine the combined complex reflection coefficient of 
both a reference step and a test discontinuity. The lo- 
cation of the disturbance to be impedance matched is then 
determined and the complex reflection coefficient is cal- 
culated. The physical dimensions and the relative loca- 
tion of a preselected matching element are calculated by 
the impedance matching computer program. The result is 
the impedance matched disturbance. If additional matching 
is required, a second matching element may be selected. 

The physical parameters of this element are again calculated 
by the impedance matching program. 

Two sets of data are taken by the comparison reflec- 
tometer at 50 MHz intervals from 7.975 GHz through 
12.425 GHz. One set of data is taken with the reference 
step terminated in a matched load, Figures 6 and 7, 

The second set of measurements is taken with the reference 
step terminated by the test disturbance. These data are 
processed by two computer programs. The first, Test 1, 






















Photograph of Comparison Reflectomete 
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determines the mean magnitude of the reflection coefficient 
as a function of distance as shown in Figure 8. This curve 
locates discontinuities relative to the reference plane. 

The second program. Test 2, calculates the complex current 
reflection coefficient as a function of frequency at a pre- 
determined location in the waveguide, Figure 9. This com- 
plex reflection coefficient characterizes the measured 
disturbance in the impedance matching program. 

The impedance matching computer program developed as 
part of this thesis determines the physical dimensions of 
a matching element and an appropriate location with respect 
to the original mismatch. The program provides the best 
impedance match over the bandwidth of interest in the root 
mean square sense. This "best" impedance match is depen- 
dent on the matching element and the method of choosing 
its parameters. The input data for the matching program 
are the results of the comparison ref lectometer program. 

Test 2. All of these programs discussed are listed in 
Appendix B. 

Theory of Operation of a Comparison Ref lectometer 

I 

Introduction 

The comparison ref lectometer first introduced by 
D. L. Hollway,! is an instrument designed to locate and 
measure the characteristics of reflections in waveguides and 
transmission-line systems. It is particularly suitable for 
measuring small reflections in microwave components up to 


nfiGNITU 
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Figure 8. Typical Locating Plot 
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one meter in length. 

The principle of operation is the comparison of the 
total magnitude of reflection coefficient of a known re- 
ference disturbance and a test disturbance with that of the 
known reference disturbance. From this comparison the lo- 
cation of the test mismatch can be determined, and its com- 
plex reflection coefficient as a. function of frequency cal- 
culated. A high speed digital computer is used to perform 
the necessary calculations. 

£> 

Governing Equations 

Consider a waveguide system shown in Figure 6. A 
component having a single reflection |r,|exp[j0] is 
connected to a reference reflection having the scattering 
coefficients shown in Figure 1©. The reflection coefficient 
of the combination may be written, 

bl/ a l = ( s 22 d” S 33^1^ +S 32^1^23^/ ^ -S 33^1^ ^ 

b i /a 1 = ^S 22 - r 1 (S 22 S 33 - S 32 S 23 )J/ (1-S 33 T 1 ) . (2) 

Generally, the reflection to be measured will be 
located at some distance from the reference plane, so 
that the value of T 1 measured at this plane will be, 

r l = 1 r -L | exp [j(6-23L 1 )] . (3) 
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Reference 



Figure 10. Scattering Representation of a Single Reflection 

Connected to a Reference Connection 

The theory and computations are simplified, and the 

accuracy of the results improved if a reference reflection 

is used having a reflection coefficient which is essentially 

constant in magnitude and phase throughout the frequency 

band at a stationary reference plane. A design which has 

12 

been found by D. L. Hollway to be superior to others in 
this respect consists of a symmetrical E-plane taper, having 
only negligibly small reflections followed by a sudden step 
back to full guide height. 

If we ignore for later simplicity the small correc- 
tion required for the step capacitance, the step reflec- 
tion, shown in Figure 11a, may be considered as a lossless 
transformer^ - , having turns ratio N, set in the guide at 
the reference plane as in Figure lib. 

Let T r be the voltage reflection coefficient of 
this reference as seen from the generator side, then 


(a) Reference Step and Typical Disturbance 



Ratio 
N : 1 


(b) Transformer Model 


7ig ure 11. A Model for the E-Plane Taper 
and a Typical Disturbance 
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r r = s 2 2 = <N 2 -1)/(N 2 +1) (4) 

and S 33 = (1-N 2 )/(1+N 2 ) (5) 

and S 32 = S 23 = 2N/(1+N 2 ) (6) 

Since r r has been chosen to have zero phase it 
will be written without a modulus sigij,^ 

The total reflection coefficient at the reference 

( 

plane is found by substituting Equations (4), (5), (6) 

into Equation (2 ) , the result being 

r = bi/a! = (r r (i+r r ri) +s 23 s 32 r 1 )/ (i+r r ri) (7) 

but S 32 • S 23 = l-£ r 

therefore r = (T r + TiJ/d+TyTi) 

= (r r +|r i |exp{j(0-2$L)}j/^l+r r |r 1 |exp{j(0-26L)}) (8 ) 

A measured curve of the total reflection coefficient 
T is shown in Figure 12 for a fabricated capacitive iris and 
the fabricated reference step. The oscillations in the 
waveform as a function of frequency is in part a function of 
separation distance L as described by Equation (8) . 




FREQUENCY IN QHZ 


Figure 12. Total Magnitude of Reflection Coefficient 
for a Capacitive Iris 27.0 Centimeters 
from the Reference Step 


Let c = 1 | cos O- 23 L) 

s = | r x | sin (6-2$L) 

Substituting Equations (9), (10) into (8), 

r = [ (r r +c+js) ]/[i+r r (c+js) ] 


r=[r r (i+|r 1 | 2 )+c(i+rrj+js(i-r r 2 )]/[i+2r c+r 2 | r x | 2 ] 
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Part of the return wave passes through the directional 

coupler C4, Figure 6, and produces a DC voltage which is 

independent of the phase of T. The output of the crystal 

2 

detector at C3 is proportional to | r | • Therefore, squar- 

ing the magnitude of (12) , 

I r | 2 = [r r 2 d+|r 1 | V + 2 cr r (i+| t 1 j 2 ) (i+r r 2 ) 

+c 2 (i+r r 2 ) 2 + s 2 (i-r r 2 ) 2 ]/(i+2cr r 

+r r 2 l r i| 2 ) 2 • (13) 

We are concerned chiefly with small reflections, 

3 3 

therefore terms such as T r | r x J are small compared to unity. 
If T r = .2 and | r ^ | = .1 then r r 3 1 r ^ | 3 is equal to 8 x 10 6 . 
Ignoring these and terms containing cos (20 -4 0L) for later 
simplicity, the denominator of Equation (13) may be written 
as : 

(i+2cr r +r r 2 |r 1 | 2 ) 2 

~ i/(i-4cr r +4r r 2 |r 1 | 2 ) , (14) 

with these approximations. 
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<i-rr 2 -|r 1 | 2 +r r 2 |r 1 | 2 )cos(e-26L) 


(15) 


In order to compute JrjJ cos (0-2BL) , (15) may be 

rearranged. 


| T 1 1 cos (0-2BL) 


2 2 2 2 2 

J r 1 - r r .+ ^,1 - 2 r r ;|r 1 l 

2 r r < 1 - r r 2 -l r il 2 +r r 2 l r i* 2) 


(16) 


Let . A = (i-r r 2 -|r 1 | 2 +r 2 |r 1 | 2 ) , (17) 

and G (v) = |r| 2 /2r r A. (18) 

Then (16) can be rewritten as, 

| r ! | cos (0-2 Bl) * G(v)-(r r 2 +|r 1 ] 2 -2r r 2 |r x | 2 )/2r r A ( 19 ) 


The term G(v) is a function of wavenumber and is to be 
determined from the readings of the digital volt meter as 
will be described in the next sections. A is a correction 
term near unity. 

Up to this point in the derivation only one reflec- 
tion has been considered in addition to the reference. In 
general, the components being measured will contain more than 
one reflection, , r 3 • • .r n , at distances 1^, L 2 L , and 
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these must be computed from the data G(v). 

Assume that interactions between adjacent test re 
flections may be neglected, then by superposition, 


G(v) 



r n | cos (6 n -4iTL n v) 


The above equation may be written as a Fourier 

series. 


G(v) 


-z 

r»=l 


[a n cos (47rnL^v)+b n sin (47rnl/v) ] 


where L n = nL' and h' is the unit length. 

The component reflections are found by taking the 
finite Fourier transforms. 


v. 


a n = 1/v' J G ( vj cos (nirv/v ^ ) dv 

v, 


and 


b n = 1/v 


v 2 

G (v) sin (n7rv/v^) dv 


where 2v" = v 2 - v lf the range of wavenumber interval in 

* 

which the G(v) is measured. And L ' = ^ 

4 

From (20) a n and b n are also equal to 


( 20 ) 


( 21 ) 


( 22 ) 


( 23 ) 
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£tn — | r n | cos0 n (24) 

b n = l r nl sin6 n (25) 

using the trig expansion, 

cos (x - y) = cos x cos y + sin x sin y . 

2 2 h 

Therefore, |r n | = (a n +bn ) 

f 

and 0 n = arctan (b n /a n ) . 

The phase angle of the total voltage reflection co- 
efficient was lost when the r.f. signals were rectified by the 
detector crystals, however, both the magnitude and phase 
angle of the measured component may be computed from the 
transforms. It is well to point out that |r| is the mag- 
nitude of the total voltage reflection coefficient when the 
reference step and the unknown device to be measured are 
attached in the guide. 

V 

Substituting L' - j into the sine and cosine argument 
of Equations (22) and (23) , 

mrv/v' = 47TL n V. (28) 


(26) 

(27) 
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Since L n = n/2(v 2 -vi) , (29) 

for an arbitrary distance L the wavenumber range v 2 - Vj. 
must be adjusted. The message here is that for an arbitrary 
distance along the waveguide, the angle traversed by the ar- 
guments of the sine and cosine terms of the Fourier transforms 
(22), (23) must be an integral multiple of 2 tt. This is im- 
portant to recognize since the computer program will be cal- 
culating the Fourier integrals for an arbitrary distance L. 
Therefore the limits of the wavenumber range must be adjusted 
to satisfy Equation (22) and (23) . 

Method of Measuring G(v) . The Hewlett-Packard Model 
416B ratio meter has been used, Figure 6. This particular 
ratio meter requires inputs of 1 KHz signals from two square 
law detectors. Two sets of measurements are taken, one with 
T ^ measured and a second with r 2 measured. The corresponding 
voltages out and E 2 , obey the equation, 

2 2 

|T 1 | /\V 2 \ = tan(CE 1 )/tan(CE 2 ) , (30) 

where C is a constant peculiar to the ratio meter used and 

can be related to a conduction angle, i.e., time lag between 

pulses generated in the ratio meter. Consider now that if we 

take one set of readings with the reference step terminated 

in a matched load, then I r I = r Then consider that the 

1 2 1 r 
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second set of measurements are taken with the reference step 
terminated by the unknown device to be measured. Using this 
approach Equation (30) may be rewritten as: 


2 2 

| r I = T tan (CE)/tan (CE ). (31) 

r r 


The quantities T r and C are known from, previous 
measurements. The quantities E and E r are measured for each 
unknown. Using Equation (31)/ Equations (22) and (23) may 
be rewritten as 


v- 


a = T /2v'h I [tan (CE)/tan (CE ) ] cos (2 ttLv) dv 


n 


(32) 


and b = T /2v'A 
n r' 


v 2 

[tan (CE)/tan (CE r ) ] sin (2 ttLv) dv. (33) 


Equations (32). and (33) are the equations which are 
used by the digital computer to calculate the reflection co- 
efficient as a function of distance. 

At a particular distance corresponding to the loca- 
tion of a discontinuity in the guide, the complex reflection 
coefficient as a function of frequency is determined by in- 
tegrating Equation (32) and (33) over small overlapping 
wavenumber ranges and selecting the frequency corresponding 
to the wavenumber in the center of each range as the fre- 
quency of the particular integration. Figure 13. Each sub- 
interval initially has ten data points or 500 MHz bandwidth 


35 


and is staggered two data points or 100 MHz. Each subin- 
terval must satisfy Equation (29) for the location down the 
waveguide. Since this location is an arbitrary number, the 
wavenumber range must be appropriately adjusted such that 
an even interger multiple of it i a equal to L n , the distance 
location. This adjustment is done from the low side of the 
subinterval. 


11.875 I 
11.925 


11.975 


12.975 112.175 


1l2. 


275 


12.025 12.125 12.225 12.325 


11.975 

12.025 


j 1 1 1 1 r- 

11. 075 1 11. 175 |ll. 275 jll.375 | 
12.125 12.225 12.325 12.425 


Figure 13. Wavenumber Subintervals Given in GHz 

As a result of the overlapping procedure and the subinter- 
val adjustment, the calculated results of complex reflection 
coefficient* are known at unequal frequency intervals of 
approximately 100 MHz. While the ref lectometer data is 
taken at equal frequency intervals, the integrals are in 
terms of wavenumber, 1/Ag,** which does not correspond- 
ingly occur at equal intervals. The method used in pro- 


*At this point the complex voltage reflection co- 
efficient F is transformed to the complex current reflection 
coefficient r j by the equation r j = -T .' 

**The usual definition, 27r/\g is not used here for 
convenience. 
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gramming the integrals (32) and (33) is a point by point 
calculation of the area under the curve formed by the argu- 
ment of the integrals. This is discussed further in 
Appendix B . 


Selected Waveguide Discontinuities 

Introduction 

Four waveguide discontinuities were selected as a 
library of matching elements. Families of curves of the 
characteristics of these elements, the symmetrical capaci- 
tive iris, the asymmetrical inductive iris, the solid metal 
inductive post, and the dielectric post, are given in Appen- 
dix A. The above elements also serve as matching elements 
for the impedance matching computer program. The equations 
which model the matching elements are programmed in the match- 
ing program and are called upon to impedance match the dis- 
continuities measured by the comparison ref lectometer . 
Symmetrical Capacitive Iris 

The symmetrical capacitive iris shown in Figure 14a 
is modeled by the circuit in Figure 14b. This iris is 
formed by obstacles of small but finite thickness with edges 
perpendicular to the electric field (H^q - mode in rectangu- 
lar guide). The equations used are given by N. Marcuvitz 
and are in terms of the physical dimensions of the iris. 

Figure 14, on the next page, shows this arrangement. 


Cross Section View 


Side View 


(a) Symmetrical Capacitive Iris 


-jB 



(b) Equivalent Circuit 


Figure 14. Symmetrical Capacitive Iris and Equivalent Circuit 
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N 


Consider now the equivalent circuit of Figure 14b. 
is given by: 


V Y q = |csc(2ttL/X ) 


(34) 


and 


where 


V Y o = V Y o + 3 tan (iTL/X g ) 

B/Y = 2b[ In secfl^.TTdVL 
1 o A_ V 2b / 2bd 


(35) 


2 


4 I TTd'\ z 

+Asin \ 2b 1 / b \ / T ^ 2 / Trd'\^\ . 4 / Trd'\ 

l + Aoos 4 16 ( V ( M COS ( 2b ) H ( 2b ) 


(36) 


and 



(37) 


Where the physical parameters b, d' , L, and d are shown in 

Figure 14a, aind A is the guide wavelength. 

y 

Equations (34) to (37) are programmed in the impe- 
dance matching program. They were also programmed in order 
to calculate the library of families of curves in Appendix A. 

Restrictions . The equivalent circuit is valid in the 
wavelength range b/A^ < 1, where Equation (34) is in error by 
less than 2 percent. Equation (36) is in error by less than 
5 percent when d'/b < 0.5 and 1/d <0.5. A comparison of 
calculated and measured results are made in Chapter V. 
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Asymmetrical Inductive Iris 

The asymmetrical inductive iris shown in Figure 15a 
is modeled by the circuit in Figure 15b. This inductive 
iris is formed by an obstacle of small but finite thickness 
with edges parallel to the electric field (H^q - mode in 
rectangular guide),. The equations used are given by N. 
Marcuvitz and are in terms of the physical dimensions of 
the inductive iris. 

Consider now the equivalent circuit of Figure 14b. 

Where: 



V z o - 

(4a/X g ) (a/7rD') 2 , irD'/X <<1, 

(38) 

and 

V Z o - 

(a/16X g ) UD^a) 2 , «1 

(39) 

where 

D' » * 

J 2 

[ 1+ H’ In ( T )]> 

(40) 


and D s (4Ld‘ 3 /3 tt) ^ , L/d' << 1. (41) 

Where the physical parameters b, d* , L and a are shown in 
Figure 15a. The free space wavelength is X and the guide 
width is a. Equations (38) through (41) are programmed as 
part of the impedance matching program. They were also 
programmed in order to calculate the library of curves in 
Appendix A. 




Cross Sectional 
View 


T T 

Top View 


(a) Asymmetrical Inductive Iris 


-jx b 



(b) Equivalent Circuit 


Figure 15. Asymmetrical Inductive Iris and Equivalent 
Circuit 
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Restrictions . The equivalent circuit is valid in the 
wavelength range a < X < 2a. For D' , < 0.2a and a < X 

Equation (38) through (40.) are estimated to be in error by 
less than 10 percent. A comparison of calculated and mea- 
sured results are made in Chapter V. 

Solid Inductive Post 

A solid metallic post of circular cross section with 
axis parallel to the electric field (H^q “ mode in rectan- 
gular guide) is shown in Figure 16a with the equivalent cir- 
cuit in Figure 16b. The equations used are given by N. 
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Marcuvitz and are in terms of the physical parameters of 
the inductive post. 

Consider now the equivalent circuit of Figure 16b. 

Where: 



X a /z 0 - X b /2Z 0 = 

■ ( 

(£) 

’ 

-/^\ 2 (S n cot 1 

( — ) 

-S.) 2 1 esc 2 ( 

7TX \ 


\2a / ° 

\ a / 

1 -> \ 

a / 

and 

x b - a_ j 
z o 

^ TTd 

l a , 

fsi„ 2 ( F ) , 



where 


S o “ ln 


4a sinux 
Lird a. 


2 sin 



(42) 


( 43 ) 


Cross Sectional View 


Top View 


Figure 16. 


(a) Solid Metal Inductive Post 



(b) Equivalent Circuit Parameters 


Solid Metal Inductive Post and Equivalent Circuit 
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The physical parameters x, d, and a are shown in Figure 16a. 
Equations (42) through (45) are programmed as part of the 
impedance matching program. They were also programmed in 
order to calculate the library consisting of families of 
curves in Appendix A. 

Restrictions . The equivalent circuit is applicable 
in the wavelength range a < X < 2a. The results are reliable 
to within a few percent for d/a < 0.1 and 0.8 > x/a > 0.2. 

A comparison of calculated and measured results will be 
made in Chapter V. 

Dielectric Post 

A cylindrical dielectric obstacle of circular cross 

section aligned parallel to the electric field, (H 1Q - mode 

in rectangular guide) , is shown in Figure 16a with the 

equivalent circuit in Figure 17b. The equations used are 

13 

given by N. Marcuvitz and are in terms of the physical 
parameters of the inductive post. 



Cross Sectional View 


Top View 


Figure 17. 


(a) Dielectric Post 



(b) Equivalent Circuit 


Dielectric Post and Equivalent Circuit 
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Consider now the equivalent circuit of Figure 17b. 


2 z 

(2a/Ag) (frd/a) sin (irx/a) 

g 2 Ji (3) 1 

Jl (a) aJ o (a)J 1 (0)-6J o (6)J 1 (a) 


Jq(B) and J^(B) are the bessel functions of the first kind 


with orders 0 and 1 respectively. 


= iTd/X , 3 = Je' ird/X. 


X a /Z Q - X b /Z Q = (a/2X g )csc 2 (Trx/a)[j o (0)/J o (o) 


*3J 0 ) J i (3) ~ gJo^3)Ji(g) 


- So + a/4 


where 


S Q = In ud 


sin MI-2 sin 2 M 


L f 1 1 • _ 2 / nrrx\ 

/n 2 -(25) 2 ’ "J ' a ' 

n=2 aJ \ a / 


where the dielectric constant e 1 = — , and e Q is the free 

space dielectric constant. For an obstacle with a complex 

dielectric constant, e'-je"= the above formulas are still 

e o 

valid provided e', X a /Z Q and X b /Z Q are replaced by e'-je" and 
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j (X a /Z Q ) , and ~j (Xj^/Zq) respectively. It should be noted 
that resonant effects occur for large values of e' with 
attendant changes in sign of the circuit elements. Only 
low loss dielectrics are considered herein and only the 
real part of the dielectric constant will be considered in 
this research. 

Restrictions . The equivalent circuit is applicable 
in the wavelength range 2a > A > a, and for the centered 
cylinder (x = a/2) in a wider range 2a > X > 2a/3. Equa- 
tions (46) through (48) are estimated to be in error by 
only a few percent in the range d/a < 0.15 and 0.2 < x/a < 
0.8 provided that neither X a /Z Q nor K^/Z Q are too close to 
resonance. Since this research is concerned with reflec- 
tion coefficient values <_ . 2 the above resonance criterion 
presents no problem. A comparison of calculated and mea- 
sured results are made in Chapter V. 

Matching Computer Program 

Introduction 

The computer program which calculates the physical 
dimensions of the matching element and its location from 
the mismatch is referred to here as the matching program. 
This matching program uses an iterating technique to deter- 
mine the desired physical dimensions. The criterion which 
has been used is the least root mean square of the resultant 
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mismatch (LRMS).* The RMS of the mismatch reflection co- 
efficient has been used as the matching criterion in this 
research to reduce the net energy reflected over the band- 
width. This criterion does not always lead to the lowest 
peak VSWR over the matched band. If it is important to 
decrease the upper limit of the VSWR obtained by using the 
RMS criterion then the criterion described in the recommen- 
dations section of Chapter VI may be used at the possible 
expense of increasing the average VSWR over the matched 
bandwidth. The term resultant mismatch is used here to 
identify the reflection coefficient observed by the gen- 
erator when a matching element is present in the waveguide . 

In order to assist in the selection of a matching 
element, a library of selected matching elements has been 
compiled and is given in Appendix A. The selected matching 
elements consist of the symmetrical capacitive iris, the 
asymmetrical inductive iris, the solid metal inductive post, 
and the dielectric post. Information is presented in the 
form of plots of magnitude of current reflection coefficient 
and phase angle as a function of frequency. Also, plots of 
shunt susceptance are given, normalized to the characteris- 
tic impedance, as a function of frequency. 


mismatch reflection coefficient) 


*RMS = ± l , 

M n= 1 

where M is the number of data points in the preselected 
bandwidth. 
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Program Description 

The impedance matching computer program is described 
with reference to the functional flow chart of Figure 18. 

A listing of the program statements is given in Appendix B. 
The impedance matching computer program is an iterative pro- 
gram which iterates in three dimensional space. The space 
has three coordinates of which two are of the physical di- 
mensions of the matching element. The third dimension is 
the location of the matching element relative to the mis- 
match. A 3 x 3 x 3 array is used in dimension space form- 
ing a total of twenty-seven points. 

The program begins by reading in data and storing 
needed constants. Nine sets of data cards are required as 
input for this program. The symbol corresponding to the 
particular data read is presented in Table 1. An explana- 
tion of the data cards read is also given by comment state- 
ments in the matching program listed in Appendix B. 

Depending on the matching element chosen, an ini- 
tial center point in dimension space is selected. This 
point is denoted as the (1,1,1) point. If the capacitive 
or inductive iris is selected, the I coordinate corresponds 
to the thickness of the iris, and the J coordinate corres- 
ponds to the height of the iris. If, however, the metal 
inductive post or the dielectric post is selected as match- 
ing elements , then the I coordinate corresponds to the 
distance from the centerline of the post to the sidewall. 
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Figure 18. Impedance Matching Program Functional 
Flow Chart 
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Table 1. Impedance Matching Computer Input Data 


Symbol 

MACHEL 

EPRIM 

LE 

FR(I) 

XEEL 

YIMAG 

FRBEG, 

FREND 


Corresponding Data 

Designates matching element, 

Value of the dielectric constant of the 
dielectric post matching element 

Number of data points of the calculated 
mismatch reflection coefficient 

Frequency in GHz corresponding to data 
points of the calculated mismatch reflec- 
tion coefficient 

Real part of the mismatched reflection 
coefficient to be impedance matched 

Imaginary part of the mismatched reflec- 
tion coefficient to be impedance matched 

Band limits .to be matched over 


NOO Number of times mismatch is to be impe- 

dance matched 

ALENGH(I) Value of the initial relative distance 

between the measured mismatch and the 
preselected matching element (the number 
of different ALENGH (I ) equal NOO) 
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and the J coordinate corresponds to the diameter of the 
post. In all cases the K coordinate corresponds to the 
distance between the matching element and the mismatch. In 
step 2 of Figure 18, the (1,1,1) point is used to generate 
a set of satellite points on either side of the (1,1,1) point 
lying on the three axes. These points are referred to as the 
nearest neighbor points and are described by the coordinates 

(1.1.1) , (3,1,1), (1,2,1), (1,1,3), (2,1,1), (1,3,1) and 

(1.1.2) , as shown in Figure 19. The nearest neighbor points 
are of the greatest interest, therefore, in order to conserve 
computer run time calculations are made only for nearest 
neighbor points. 

Step 2 of Figure 18, therefore consists of incre- 
menting the (1,1,1) point to generate the 7 principal points 
of the 27 point array. In step 3 the matching element is 
selected according to the information given by the first 
data card. Table 2 shows the initial fconditions used in the 
matching program.* In step 4 the reflection coefficient of 
the matching element chosen is calculated over the desired 
bandwidth (given by the last data card) . This is done for 


*The initial conditions were selected to reduce the 
run time of the matching program by selecting them as rea- 
sonably as possible. Large deviations of the physical di- 
mensions of the matching elements were made and in all cases 
the resulting dimensions returned to the same value. There- 
fore, there was an absence of secondary minima as a function 
of matching element physical dimensions. 



( 1 , 1 , 1 ) ' ( 3 , 1 , 1 ) 


( 1 , 2 , 1 ) 


Space Coordinates , Showing 
Nearest Neighbor Points 
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all seven points in dimension space. Step 5 calculates the 
resultant mismatch between the disturbance in the waveguide 
(given by data) and the calculated matching element for all 
seven nearest neighbor points in dimension space. The equa- 
tion used is referred to as the matching equation derived in 
Appendix C and given here as: 


n Ia exp t-j26L]+r Ib -2 r Ia r Ib exp t-j23L] 

n ti -n r -1 o n t 1 


r ia is the current reflection coefficient of the disturbance 
at point A in an otherwise matched transmission line, r ib 
is the current reflection coefficient of the disturbance at 
point B in an otherwise matched line, and L is the distance 
between points A and B in Figure 20.* Equation (49) gives 
the total current reflection coefficient as seen by the gen- 
erator as a result of two disturbances located in the line at 
points A and B. The equation is in. terms of the individual 
current reflection coefficients that would be generated if 
each disturbance were placed individually in an otherwise 
matched transmission line. The program is implemented such 
that either r Ia or r Ib can play the roll of the matching 
element. If the disturbance is matched by placing the 

*The complex current reflection coefficient Tj is 
equal to the negative of the complex voltage reflection co- 
efficient T, i.e. Tj = -r. 
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matching element on the generator side, then the matching 
element is represented by If matching is accomplished 

on the load side of the disturbance, then the matching ele- 
ment is represented by Tj a . 

Table 2. Matching Element Initial Conditions on Dimensions 


Element 

Thickness 

Height 

Capacitive Iris 

0.05 cm 

0.32 cm 

Inductive Iris 

0.05 cm 

0.32 cm 

Element 

Diameter 

Sidewall Distance 

Metal Inductive Post 

0.16 cm 

0.47 cm 

Dielectric Post 

0.16 cm 

0.47 cm 


At this point in the program the resultant mismatch 
as a function of frequency at intervals of approximately 100 
MHz has been calculated for each point in dimension space. 
Step 6 of Figure 18 determines the root mean square (RMS) 
value of the current reflection coefficient over the fre- 
quency band of interest, i.e., the seven RMS mismatch values 
corresponding to the seven points in dimension space. Step 
7 of Figure 18, selects the smallest of these RMS values 
or least root mean square (LRMS) and determines its corres- 
ponding coordinates in dimension space. 
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Step 8 is a branch decision. If the point selected 
in step 7 is not equal to the (1,1,1) point then the branch 
is made to step 11 and the values of the dimensions of the 
point (I,J,K) corresponding to the LRMS are assigned to the 
(1,1,1) point. In effect this moves the (1,1,1) point in 
dimension space. The computer returns to step 2 and calcu- 
lates new nearest neighbor points around the new (1,1,1) 
point. If, however, the point corresponding to the LRMS in 
step 8 is in fact the (1,1,1) point, then a regional minimum 
has been determined and a check is made on the size of the 
increment. If the increment size has not been reduced to 
1/8 the original size at the beginning of the program the 
computer reduces the existing increment size by a factor of 
2 in step 12 and returns to step 2 to determine a better 
representation of the minimum mismatch. On the other hand, 
if the increment size has been reduced by a factor of 8 then 
the computer exits step 9 and prints, plots and punches re- 
sults. The initial increment size used in the oroaram 'is 
shown in Table 3 and was determined empirically to avoid 
large differences in reflection coefficients from neighbor- 
ing points . 

Program Output . The matching program output con- 
sists of printout, plots, and punched cards. A typical 
printout is given in Figure 21. The type of matching ele- 
ment, the bandwidth matched over, the physical dimensions 
of the matching element and the location of the matching 


CAPACITIVE IRIS IS THE MATCHING ELEMENT 


NUM3tR OriNCREMENT = 8. 


MATCHED BANG IS FROM 9.795 GHZ* T O 10*99 5 GHZ* __ 

“TffIC'<NtSS= .0438 IN CM HEIGHTS 71600 IN CM &T5fANCE TOWARD THE GEn 7“FR0M Ml §MAT£ h= 73375 CM 


RESULTANT MISMATCH 


— MATCHED 


UNMATCHED ~ 


FREQJENCV IN GHZ. 

MAG. REFL. COEFfT 

vswft 

mag.reFl.Coe 

:fF.DIS~urb. measured 


8.195 

.07719 

1.1673 

•06696 



0.395 

.07506 

1.1623 

.07138 



8.595 

•07186 

1.1549 

.07563 



0.795 

.06767 

1.1452 

.07973 



8.995 

.06255 

1.1334 

.08373 



9.195 

.05655 

1.1199 

.08761 



9.395 

.04973 

1.1047 

.09140 



9.595 

•04214 

1.0880 

.09511 



9.795 

•03381 

1.0700 

.09875 



9.995 

.02480 

1.0509 

.10232 



10.195 

•01516 

1.0308 

.10504 



10.395 

.00510 

1.0103 ' 

.10931 



10.595 

•00638 

1.0128 

.11272 



10.795 

.01758 

1.0358 

•11609 



10.995 

♦02941 

1.0606 

•11942 



11.195 

.04171 

1.0670 

.12271 



11.395 

•05440 

1.1150 

.12596 



11.595 

.06743 

1.1446' 

.12919 



11.795 

•08076 

1.1757 

•13238 



11.995 

•09434 

1.2083 

.13554 



12.195 

•10812 

1.2424 

•13868 



12.395 

•12205 

1.2780 

.14179 




Figure 21. Printout of Matching Program Results 
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element from the mismatch toward the generator is given. 

The modulus of the reflection coefficient of both matched and 
unmatched discontinuity as well as the VSWR of the matched 
discontinuity are given. 

A typical plot is given in Figure 22. The VSWR of the 
original discontinuity as well as the VSWR of the impedance 
matched discontinuity are given on the same plot as a 
function of frequency. 

In order to add more flexibility to the impedance 
matching computer program, the output also consists of a set 
of punched cards which are punched according to the format 
of the data cards for the matching program. The real and 
imaginary parts of the resultant mismatch are punched on 
these cards. Therefore, if it is desirable to improve on 
the existing impedance match, the output cards may be used 
as data cards for the matching program and the type of 
matching element changed to provide, in many cases, an 
improved impedance match. The result is the original 
discontinuity matched by two matching elements. The lo- 
cation of the second matching element calculated by the 
matching program is referenced to the element nearest the 
generator, resulting from the first impedance match. This 
can be done as many times as may be deemed practical. 
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Table 3. Initial Increment Size of Physical Dimensions 


Dimensions* 

Capacitive and 
Inductive Iris 

Dimensions* 

Metal and 
Dielectric Post 

Relative 

0.1 

Relative 

0.1 

Distance 


Distance 


Thickness 

0.01 

Sidewall 

0.01 



Distance 


Height 

0.1 

Diameter 

0.05 


*in centimeters 



Figure 22. Typical Example of the Matching Program Output 
Plot of the VSWR of an Unmatched and Matched 
Disturbance 



CHAPTER IV 

EXPERIMENTAL EQUIPMENT 
Introduction 

The comparison ref lectometer system and various test 
waveguide discontinuities constitute the experimental equip- 
ment. The comparison ref lectometer shown in Figure 23 is 
comprised of various microwave components including a back- 
ward wave oscillator (BWO) , a coherent synchronizer, a 

ratio meter, and directional couplers. 

Four waveguide discontinuities were fabricated and 
measured to demonstrate the validity of Equations (III-34) 
through (III-48) . A shunt slot radiator was also fabricated 
to demonstrate the validity of the impedance matching tech- 
nique. 


Comparison Ref lectometer System 
Functional Diagram . A block diagram of the comparison 
ref lectometer system fabricated for this research is shown in 
Figure 23. A photograph of the system is given in Figure 7. 

Microwave energy is generated by the backward wave 
oscillator (BWO) and is directed through the waveguide to 
the pin diode modulator. The modulator amplitude modulates 
the microwave energy at the audio oscillator frequency of 
1 KHz. While the audio oscillator generates a sine wave. 


Scope to Monitor 
Phase Locks 


Punch Signal 


Push Button 
to take Reading 



\ 

\ 

\ 



Figure 23. Comparison Ref lectometer Block Diagram 
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the modulator is driven to the extent that the output of the 
modulator is an amplitude modulated square wave. The iso- 
lator serves to prevent reflected energy from pulling the 
BWO frequency. The wave is next partially reflected by the 
reference step. This step reflects a known ratio of inci- 
dent energy from the reference plane. The amplitude and 
phase of this reflection coefficient is constant over the 
frequency band. Next either the matched load is attached, 
in which case essentially all of the remaining wave is ab- 
sorbed, or the unknown is attached in which case part of the 
energy is reflected from the unknown and the remaining en- 
ergy is absorbed in a matched load. 

The two directional couplers Cl and C2 provide feed- 
back necessary to ensure a phase locked, constant power out- 
put from the BWO. Couplers C3 and C4 are the essential el- 
ements of the ref lectometer . The incident wave traveling 
toward the load is coupled out by C3. The reflected wave is 
coupled out by C4. Couplers C3 and C4 are terminated in a 
square law detector. The output of each detector is a 
square wave whose amplitude is proportional to the square of 
the electric field intensity incident on the respective 
detector. The ratio meter takes the ratio of these two in- 
put signals and provides a voltage output that is propor- 
tional to the square of the total reflection coefficient. 
This reflection coefficient can either be the total re- 
flection coefficient of the reference step and the measured 


63 


device or it can be that of the reference step terminated 
by a matched load. This voltage output E is measured by the 
digital voltmeter, and the value measured is punched on 
paper tape at the command of the push button. 

The BWO is phase locked every 50 MHz by the coher- 
ent synchronizer from 7.975 GHz to 12.425 GHz. At each of 
these 89 phase lock points a measurement is taken by de- 
pressing the push button. Figure 23. One complete set of 
measurements is made with the reference step terminated in 
a matched load and a second complete set is taken with the 
reference step terminated by the unknown device. During 
these measurements the gain of the ratio meter is left un- 
disturbed. From the data collected the magnitude of the 
total reflection coefficient of the reference step termina- 
ted in the unknown device may be calculated by the digital 
computer. Detailed operating instructions for the comparison 
ref lectometer are given in Appendix C. 

Backward Wave Oscillator . The RF source used in 
this research is an Airborne Instruments Laboratory, Sweep 
Oscillator type 210. This BWO is used with the plug in 
unit which provides a frequency sweep from 7.975 GHz to 
12.425 GHz. The oscillator frequency can be controlled 
by an external sweep voltage or by an internal sweep gen- 
erator. The internal sweep can be set to a repetitive 
sweep with a variety of sweep rates or can be set up to 
sweep once through the band and return to the lowest 



64 


frequency, remaining there until activated by a front 
panel push button. The rate of this sweep is an adjustable 
parameter. In the reflectometer setup of Figure 23 the BWO 
is operated in this internal signal sweep mode. The RF 
signal is then blanked out while the oscillator is return- 
ing to the lowest frequency. The oscillator then remains 
at this starting frequency until again triggered by the 
front panel push button. In order to enable the coherent 
synchronizer to phase lock and a voltage reading to be 
taken, the oscillator sweep is slowed down by adding an 
external 200 micro farad capacitor to the connector pro- 
vided on the back panel of the unit. 

When the oscillator frequency is in a capture re- 
gion of the coherent synchronizer a correct error signal to 
the helix of the BWO will result in the phase locking of 
the BWO to the crystal controlled oscillator of the co- 
herent synchronizer. The frequency is then held within 7 
parts in 10 . When the internal sweep voltage increases 
sufficiently to overcome the error correction signal on the 
helix of the BWO, the phase lock is broken. The error sig- 
nal returns to zero, and the sweep signal resumes control of 
the RF output. The RF frequency, therefore, continues to 
increase until it again is within the lower bound of the 
2.5 MHz capture region of the 25 MHz IE of the coherent 
synchronizer. In this way, the BWO is phase locked every 
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50 MHz from 7.975 to 12.425 GHz. This is done automatically 
with perfect repeatability when the system is properly set 
up. The entire run of 89 measurements requires 10 to 15 
minutes depending on the selected sweep rate. 

Coherent Synchronizer . The instrument used to obtain 
desired phase locks is the Sage Laboratories, Inc. model 244 
coherent synchronizer. This device is generally used to 
phase lock a source at a single frequency. Because of this, 
the unit is equipped with five crystals in the basic os- 
cillator. Depending on the frequency at which the phase 
lock is desired, one of the crystals operating at 5000.000 
KHz, 5012.500 KHz, '5006. 250 KHz, 5000.00 KHz, or 4993 . 750 KHz 
is selected. One of the 5000.000 KHz crystals is used in 
the fixed mode. The other four crystals are used in the 
variable mode. In the variable mode,, by selecting the 
correct crystal and pulling the frequency of this crystal by 
a variable front panel control, any frequency from 100 MHz 
to 18 GHz can be phase locked. However, in the fixed mode 
the basic oscillator operates at a crystal controlled fre- 
quency of 5.000000 MHz. This signal is multiplied by a 
frequency multiplier chain to provide a 100 MHz signal. 
Figure 24. Functionally, the 100 MHz signal excites a 
varactor multiplier which in turn generates harmonics from 
100 MHz to 18 GHz. Each harmonic is 100 MHz apart. This 
signal is then applied to the mixer which also has as an 
input the RF signal from the BWO oscillator to be phase 


5MHz 
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Figure 24. Simplified Block Diagram of Coherent Synchronizer 







67 


locked. The output of the mixer is the sum and difference 
of the two inputs. This output is the input to the 25 MHz 
amplifier which has a 2.5 MHz 3 dB bandwidth. If there 
exists a signal within approximately 1.5 MHz of 25 MHz 
this signal is amplified. Other signals are not amplified. 
The level meter gives an indication of the presence of a 
25 MHz signal. Therefore, when a phase lock has occurred 
there is positive indication on the level meter. The output 
of the 25 MHz amplifier drives a phase detector which com- 
pares the phase of the 25 MHz reference signal tapped from 
the multiplier chain with the output of the 25 MHz amplifier. 
The output of the phase detector is an error signal which, 
after being amplified, is applied to the helix of the back- 
ward wave oscillator. This error signal can be monitored 
by the phase meter. 

Because the comoarison ref lectometer phase locks 
at 50 MHz intervals, the reference oscillator of the coher- 
ent synchronizer is set in the fixed position. Therefore, 
the output of the harmonic generator is a frequency comb 
with harmonics at intervals of 100 MHz from 100 MHz to 
18 GHz. A phase lock will occur when the RF signal input 
to the mixer is 25 MHz on either side of a harmonic on the 
comb, Figure 25, since the difference of the harmonic and 
the input RF from the BWO is 25 MHz at these points. 
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Figure 25. Partial Frequency Comb of Harmonic Generator 


It was particularly important for this research to 
optimize conditions of the harmonic comb of the harmonic 
generator. In order to make a rapid set of measurements it 
is important to phase lock automatically from one desired 
frequency to the next. This requires essentially a single 
adjustment to the coherent synchronizer controls. It was 
discovered initially that this was impossible to do and in 
fact some of the desired phase lock points were not achiev- 
able. Further investigation pointed out that the spectrum 
of the harmonic comb from the harmonic generator was not at 
a constant level. In fact, some harmonics were apparently 
not present. 

While it is convenient to consider the harmonic gen- 
eration and the mixing as two separate operations , in the 
model 244 both are accomplished in one step by a single 
varactor. It was discovered that this varactor was not bi- 
ased correctly to provide a relatively constant 25 MHz IF 
signal when the BWO was swept through the frequency band. 
This is equivalent to saying that the harmonic generator 
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was not providing a flat spectrum throughout the bandwidth. 
Adjustments were made to achieve this flat spectral behav- 
ior by observing the 25 MHz signal at the BNC connector 
"TP1" provided on the 25 MHz amplifier. Details of this 
adjustment are given in Appendix D. 

A single adjustment to the error signal potentio- 
meter enables consistent automatic phase locks as the BWO 
sweeps slowly through the band. Because the IF levels are 
relatively constant the time duration of the phase locks 
are also relatively constant. Further, this phase lock 
time can be adjusted by increasing or decreasing this max- 
imum allowable error to the helix of the BWO. By increas- 
ing the maximum error voltage a given phase lock exists for 
a longer period of time. 

Ratio meter . The ratio meter, Hewlett-Packard model 
416B was used to determine the ratio of the incident and re- 
flected electric field intensities. Figure 23. This ratio 
meter has two inputs, both are 1 KHz square wave signals 
whose amplitudes are proportional to the square of the inci- 
dent and reflected electric field intensities of the output 
of the 10 dB directional couplers. Two Hewlett-Packard model 
424A crystal detectors operating in the square law region 
were used, at the output of the two directional couplers C3 
and C4, Figure 23. 

The output voltage of the ratio meter is related to 
the reflection coefficient measured by the relation: ^ 
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2 

| T | = £ tan 0 , (1) 

where ^ is a proportionality constant depending on the ref- 
erence adjustment, and the angle 0 is expressed in terms 

12 

of the output voltage measured: 


0 = yE/E 


max ' 


( 2 ) 


where E is the voltage measured, E is the voltage meag- 
ured when the meter reading is at the uppermost point of 
the scale, and y is defined as a conduction angle. The con- 
duction angle Y was determined by programming Equation (3) 
on a desk calculator. A discontinuity with a known value 
of | r ( f ) | was selected and measured at two different fre- 
quencies such that, 


lilij = tan 

l r 2 l tan tyE,/^} 


(3) 


E 1# E 2 and Ejjj^ were measured and | F -l | and | T 2 1 were known. 
The unknown y was calculated by an iterative technique to 
satisfy Equation (3) . In most cases the value of E m , v for 
the particular ratio meter used is 6.928 volts, giving a 
value for the proportionality constant, 
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Y/Ejjjax = (56. 6 tt/180) (1/6.928) 

= 0.1426 radians (4) 

volt 

Digital Voltmeter . The voltage recording system 
used is a Dymec Model DY-5552A voltage measuring and re- 
cording system. The system consists of a Model DY-2210 
voltage-to-frequency converter, a Hewlett-Packard model 
5211 A/E electronic frequency counter, a model DY-2540 
scanner/ coupler and a motorized tape punch. The voltage- 
to-frequency converter converts the voltage output of the 
ratio meter to a proportional frequency. This frequency is 
sampled and counted by the frequency counter. The display 
of this counter is the voltage accurate to the millivolt. 

The output of the counter is then scanned by the scanner/ 
coupler, and the paper punch records the voltage. The 
counter is triggered with a remotely located push button. 

Reference Step . The reference step. Figure 26, was 
fabricated for the comparison ref lectometer according to the 
design given by D. L. Hollway. The reference step is a 
gradual symmetrical E-plane taper. At the reference plane 
the guide suddenly resumes the standard guide dimensions re- 
sulting in a reflection at the reference plane. The mag- 
nitude and phase of the reflection coefficient of the ref- 
erence step are essentially constant over the band 8.0 to 
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Figure 26. Reference Step 

12.4 GHz. This was verified by the National Bureau of 
Standards over the band 8.2 GHz to 12.4 GHz, Figures 27 
and 28. The reference step was measured in a line termi- 
nated in a precision matched load, the two separated by a 15 
centimeter length of precision waveguide. Frequency marker 
pips are provided at 9, 10 and 11 GHz. 

A detailed error analysis of the sweep frequency cal- 
ibration results has not undergone full NBS review. How- 
ever, the uncertainty in the return loss magnitude is be- 
lieved to be within ± 2%. This is an uncertainty of less 
than 0.006 in reflection coefficient. The uncertainty in 
reflection coefficient angle is believed to be within ± 2.5 
degrees. These uncertainties were obtained by swept fre- 
quency measurements on a calibrated sliding load by NBS. 

The magnitude of the reflection coefficient for the 
reference step is calibrated at . 2 ± .011 (VSWR =1.5 ± 

•034) according to the results of Figure 27. In order to 



Figure 27 • Reference Step, Return Loss 


u> 
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simplify programming a constant value (0.2) has been used 
for the reference step. Figure 28 shows the phase angle va- 
riation to be negligible, therefore, the angle has also been 
considered a constant, and the phase of the disconuity meas- 
ured by the comparison ref lectometer is referenced to it. 

The variation present in the curve of Figure 28 is due to a 
phase variation of 2(31 where L is the separation between the 
HBS reference plane and the step measured. 

Typical Recorded Voltages . Voltages recorded by the 
comparison ref lectometer setup consist of two sets. The 
first set refered to as the reference set is taken with the 
reference step terminated with a matched load, Figure 29a. 

A second set of voltage readings is taken with the reference 
step terminated in the device to be measured, Figure 29h. a 
typical set of measurements for a symmetrical capacitive 
iris at 50 MHz intervals from 7.975 GHz to 12.425 GHz is 
given in Table 4 . 

The ref lectometer computer program. Appendix B, uses 
Equation (III-31) to calculate the total magnitude reflection 
coefficient of the reference step and the discontinuity to 
be measured as shown in Figure 30. The data points of 
Figure 30 are connected by a continuous curve, however 
| r | is only known at frequency points separated by 50 MHz, 
starting at 7.975 GHz. | r ( f ) j of Figure 30 is the function 
that is Fourier transformed from the wavenumber to the dis- 
tance domain using Equations (III-22, 23), where G(v) is 


^ 'BAAJ'b %£££££ A c £ /S& S Ttjr'* f<?SJ6ZL 


Figure 28. Reference Step, Phase Angle of the Reflection Coefficient 
Referenced to a Short 
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Reference I Precision Matched 

Step Guide Load 


(a) Reference Step Terminated by a Matched Load 


Precision Guide 



(b) Reference Step and Test Element Terminated by a 
Matched Load 


Figure 29. Measured Setup of Reference Step 



Table 4 


Typical Voltage Measurements, Read from Left to Right, and Top to Bottom 





— HE 

TZXZhCT 

VOlTASE 

DATA 









2.939 

2*051 

3.171 

2.775 

3.691 

3.176 

3.140 

2.830 

2.277 

2.719 

2.310 

2.916 

2.885 

2.713 

3.157 

2.561 
3 • 1 3 0 

2.575 

3.092 

T;999 

2.979 

“2.737* 

3.136 

3,377 

3.027 

3.236 

2.996 

3.161 

3.103 

3.463 

3.004 

3.017 

3.050 

3.010 

3.171 

3.161 

3.128 

2.856 
3 .252 

3.064 

3.352 

3.180 

3.205 

2 . 992 
3.253 

3.277 

3.159 

3.0 39 
2.999 

3.144 

3.041 

3.191 

3.200 

3.112 

3.135 

3.302 

3.180 

3.368 

3.358 

3.296 

3.306 

3.418 

3.247 

3.442 

3.350 

3.268 

3.314 

3.279 

3.264 

3.292 

3.338 

3.090 

3.315 

3.108 

3.192 

3*197 

3. 132 

3.025 

3.023 

'3.035 

2.976 

3.069 

3.194 

3.183 

3.219 

3.371 

3.335 

3.244 

3.289 






voltage data wITH UNKNOWN IN 
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1*10 5 
1.296 

1.870 

2.945 

2.328 

4.378 

4.400 

5.037 

5.359" 

4.776 

" 3 • 7 0 6 

3.690 

2.080 

1.990 

1.155 

1.028 

1.112 

2.263 

3.107 

4.291 

4.627 

5.329 

4.821 

5.471 

4.294 

4.876 

2.701 

3.161 

.980 

.991 

.'323 

2*350 

2.454 

.594 

4.156 

.909 

5.375 

3.244 

5.771 

5.235 

5.099 

6.093 

3.317 

6.159 

1.446 

5.601 

.729 

4.045 

1.553 

1.601 

3.758 

.382 

5.591 
1 .524 

6.167 

3.661 

5.685 

5.311 

4.339 
6.131 _ 

6.279 
3 • r>bl 

5.459 

5.540 

3.538 

5.399 

1.204 
6 • 4b4 

.281 

5.897 

1.510 

4.385 

4.002 

1.948 

5.880 

.162 

6.617 

.614 

6.516 

2.956 

5.646 

5.286 

3.899 

6.530 

1.519 

6.900 

.184 

6.565 

1.213 
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given by Equation (III-18) . 

Description of Instrumental Errors . All ref lectometers 
suffer from small imperfections such as a lack of perfect 
directivity in the directional couplers and discontinuities 
in the couplings and imperfections in the waveguides. Gen- 
erally, the two crystal detectors have different frequency 
responses, and these can cause changes in the output re- 
sembling those from reflections. The coupling coefficients 
of the directional couplers also vary slightly with fre- 
quency. 

In the comparison ref lectometer all measurements are 
made by comparison with a single known standard, the ref- 
erence reflection Y , and therefore, all of the instrumen- 

r 

tal defects are nearly cancelled out. 

Consider a test component having no internal re- 
flections. When the reference readings are taken, all the 
instrumental defects show up as a variation in E r with fre- 
quency. However, the test readings E will be equal to 
E r at every frequency and from (III-31) , G(v) is constant 
with frequency and a n and b n in (III-22, 23) are zero. 
Therefore, the instrumental errors have been cancelled out 
completely . 

When a test component includes reflections, a high 
degree of cancellation still exists. The instrument defects 
correspond to a pattern of reflections spaced at different 
distances from the reference. Because these contribute to 
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both readings, only a small residue can appear in G(v), and 
then only when a reflection in the test component T c is in 
the vicinity of an instrumental reflection In the worse 

case, when the two coincide in position and phase, it has 
been shown by Hollway that? 



(measured) 



(true) 



(true) 


(5) 


Since r.j is small this error usually may be neglected. 

A very serious error occurs should a phase lock be 
missed, and special care must be taken to ensure that this 
does not occur. If a phase lock is missed the entire set of 
measurements must be repeated, unless the value of the re- 
flection coefficient can be determined at the point that was 
missed. 

The method of taking measurements cancels out a great 
many errors; however, the drift with time of the ratio meter 
is not cancelled out. Therefore; it is important to take the 
two sets of measurements as quickly as possible to minimize 
this drift error. The effect of this drift is demonstrated 
in Figure 31a. Two sets of measurements were taken, each 
with the reference step terminated in a matched load sep- 
arated by precision guide. Equation (III-22, 23, and 26) 
were applied to the data of Figure 31a to obtain the modulus 
of reflection coefficient as a function of distance given in 
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Figure 31. Magnitude of Reflection Coefficient of (a) Plain Guide and Reference 
Step, (b) Noise Due to Ratiometer Drift 
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Figure 31b. The noise level of the system was always found 
to correspond to a reflection coefficient variation of less 
than .00005 as a function of distance. 

After calibration of the equipment, it was found that 
point reflections could be located in the waveguide within 
a few tenths of a millimeter and measured with an accuracy 
of ±3 percent in magnitude and ±5 degrees in phase angle; 
these results agree with those of Hollway.^ This phase 
error has been verified by comparing the phase of the meas- 
ured sample capacitive iris, inductive iris, metal induc- 
tive post, and dielectric post referenced to the unit con- 
ductance circle. These results are detailed in Chapter V. 


CHAPTER V 
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EXPERIMENTAL VERIFICATION OF COMPUTED RESULTS 

Introduction 

A set of experiments were conducted to validate the 
theoretically determined values of reflection coefficient and 
resultant mismatch. Computed results and measured results 
are presented in the following pages. As a verification of 
the validity of the waveguide element equations given in 
Chapter III a symmetrical capacitive iris, asymmetrical in- 
ductive iris, metal inductive post, and a dielectric post 
have been fabricated and were measured using the comparison 
ref lectometer and also using a slotted line. A standard 
waveguide mismatch was also measured by the comparison re- 
f lectometer and by the slotted line. In all cases the 
slotted line measurements, comparison ref lectometer measure- 
ments and calculated results are within the uncertainty of 
the comparison ref lectometer . 

The utility of the impedance matching program was 
tested in two demonstrations. In the first a capacitive 
iris was matched with a second capacitive iris as the match- 
ing element, and the resulting impedance match agreed with 
the expected result for the simple combination of waveguide 


elements . 
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In order to demonstrate the power of the impedance 
matching technique a shunt slot radiator was fabricated, 
measured, and impedance matched. A comparison of predicted 
and measured resultant reflection coefficients is pre- 
sented at selected frequencies. 

Reflection Coefficient of Selected Waveguide Elements 

The accuracy of the comparison ref lectometer was dem- 
onstrated with corroborative measurements of several ele- 
ments taken with both the comparison ref lectometer and a 
slotted line setup. The measurements also demonstrate the 
validity of Equations (III-34) through (III-48) as models 
of the waveguide matching elements in this research. 
Waveguide Standard 

An X-band standard disturbance, Ga. Tech Model 
SR120X, with a VSWR of 1.19 was measured by the comparison 
ref lectometer . Voltage measurements were recorded with the 
standard located 30.2 centimeters from the reference plane 
of the comparison ref lectometer . Calculations of the mag- 
nitude of the reflection coefficient as a function of dis- 
tance were made using Equation (III-31) and results have 
been plotted in Figure 32. The plot of Figure 32 appears 
as a continuous curve; however, only the discrete points at 
50 MHz intervals from 7.975 GHz to 12.425 GHz are precisely 
known. This curve shows the expected periodic dependence 
on frequency as a result of the constructive and destructive 
interference of waves reflected by the reference step and by 
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Figure 32. Magnitude of the Total Reflection Coefficient 
of a Waveguide Standard and Ref lectometer 
Reference Step 


the measured discontinuity. 

The magnitude of the average reflection coefficient 
generated as a function of distance was calculated by the 
computer program Test 1, listed in Appendix B. These cal- 
culations were made using the results shown in Figure 32 
and Equations (111^22), ( 111 - 23 ) and (III-26) , A curve of 
the reflection coefficient generated as a function of dis- 
tance is presented in Figure 33. An accurate electrical lo- 
cation was determined by measuring the distance at the 3 dB 
(half power points) and averaging the two distances. The 
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electrical location as found by this procedure was 30.16 
centimeters from the reference plane of the ref lectometer . 
The difference between .this apparent electrical location 
and the physical location is .04 centimeters. 



Figure 33. Electrical Location of a Waveguide Standard 


The complex reflection coefficient as a function of 
frequency is calculated at the electrical center by Equa- 
tions (III-22), (III-23) and (III-26) using overlapping in- 

tervals of integration as described in Chapter III and Appen- 
dix B. The computer program Test 2, Appendix B, implements 
the above equations to obtain the complex current reflection 



coefficient as indicated. Figure 34. Slotted line measure- 
ments were taken at the selected frequencies of 9.0 GHz, 

9.5 GHz, 10.0 GHz, 10.5 GHz, 11.0 GHz, and the maximum dif - 
ference between slotted line measurements and ref lectometer 
measurements is only 5%. This difference is within the un- 
certainty of the slotted line measurement. The phase of the 
waveguide standard is shown in Figure 34b on an expanded 
scale. The small variation of one to eleven degrees in 
phase angle as a function of frequency is characteristic of 
the waveguide E-plane step. 

Waveguide Matching Elements 

A symmetrical inguide capacitive iris, asymmetrical 
inguide inductive iris, solid inguide inductive post, and an 
inguide dielectric post were fabricated and their reflection 
coefficients measured. The results of these measurements 
were compared with the calculated reflection coefficients 
using Equations (III-34) through (III-48) for the respec- 
tive waveguide elements. 

Table 5 provides a comparison of the physically meas- 
ured location and the apparent electrical location of the 
various elements as measured by the comparison ref lectometer . 


MAGNITUDE OF RFH.FCU0N COEFFICIENT 




(a) Magnitude of Reflection Coefficient (b) Phase of Reflection Coefficient 

o Slotted Line Measurement 


Figure 34. Magnitude and Phase of the Reflection Coefficient of a Waveguide 
Standard E-Plane Step, Ga. Tech Model SR120X 




Table 5. Location of Fabricated. Waveguide Elements 


Elements 


Waveguide 
Location in 
Centimeters 

Ref lectometer 
Result Location 
in Centimeters 

Difference 

Capacitive 

Iris 

27.0 

26.94 

-0.06 

Inductive 

Iris 

27.0 

27.16 

+ 0.16 

Inductive 

Post 

27.0 

27.24 

+ 0.24 

Dielectric 

Post 

27.1 

27.02 

-0.08 

Waveguide 

Standard 


30.2 

30.16 

-0.04 


It is interesting to note that the electrical loca- 
tion of the capacitive iris appears .06 centimeters on the 
generator side of the geometrical location. The. waveguide 
standard, also capacitive, is located electrically .04 cen- 
timeters on the generator side of the geometrical location. 
This trend is also followed by the dielectric post which 
appears electrically to be .08 centimeters on the generator 
side of its geometrical center, and is also capacitive. 

On the other hand, the inductive post appears to be located 
.16 centimeters on the load side of the geometrical location. 
Therefore, it could be concluded that a capacitive element 
appears electrically to be in front of its physical posi- 
tion while an inductive element appears electrically to be 
behind its physical location. This interesting observation 
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has little consequence as far as impedance matching is 
concerned as long as the complex reflection coefficient is 
computed at the electrical center. 

The phase angle of the complex reflection coeffi- 
cient at the apparent electrical location is different from 
the phase at the physical location by a factor of 2gL 
radians where L is the difference in the two locations. Re- 
flectometer calculations of the complex reflection coeffi- 
cient were made for the apparent electrical location of each 
sample element; theoretical calculations (Equations (III-34) 
to (III-48)) were made for the respective physical location. 
The results obtained by the comparison ref lectometer should 
have the phase correction (2$L) applied before comparing 
with theoretically predicted results. In the following, 
comoarisons are made between the predicted calculated 
phase at the physical location with that of the measured 
phase at the electrical location. 

The fabricated symmetrical capacitive iris, Figure 14, 
(Chapter III) was measured by the comparison ref lectometer . 

A plot of the average magnitude of the reflection coeffi- 
cient generated as a function of distance appears as Figure 
35. The complex reflection coefficient was calculated at 
the electrical location of 26.94 centimeters.' This complex 
current reflection coefficient is displayed by the admit- 
tance Smith chart plot shown .in Figure 3 s . All Smith chart 
plots presented in this thesis are referenced to a matched 
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Figure 35. Distance Plot of a Symmetrical Capacitive Iris 

load termination. As the observer moves away from the loca- 
tion of the disturbance the phase of the complex reflection 
coefficient varies according to exp (-j2BL) where 3 is given 
by Equation (1-1) and L is the distance between the location 
of the disturbance and the observer. For a frequency of 
10.0 GHz the phase sensitivity is 182.88 degrees per centi- 
meter distance. The departure of the measured points in 
Figure 36 from the unit conductance circle can be attributed 
to the 0.06 centimeter difference between the physical center 
and the electrical center. This difference in distance re- 
sults in a phase angle difference of approximately 11 
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(a) Smith Chart Display of Current Reflection Coefficient 
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(b) Symmetrical Capacitive Iris 


Figure 36. Complex Current Reflection Coefficient of the 
Fabricated Capacitive Iris on Expanded 
Admittance Plot at the Electrical Location 
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degrees at 10.0 GHz. The actual phase difference measured 
from Figure 36 is approximately 5 degrees. This difference 
in phase is within the predicted accuracy of the comparison 
ref lectometer . ^ 

A comparison of the measurements taken of the fabri- 
cated capacitive iris and the calculated results are shown 
in Figure 37. In Figure 37, a curve of magnitude of reflec- 
tion coefficient is given as a function of frequency. Su- 
perimposed on this curve are the theoretical points and the 
slotted line measurements, both at selected frequencies. 
There is good agreement between the three sets of results. 
The phase, in Figure 37b, is compared with the calculated 
phase at selected frequencies. The measured results are 
in good agreement with the predicted characteristics and it 
can be concluded that the capacitive iris is adequately 
modeled by Equations (III-34) through (III-37) in the im- 
pedance matching computer program listed in Appendix B. 

The fabricated asymmetrical inductive iris of Figure 
15 was also measured by the comparison ref lectometer . A 
plot of the average magnitude of the reflection coefficient 
as a function of distance is given in Figure 38. 

The complex reflection coefficient was calculated 
at the apparent electrical location of 27.16 centimeters 
with respect to the ref lectometer reference plane. This 
complex reflection coefficient is given by the admittance 
Smith chart plot of Figure 39. The results of Figure 39 are 
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(a) -Ref lectometer Measurement 
+Calculated Result 
o Slotted Line Measurement 


Figure 37. Reflection Coefficient of a 
Thickness 0.081 cm and Iris 
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(b) -Ref lectometer Measurement ' 

+Calculated Result | 
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Fabricated Capacitive Iris with Iris 
Height = 0.145 cm 
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Figure 38. Distance Plot of an Asymmetrical Inductive Iris 

within the uncertainty of the comparison ref lectometer when 
electrical and physical center separation are taken into 
account. A comparison of the measurements taken of the fab- 
ricated inductive iris and the calculated results are shown 
in Figure 40. The measured magnitude of the reflection co- 
efficient is plotted as a function of frequency in Figure 
40a, and compared to theoretical points and slotted line 
measurements, both at selected frequencies. There is good 
agreement between the three sets of data. 

The phase in Figure 40b, is compared with the cal- 
culated phase at selected frequencies. The 29 degree dif-? 
ference in phase between the measured and calculated results 
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(a) Smith Chart Display of Current Reflection 
Coefficient 



(«? — 0 . 3 0 cm 


0.081 cm 







Cross Sectional View Top View 

(b) Asymmetrical Inductive Iris 

Figure 39. Complex Current Reflection Coefficient of the 
Fabricated Inductive Iris at the Electrical 
Location on Expanded Admittance Plot 
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(a) -Ref lectometer Measurement 
+Calculated Result 
o Slotted Line Measurement 


(b) -Ref lectometer Measurement 
+Calculated Result 


Figure 40. Reflection Coefficient of a Fabricated Inductive Iris with Iris 
Thickness = 0.081 cm and Iris Height = 0.30 cm 
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is largely due to the difference of 0.16 centimeters (equiv- 
alent to 29.8 degrees at 10 GHz) between the electrical and 
geometrical locations. Therefore, it can be concluded that 
the measured result is in good agreement with the predicted 
characteristics . 

A sample of the metal inductive post of Figure 16 , 
was given a similar exDerimental check. A plot of the aver^ 
age magnitude of the reflection coefficient as a function 
of distance is given in Figure 41. 



Figure 41. Distance Plot of a Metal Inductive Post 
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The complex reflection coefficient was calculated at the 
apparent electrical location of 27.24 centimeters. This com- 
plex reflection coefficient is displayed, Figure 42, by the 
admittance Smith chart plot. While the difference in phase 
between the calculated and measured results is 44 degrees at 
10. GHz, the difference in phase between the geometrical 
location and the electrical location is 43.89 degrees cor- 
responding to the 9.24 centimeters of Table 5. Therefore, 
there is a difference of only 0.11 degrees which is well 
within the accuracy of the ref lectometer . 

A comparison of the experimental and theoretical data 
for the fabricated inductive post is shown in Figure 43. In 
Figure 43a, a curve of magnitude of reflection is given as a 
function of frequency. Superimposed on this curve are the 
theoretical points and the slotted line measurements, both 
at selected frequencies. Again good agreement is found 
among the three sets of results. 

The dielectric post sample was fabricated from 
hot-pressed boron nitride, HD-0092. Boron nitride has a 
relative dielectric constant of 4.07 and a loss tangent of 
0.0003. 

The locating plot of the average magnitude of the 
reflection coefficient as a function distance is given in 
Figure 44. 
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(a) Smith Chart Display of Current Reflection Coefficient 
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(b) Metal Inductive Post 


Figure 42. Complex Current Reflection Coefficient of the 
Fabricated Metal Inductive Post at the 
Electrical Location on Expanded Admittance Plot 
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Figure 43. Reflection Coefficient of a Fabricated Inductive Post with Diameter 
0.137 cm and Sidewall Distance = 0.260 cm 
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Figure 44. Distance Plot of an Inguide Dielectric Post 


The complex reflection coefficient was calculated at 
the electrical location of 27.02 centimeters. For this ele- 
ment, the expected difference in phase resulting from the 
lack of coincidence of geometrical and electrical centers is 
14.6 degrees. The actual phase difference between the meas^ 
ured response and the unit conductance c.ircle is 16 degrees, 
shown in Figure 45. Therefore, the difference in predicted 
and measured phase is only 1.4 degrees. 

In a test of resolution of the technique a fabricated 
capacitive iris and an inductive iris were measured at the 
same time by the comparison ref lectometer . The capacitive 





Cross Sectional View 
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(b) Dielectric Post with Dielectric Constant s' = 4.07 


Figure 45. Complex Current Reflection Coefficient of the 
Fabricated Dielectric Post at the Electrical 
Location on Expanded Admittance Chart 
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(a) -Ref lectometer Measurement 
+Calculated Result 
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(b) -Ref lectometer Measurement 
+Calculated Result 


Figure 46. Reflection Coefficient of a Fabricated Dielectric Post with Diameter 
0.23 cm. Sidewall Distance = 0.80 cm, and Dielectric Constant = 4.07 
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iris was located 27.00 centimeters from the reference plane - 
and the inductive iris 51.00 centimeters from the reference 
plane. A curve of the total magnitude of the reflection co- 
efficient of the reference step, capacitive iris, and induc- 
tive iris is given by Figure 47a. From this information the 
average magnitude of the reflection coefficient generated 
as a function of distance is calculated. Figure 47b. The 
electrical location for the capacitive iris is calculated to 
be 26.98 centimeters, and that of the indhctive iris 51.16 
centimeters. The difference in the geometrical and electri- 
cal location of the capacitive iris is .01 centimeters. In 
Table 5, for a single capacitive iris in the guide the dif- 
ference is 0.^6 centimeters. 

The distance between the electrical and geometrical 
locations for the inductive iris is 0.16 centimeters and 
compares with the 0.16 cm difference of Table 5. This good 
agreement between the single disturbances and multiple dis- 
turbance tests made by the ref lectometer tends to justify 
the approximations made in Chapter III, Equation (14). 

Impedance Matching 

Impedance Matching a Capacitive Iris 

In a trial test of the impedance matching program a 
single capacitive iris was used as the waveguide disturbance. 
A second capacitive iris was selected as the matching ele- 
ment. The disturbance was impedance matched over three dif- 
ferent bandwidth s to demonstrate the adaptability of the im- 
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pedance matching program. The iris selected as the mis- 
match has a thickness of 0.081 centimeters and an iris 
height of 0.145 centimeters. 

If an impedance match were calculated using a Smith 
chart as a tool, generally, a capacitive iris of the same 
dimensions would be placed a quarter wavelength from the 
first iris, where the wavelength corresponds to the center 
frequency of the matched bandwidth. Using this procedure 
there would be no control of the resulting mismatch over the 
rest of the band. The impedance matching computer program, 
Appendix B, calculates the combination of physical dimen- 
sions and relative distance between elements which results 
in the least root mean square mismatch over the bandwidth 
to be matched. Therefore, the technique developed in this 
thesis provides some control over the resulting mismatch 
over the entire bandwidth. This can be seen by considering 
Figure 48. In Figure 48 the VSWR of the initial mismatched 
capacitive iris is given together with three other curves 
of VSWR, one for each match band. 

The narrowest band provides the smallest mismatch 
at the center of the band, but the mismatch increases more 
rapidly than that of the broader band impedance match. The 
narrowest band match is from 10.195 to 10.595 GHz, with a 
center frequency of 10.395 GHz. At the center frequency the 
VSWR is equal to 1.0006. The next broader bandwidth is from 
9.795 GHz to 10.995 GHz, also with a center frequency of 
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FREQUENCY IN GHZ . 


(1) Bandwidth 10.195 - 10.595 GHz 

(2) Bandwidth 9.795 - 10.995 GHz 

(3) Bandwidth 8.595 - 11.595 GHz 


Fxgure 48. 5 h f r u C ^ eristics of a Ca Pacitive Iris Impedance 
Matched Over Three Bandwidths 
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10.395 GHz. At the center frequency the VSWR is equal to 

1 . 022 . 

While the effectiveness of the broad band impedance 
matching program is demonstrated in these test runs , the 
restriction of the matching element to a second capacitive 
iris makes the change in resultant mismatch small. 

Shunt Slot Radiator 

The shunt slot radiator has been selected to demon- 
strate the power of this computer aided impedance matching 
technique. Interest in the slot was primarily motivated by 
by the need to broadband impedance match this element when 
it is used in a broad side steerable array. While this re- 
search is concerned primarily with inguide disturbances, 
the radiating slot provides an interesting application of 
this matching technique. 

The ref lectometer measured the slot apparent electri- 
cal location as 28.68 centimeters from the reference step 
as shown in Figure 29. The geometrical center of the slot 
is located 27.0 centimeters from the ref lectometer refer- 
ence step. . Therefore the electrical center is located 1.6 8 
centimeters on the load side of the physical center of the 
slot. 

The shunt slot is a radiating device, and there are 
surface currents on the outer surface of the waveguide in 
the vicinity of the slot. These surface currents together 
with the disturbed surface currents on the inside wall of 


no 


tne slot, determine the electrical location of -the slot as 
measured by the comparison ref lectometer . 



Figure 49. Distance Plot of the Shunt Slot Radiator 


Unlike the nonradiating elements, the shunt slot 
radiator showed a marked separation between the electrical 
location and the geometrical location. 

In Figure 50a, and 50b the magnitude and phase of 
the complex reflection coefficient at the apparent elec- 
trical location are given. These results were calculated 
by the comparison ref lectometer programs, Test 1 and 2. 

It is possible to calculate the complex reflection coeffi- 
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Figure 50. Magnitude and Phase of the Reflection Coefficient of a Shunt Slot Radiator 
at the Apparent Electrical Location 
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cient two ways at the physical location. The complex re- 
flection coefficient can be calculated at the apparent elec- 
trical location and multiplied by exp (-j28(1.68 cm)) re- 
sulting in the correct magnitude and phase at the physical 
center of the slot. On the other hand, the complex reflec- 
tion coefficient can be calculated directly at the physical 
location by the program Test 2 . This calculation at the 
geometrical location results in an erroneous result. The 
error stems from the fact that the locating nature of the 
transforms of Equation (III-22) and (III-23) attenuate the 
magnitude of the reflection coefficient at distance points 
not equal to the electrical location, shown in Figure 49. 

To test for the existance of secondary regional 
matching minima the slot radiator was matched by a capaci- 
tive iris using several different initial positions for the 
matching iris. Ideally, the matching program should find 
the same match point (location, iris dimension and RMS 
mismatch) independent of the starting point, in the ab- 
sence of secondary minima. These results are tabulated in 
Table 6. The relative starting locations are the electri- 
cal center, (0,0), and ±0.5566, ±1.1133, ±1.6699 and 
±2.2265 centimenters from the electrical center. These 
distances correspond to 1/8 wavelength, 1/4 wavelength, 3/8 
wavelength, and 1/2 wavelength magnitudes respectively at the 
frequency of 9.395 GHz. The frequency of 9.395 GHz was 
chosen for convenience. 
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Table 6. Capacitive Iris Matching of the Slot Radiator 
Referenced to the Slot Electrical Center 


Initial 

Distance 

From 

Electrical 
Center 
of Slot* 

RMS 

Mismatch 

Iris 

Thickness* 

Iris 

Height* 

Iris 

Distance 

From 

Electrical 
Center 
of Slot* 

2.2265 

0.105 

0.054 

0.154 

,2.177 

1.6699 

0.105 

0.041 

0.160 

2.182 

1.1133 

0.064 

0.044 

0.185 

0.013 

0.5566 

0.063 

0.095 

0.160 

0.032 

0.0000 

0.063 

0.095 

0.160 

0.025 

-0.5566 

0.063 

0.095 

0.160 

0*031 

-1.1133 

0.064 

0.044 

0.185 

0.012 

-1.6699 

0.098 

0.050 

0.160 

-2.182 

-2.2265 

0.098 

0.050 

0.160 

-2.176 


*In centimeters 


At any selected frequency point the complex reflec- 
tion coefficient varies in phase as a function of distance 
from the generating disturbance. This variation can be 
expressed in the form of Equation (1) : 


r(f,L) = r c (f)exp[-j2$L] , 


( 1 ) 


- JJ.4 


where 3 is -the propagation constant and F c (f) is the com- . 
plex reflection coefficient of the disturbance when L = 0. 
The greater the distance L between the matching element and 
disturbance the more rapid the change in phase with fre- 
quency . 

As a result of this variation in phase of the com- 
plex reflection coefficieint as a function of distance, it 
is reasonable to expect a least root mean square mismatch 
(LRMS) to occur for matching iris placements at intervals 
of Ag/2 distance as shown in Figure 51, where Xg is the 
guide wavelength of the center frequency of the matching 
bandwidth. However, because the RMS is taken over the en- 
tire bandwidth of interest, this distance is not expected 
to be exact. As the magnitude of the distance L increases 
so does the phase change with frequency. Therefore, the 
quality of the impedance match is expected to change as the 
matching element is located at distances that are integral 
multiples of X^/2. It is important to notice the absence 
of large secondary minima when initiating calculations 
from different relative separation distances L between the 
disturbance and the matching element. All five points of 
Table 6 that are 1/4 wavelength or less from the electrical 
center resulted in virtually the same match point. At each 
end the two most distant starting points again resulted in 
essentially identical match points each located approximate- 
ly Xg/2 from the central point. 
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• Regional Minima 
O Starting Points 



A_/2 {< Separation 

Distance 


Figure 51. Schematic Plot of RMS Mismatch as a Function of 
Separation Distance L, Showing the Absence of 
Secondary Minima Within Ag/2 Intervals Centered 
at Regional Minima, Initial Starting Position of 
Table 6 Are Sketched 


Better RMS impedance matches are accomplished by plac- 
ing the capacitive iris on the load side of the geometrical 
center of the slot. In fact, the lowest PMS impedance match 
is achieved by placing the capacitive iris approximately 
0.03 centimeters on the generator side of the apparent elec- 
trical center of the slot. 

In all of the broad band impedance matching done in 
this research, the best broad band matches were achieved with 
the matching element placed in the immediate vicinity of the 
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electrical center. This result is to be expected, since a 
large separation distance L between disturbance and matching 
element results in an interference pattern which is generally 
dominated by the exp(-j26L) phase factor of Equation (HI- 
49), and is narrow band in nature. 

In Table 6 the starting positions of 0.5566, 0.0000, 
and -0.5566 centimeters results in an impedance match which 
places a capacitive iris of equivalent dimensions at prac- 
tically the same location. The slight variation in the fin- 
al location is due to the final increment size of the rela- 
tive distance, ALENGH in the impedance matching program de- 
scribed in Appendix B. Starting positions of 1.6699 and 
2.2265 centimeters in Table 6 result in the placing of an 
iris at the same location of approximately 2.18 centimeters 
resulting in an RMS mismatch of 0.105. This is not physi- • 

cally realizable, however, since the iris would be located in 
the slot. All other entries in Table 6 are physically real- 
izable because they are either located at the apparent elec- 
trical location or 2.18 centimeters on the load side of the 
apparent electrical location. From Table 6 it is easily 
recognized that when impedance matching the slot by a capa- 
citive iris , superior results are achieved by impedance 
matching in the vicinity of the electrical center. Table 
6 agrees with Figure 51 and demonstrates the absence of 
large secondary minima within relative distances which are 
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multiples of A g /2. 

Impedance Matching the Slot Radiator 

For completeness the characteristics of the shunt 
slot radiator were impedance matched by each of the four 
matching elements at the apparent electrical location* 

Table 7 provides a tabulation of the calculated results of 
the slot relative to the apparent electrical location. The 
relative quality of the impedance match is indicated by the 
RMS mismatch achieved using each matching element. The 
physical dimensions of each matching element are given, and 
the placements toward the generator are given. The two 
best matches occur with the inductive metal post and the 
dielectric post with a dielectric constant of 4.07. The 
metal post located 1.175 centimeters toward the load matches 
the slot to an RMS mismatch of 0.056, Figure 52a. The di- 
electric post located 0.012 centimeters toward the load 
matches the slot to an RMS mismatch of 0.058, Figure 52b. 
From Figure 52a the resultant VSWR at 8.551 GHz is 1.177, 
and at the upper limit of the matching band, 10.560 GHz, 
it is 1.184. In all cases the slot radiator has been im- 
pedance matched over 2 GHz bandwidth, 8.551 GHz to 10.560 
GHz. From Figure 52b the VSWR at 8.551 GHz is 1.09 and at 
the upper limit of the band, 10.560 GHz, it is 1.188. Both 
the metal inductive post and the dielectric post provide 
good impedance matches, less than 1.2 VSWR over the chosen 
2 GHz bandwidth. It is interesting to note that the metal 
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inductive post matches the slot radiator to less than 1.23 
VSWR over the entire band from 8.551 GHz to 11.934 GHz. 


Table 7. Slot Radiator Impedance Matched by a Single 
Matching Element Referenced to the Apparent 
Electrical Center 


Distance 


Matching 

Element 

RMS 

Mismatch 

Physical 

Dimension 

Physical 

Dimension 

Toward 

Generator 

in 

Centimeters 

Capacitive 

Iris 

0.063 

0.095* 

0.160** 

0.025 

Inductive 

Iris 

0.072 

0.051* 

0.335** 

1.187 

Inductive 

Post 

0.056 

0.010+ 

0.315++ 

-1.175 

Dielectric 

0.058 

0.297+ 

0 . 477++ 

-0.025 


* Thickness in Centimeters 
** Height in Centimeters 
+ Diameters in Centimeters 
++ Sidewall Distance in Centimeters 


The above results can be improved by use of the impe- 
dance matching program to successively match the slot with 
two elements. The slot is first matched with a single match- 
ing element, and the results of this impedance match are 
then recorded on data cards in the form of the complex re- 
flection coefficient of the resulting mismatch. If the 
matching element is located on the load side of the distur- 
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(b) Slot Matched by Dielectric Post 


Figure 52. Slot Radiator Impedance Matched by (a) Metal Inductive Post 
(b) Dielectric Post 
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bance then the reference plane of the reflection coefficient 
of the resultant match is that of the original disturbance. 
When the matching element is located on the generator side 
of the disturbance, the reference plane of the resultant mis- 
match is that of the matching element. 

The results of the first match are then used as data 
cards, i.e., the disturbance to be matched, and an addition- 
al matching element is calculated by the impedance matching 
program. Following this strategy, the results of the four 
impedance matches of Table 7 were matched with an addition- 
al element. The results of this impedance o match are tabu- 
lated in Table 8 as the RMS mismatch over the 2. GHz band- 
width of interest. The element used to match the slot the 
first time is listed in the first column. The second 'match- 
ing element used is listed by the first row. For example the 
first number in the second column of numbers, 0.063 is the 
RMS mismatch resulting from matching the slot with both a 
capacitive iris and an inductive iris located in the guide. 
The first matching element is the capacitive iris and 
the second is the inductive iris. There are four com- 
binations which produce good matches as seen in Table 8. 

These combinations are the inductive post-capacitive iris 
combination with an RMS value of 0.047, inductive post- 
inductive iris combination with an RMS of 0.056, inductive 
post-dielectric post combination with an RMS of 0.049 and 
dielectric post-capacitive iris combination with an RMS of 
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0. 056. 


Table 8. Slot Radiator Impedance Matched by Two Matching 
Elements 


Second Capacitive Inductive Inductive Dielectric 

\ Match Iris Iris Post Post 


FirstV 
Match ^ 

• 




Capacitive 

Iris 

xxxxx 

0.063 

0.093 

0.061 

Inductive 

Iris 

0.071 

xxxxx 

0.078 

0.072 

Inductive 

Post 

0.047 

0.056 

xxxxx 

0.049 

Dielectric 

Post 

0.056 

0.058 

0.088 

xxxxx 


The curves of VSWR as a function of frequency are 
given for the above combinations in Figures 53, 54, 55 and 
56, respectively. The relative location of the matching 
elements with respect to the slot is given by the (a) part 
of the respective figures, and the VSWR is given by the (b) 
part of the respective figures. 

The dielectric post-capacitive iris matching combina- 
tion was fabricated to compare the computer predicted result- 
ant mismatch with an experimentally measured resultant mis- 
match. The boron nitride was machined into a post with a 
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(a) Location of Matching Elements 
Post Diameter = 0.01 cm 
Post Sidewall Distance = 0.3150 cm 
Iris Thickness = 0.120 cm 
Iris Height = 0.060 cm 


(b) Response of a Slot Matched by a 
Metal Post and by a Metal Post 
Capacitive Iris Combination 


Shunt Slot Radiator Impedance Matched by a Metal Post and by a Metal Post 
Capacitive Iris Combination 
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(a) Location of Matching Elements 
Post Diameter - 0.01 cm 
Post Sidewall Distance = 0.3150 cm 
Iris Thickness = 0.0375 cm 
Iris Height = 0.060 cm 


Figure 54. Shunt' Slot Radiator Impedance 
Post-Capacitive Iris Corabinat 
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(b) Response of a Slot Matched by a 
Metal Post and by a Metal Post- 
Capacitive Iris Combination 
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(a) Location of Matching Elements 
Metal Post Diameter . 01 cm 
Dielectric Post Diameter .1538 cm 
Metal Post Sidewall Distance .3150 cm 
Dielectric Post Sidewall Distance .4637 


(b) Response of Slot Matched by a Metal 
Post and by a Metal Post-Dielectric 
Post Combination 


Figure 55. Shunt Slot Radiator Impedance Matched by a Metal Post and by a Metal Post- 
Dielectric Post Combination 
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(a) Location of Matching Element 

Dielectric Post Diameter .2975 cm 
Dielectric Post Sidewall Distance 
.4775 cm 

Capacitive Iris Thickness .035 cm 
Capacitive Iris Height .060 cm 


Figure 56. Shunt Slot Radiator Impedance Matched by a Dielectric Post and by a 
Dielectric Post-Capacitive Iris Combination 


(b) Response of a Slot Matched by a 
Dielectric Post and by a 
Dielectric Post-Capacitive Iris 
Combination 
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diameter of 0.117 ± .002 inches and inserted between the 
walls of the waveguide to the position given in Figure 56. 
While the tolerance on the post location was difficult to 
measure, it is estimated to be within 0.002 inches of the 
desired location. The symmetrical capacitive iris was ma- 
chined to fit two slots sawed through the broad side of the 
waveguide. The slots are 0.014 ±.002 inches wide and 0.074 ± 

.002 inches deep. The dielectric post was positioned in the 
guide, and the capacitive irises were fitted into the respec- 
tive slots. A silver based conductive paint was painted on 
the outside wall of the waveguide over the ends of the irises 
and on the intersection of the iris and the inside waveguide 
walls . 

The impedance matched shunt slot was then measured 
by the slotted line at 100 MHz intervals starting at 8.575 
GHz through 10.875 GHz. The results of these measurements 
are tabulated in Table 9. The results are also given in 
Figure 57 where the calculated impedance match is plotted 
as a function of frequency and the measured results are 
plotted at discrete points. The excellent agreement dis- 
played between the computer calculated VSWR and the experi- 
mentally measured VSWR, validate the impedance matching com- 
puter program and in turn the entire computer aided broad 
band impedance matching technique. 

The Effect of Fabrication Tolerances 

While mechanical tolerance is of little importance in 
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Table 9. VSWR Measurements of the Fabricated Impedance 
Matched Shunt Slot Radiator 


Frequency in GHz 

Measured VSWR 

8.575 

1.07 

8.675 

1.08 

8.775 

1.07 

8.875 

1.12 

8.975 

1.14 

9.075 

1.155 

9.175 

1.16 

9.275 

1.14 

9.375 

1.145 

9.475 

1.115 

9.575 

1.07 

9.675 

1.04 

9.775 

1.02 

9.875 

1.015 

9.975 

1.02 

10.075 

1.04 

10.175 

1.07 

10.275 

1.08 

10.375 

1.12 

10.475 

1.13 

10.575 

1.13 

10.675 

1.15 

10.775 

1.17 

10.875 

1.19 


V0LTRGF. STANDING WAVE RAT 1 0 


128 



u.1.992 

0.98 cm k c ^} 



Figure 57. Comparison of Calculated Mismatch 
and Measured Mismatch 
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theoretical studies, the allowable tolerance greatly influ- 
ences the cost of a practical device. Certainly the spec- 
ified tolerances should depend on the effect of the dimen- 
sional variations on the resultant mismatch. This is true 
because generally closer tolerances increase the expense of 
fabrication. 

In Table 10 the RMS mismatch values for the final 
nearest neighbor points for the slot impedance matched by a 
dielectric post are given. The final increment size for 
the sidewall distance of the dielectric post is 0.00125 
centimeters, for the post diameter it is 0.00625 centimeters, 
and for the relative separation distance it is 0.0125 centi- 
meters. Table 10 shows the variation in the RMS mismatch 
as a result of a change in a single dimension holding the 
remaining two dimensions equivalent to the LRMS or (1,1,1) 
point dimensions. Since the unequal increment sizes were 
selected so as to produce equivalent changes in RMS mis- 
match, it is difficult to determine the relative sensitivi- 
ty of this mismatch to changes in a single parameter. 
Therefore, it is of further interest to examine the sensi- 
tivity of the LRMS mismatch as a function of equal incre- 
ment variations of each physical dimension, since this sen- 
sitivity determines the importance of machining tolerances. 

An examination of the effect of machining tolerances 
on the resultant mismatch was conducted in two steps. First 
the slot was theoretically impedance matched by a dielectric 


Table 10. Final RMS Mismatch Values for the Nearest 
Neighbor Points of the Dimensional Space 
Array for the Case of the Slot Matched by 
a Dielectric Post 


Coordinates* » 
(I,J,K) 

Sidewall 

Distance 

Diameter 

Separation 

Distance 

RMS 

Mismatch 

(1,1,1) 

0.47750 

0.29750 

-0.0125 

0.058713 

(1,1,2) 

0.47750 

0.29750 

-0.025 

0.058813 

(1,1,3) 

0.47750 

0.29750 

+0.000 

0.058982 

(1,2,1) 

0.47750 

0.29125 

-0.0125 

0.058790 

(1,3,1) 

0.47750 

0.30375 

-0.0125 

0.059385 

(2,1,1) 

0.47625 

0.29750 

-0.0125 

0.058715 

(3,1,1) 

0.47875 

0.29750 

-0.0125 

0.058717 


* I Corresponds to sidewall distance of the dielectric 
post. 

J Corresponds to diameter of the dielectric post. 

K Corresponds to relative separation distance. 
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post with various dimensions. The dimensions consisting of 
the optimum dimensions determined by the impedance matching 
program and variations of these dimensions by four thou- 
sands of an inch. Because the diameter of the post is the 
easiest dimension to machine to ± 0.001 inches, the diameter 
was held constant and the sidewall distance and relative 
distance between disturbances was varied. The dielectric 
post material used is boron nitride grade HD^0092 hot 
pressed. The average reflection coefficient over the 2. GHz 
bandwidth from 8.551 to 10.560 GHz was used to compare the 
quality of the resulting matches obtained. Table 11 gives a 
tabulation of the dimensional combinations used and the 
average theoretical reflection coefficient obtained. From 
Table 11 the largest average matched reflection coefficient 
0.0544 has a corresponding diameter of 0.2975 centimeters, 
sidewall distance of 0.4673 centimeters, and a relative 
displacement of 0.0023 centimeters. This combination will 
be retained and used in the next step of the perturbation 
examination, since the strategy is to obtain an indication 
of the worst effect of tolerances. 

The second phase of this examination consists 
of using the combination of dimensions which produce the 
largest average matched reflection coefficient and perform- 
ing a second match on this result using a capacitive iris as 
the matching element. The effect of tolerances is studied 
by varying the capacitive iris height and relative distance 
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from the first match determined in the first step of this 
study. The thickness of the capacitive iris is held con- 
stant in this study since this is generally determined by a 
commercially available shim stock and the match is relative- 
ly insensitive to ± 0.002 inch variation. The physical di- 
mensions used for the capacitive iris are iris thickness 
0.0350 centimeters, iris height 0.060 ± .01016 centimeters 
and iris distance from electrical location of slot 0.70 ± 
.01016 centimeters. Table 12 illustrates the average total 
matched reflection coefficient and the respective physical 
dimensions . 

Table 11. Average RMS Mismatch Over the Band of 8.55 GHz 
to 10.56 GHz Resulting from the First Phase of 
the Machining Tolerance Study, Dielectric Post 
Diameter =. 0.2975 Centimeters, e' = 4.07 


.. Relative 

Displacement* 

-0.0023 

-0.0125 

-0.0227 

Sidewall^ 

Distance* 




0.4673 

0.0544 

0.0537 

O'. 0533 

0.4775 

0.0542 

0.0536 

0.0532 

0.4877 

0.0542 

0.0535 

0.0532 


*In Centimeters 
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Table 12. Average RMS Mismatch Over the Band of 8.55 
GHz to 10.56 GHz Resulting from the Second 
Phase of the Machining Tolerance Study, 
Capacitive Iris Thickness = 0.0350 Centimeters 


Iris 

Height* 

Relative 

Displacement* 

\ 

0.6898 

0.7000 

0.7101 

0.0498 


0.058 

0.057 

0.054 

0.0600 


0.055 

0.054 

0.054 

0.0701 


0.055 

0.055 

0.055 

*In Centimeters 





The largest mismatch 

listed in Table 

12 results 

from 


the capacitive iris dimensions of iris thickness 0.0350 
centimeters, iris height 0.0498 centimeters and relative 
displacement 0.6898 centimeters. Figure 58 gives curves of 
the first impedance match by the dielectric post resulting 
in an average mismatch reflection coefficient of 0.054 and 
the results of matching the first match by a capacitive iris 
giving an average reflection coefficient of 0.058. 

It is interesting to note that while the dielectric 
post-capacitive iris combination match has a higher aver- 
age reflection coefficient than the dielectric post match, 
the combination match is a better broad band match over a 
larger bandwidth of interest. 

From Table 11 the most sensitive dimension to change 
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Figure 58. Comparison of the VSWR of a Slot Radiator 
Impedance Matched by a Dielectric Post and 
by a Dielectric Post-Capacitive Iris 
Combination 


is the relative distance along the waveguide. In the second 
row first column of numbers in Table 11 the relative dis- 
placement is -0.0023 centimeters, that is, 0.0023 centimeters 
toward the load from the electrical location of the slot 
radiator. The average reflection coefficient is 0.054. 

From Table 12 the most sensitive parameter is the iris 
height. In the second row of numbers, reducing the iris 
height to 0.0498 centimeters increases the average reflec- 
tion coefficient to 0.0579. Therefore, close tolerances 
are indicated for the relative distance between the dielec- 
tric post and the electrical location of the slot radiator. 
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Close tolerances are also indicated for the capaci- 
tive iris height. However, the overall results indicate 
that the entire matching process is relatively insensitive 
to machining tolerances of ± 4 thousanth of an inch or 
± 0.01016 centimeters. 
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CHAPTER VI 

CONCLUSIONS AND RECOMMENDATIONS 
Introduction 

Several important conclusions can be drawn as a 
result of this study. They deal with the experimental 
equipment and with the computer aided matching program. The 
following paragraphs highlight these conclusions 

Conclusions 

As a result of this research, an economical, fast, 
and reliable impedance matching technique has been estab- 
lished which can provide broad band impedance matches. 

The measured results of the comparison ref lectometer 
in resolving waveguide disturbances and in the determination 
of their complex reflection coefficients have been excellent. 
Ref lectometer measurements have been compared with both 
theoretical calculations and slotted line measurements. 

Measurements of the sample waveguide discontinuities 
have pointed out the difference in the apparent electrical 
location and the geometrical location. While the two lo- 
cations were closely spaced for the inguide disturbances , 
measured as described in Chapter V, there was a consider- 
able difference in the locations for the shunt slot radiator. 
The shunt slot radiator was impedance matched relative to 
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the apparent electrical location both theoretically and 
experimentally over a 2 GHz bandwidth as shown in Figure 56. 
The excellent agreement between the theoretical and experi- 
mental impedance matches demonstrates the validity of this 
research. 

The impedance matching computer program demonstrated 
an ability to calculate the physical dimensions of a pre- 
determined matching element in order to achieve an impedance 
match which is best over a specified bandwidth in the RMS 
mismatch sense. Chapter III. 

It has been shown that the computer matched solu- 
tions are relatively free of secondary RMS minima. Region- 
al minima are found at intervals of approximately half 
guide wavelength separations between the measured distur- 
bance and the matching element, described in Chapter V, 

Table 6. 

For the cases investigated the impedance match im- 
proved in the vicinity of the apparent electrical center of 
the test disturbance. The shunt slot radiator was impe- 
dance matched to a VSWR of 1.16 to 1.0 over a 2 GHz band- 
width. This match was fabricated to demonstrate the valid- 
ity of the impedance matching program. The agreement be- 
tween the predicted and measured results for the matched 
slot was very good. 

It was found that while machining tolerances affect 
the calculated impedance match, for a reasonably good 
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matching tolerance of ± 0.002 of an inch the results are 
relatively insensitive. This conclusion increases the 
practical significance of this research, by providing a 
technique that is practically realizable. 

Recommendations 

There are several areas in which this research can be 
extended. Extension of the theory of the comparison 
ref lectometer , to enable it to measure larger disturbances, 
both in magnitude of reflection coefficient and in physi- 
cal length, can be made. The impedance matching computer 
program can be extended by both adding additional matching 
elements and enabling the program to calculate the location 
and physical parameters of more than one matching element 
simultaneously. A modification of the matching criterion 
could be implemented. The program could also be extended 
to select the most favorable type of matching element. 
Extending the Comparison Ref lectometer 

The comparison ref lectometer fabricated for this re- 
search is capable of locating disturbances resulting in re- 
flection coefficients less than or equal to 0.2. This limita- 
tion is in part due to the approximations made in the develop- 
ment of the theory and in part due to the reference step fab- 
ricated, since the development assumes that the tested reflec- 
tion coefficient does not exceed the reference reflection coef 
ficient. Secondly an equipment limitation is present in the 
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ratio meter used. The combined magnitude of the reflection 
coefficient of the reference step and the test disturbance 
of 0.2 reflection coefficient results in a full scale de- 
flection of the ratio meter. Measuring reflection coeffi- 
cients greater than 0.2 would result in off scale deflec- 
tions. For single frequency measurements the range of 
the ratio meter may be adjusted to prevent overdriving the 
unit. This is not possible in the comparison reflectom- 
eter system, however, since the two sets of measurements 
taken are compared to reduce instrumental errors. 

The length of line which can be measured is dependent 
on the sample intervals. At present the comparison reflec- 
tometer is phased locked, and a sample is taken every 50 MHZ. 
As such, the approximate limit on the line length measured 
is one meter. This can be increased to 5 meters by phase 
locking and taking a reading every 10 MHz. The closer the 
samples are taken in the wavenumber domain the greater the 
relative allowable distance between the reference step and 
test element. 

The Impedance Matching Program 

The impedance matching program was developed to match 
an arbitrary waveguide element as measured by the comparison 
ref lectometer . The impedance match is performed by calculat- 
ing the physical dimensions of a single preselected matching 
element. The criterion used by the matching program is the 
LRMS reflection coefficient resulting from the combined 
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wave reflected by the disturbance and the matching element. 

It may be desirable in some applications to use a 
maximum point VSWR criterion. That is, the matching pro- 
gram determines the maximum VSWR over the matching band- 
width for each of the dimensional space nearest neighbor 
points. The point selected as the (1,1,11 point would be 
that nearest neighbor point which has the smallest maximum 
VSWR. 

The matching program could also be extended to auto- 
matically select an element among the programmed matching 
elements which would provide the best broad band impedance 
match. This could be accomplished by assuming a broad band 
match and solving Equation (III-49) for the complex re- 
flection coefficient of the required matching element. 

The prospective matching element that compares most favora- 
bly with the required complex reflection coefficient could 
then be selected. This selection could be made by loosely 
comparing the slope and phase of each available matching 
element with the necessary characteristics calculated from 
Equation (III-49). 

The function of the matching program implemented in 
this work is to minimize the VSWR resulting from the test 
element and a single matching element using the RMS cri- 
terion. To this end the program functions properly. In 
some cases it is important to achieve a lower VSWR mis- 
match than is obtainable with the programmed matching ele- 
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ments. For these cases it is desirable to use two or more 
matching elements to achieve the desired impedance match. 
This, in fact, has been demonstrated in Chapter V when 
matching the shunt slot radiator. 

It is important to recognize that the first impedance 
match using a dielectric post was calculated as the LRMS 
over the selected 2 GHz bandwidth, subject to increment size. 
The result of the procedure was then improved upon by the 
capacitive iris by again impedance matching this initial 
result using the LRMS criterion. The slot and two matching 
elements are not a true optimum, however, since the two 
matching elements were not determined simultaneously. 

To demonstrate that a lower overall reflection co- 
efficient is possible, the mismatch resulting from the di- 
mensional combinations of Table 11 were used as the input 
data to the impedance matching program. A capacitive iris 
was sized for each case providing an LRMS mismatch. The 
results of this study are given in Table 13. The first 
row listed in Table 13 corresponds to the original com*' 
bination of Figure 55a. This is also the combination of 
elements which has been fabricated, and a comparison of 
theoretical and experimental results are shown in Figure 56. 
There are six combinations in Table 13 , resulting from per- 
turbations of the dielectric post, which provide a lower RMS 
mismatch than that originally obtained. A comparison of 
the response of the initial combination and the combina- 
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tion of matching elements resulting in an RMS mismatch of 
0.054 is given in Figure 59. 

While the improvement observed is small it must be 
remembered the dimensional changes were constrained to be 
small. The results do indicate clearly that if more than 
one element is needed to provide a desired broad band imped- 
ance match an extension of the program to simultaneously 
adjust more than one element would be desirable. 

Table 13. RMS Mismatch Using a Capacitive Iris 
as the Second Matching Element 


Relative Distance 


RMS 

Mismatch 

Iris 

Thickness 

Iris 

Height 

from Electrical 
Center 

0.056 

0.0350 

0.0600 

0.7000 

0.055 

0.0575 

0.0600 

0.6875 

0 . 056 

0.0300 

0.0600 

-0.6275 

0.055 

0.0512 

0.0600 

0.7625 

0.057 

0.0412 

0.0538 

0.6375 

0.054 

0.0737 

0.0600 

0.7375 

0.056 

0.0400 

0.0600 

-0.6875 

0.056 

0.0450 

0.0600 

0.6125 

0.057 

0.0300 

0.0600 

-0.5625 
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Figure 59. Comparison of Initial and Improved VSWR 
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APPENDIX A 

LIBRARY OF THE CHARACTERISTICS OF SELECTED MATCHING ELEMENTS 

The. following families of curves are included to 
provide engineering insight into the type of response 
associated with each element. Curves also indicate the 
dependence on parameters of the discontinuity. 

The families of curves provide the current reflection 
coefficient and the normalized shunt admittance as a func- 
tion of physical dimensions and frequency, for the capaci- 
tive iris, the inductive iris, the solid metal inductive 
post and the dielectric post located in an otherwise 
matched waveguide. A waveguide description and an equiva- 
lent circuit are given for each type of discontinuity. 

The following results are based on an evaluation of 

13 

the respective equations given by N. Marcuvitz, and 
reproduced in Chapter III, Equations (34) to (48). The 
following curves were calculated for propagation in rec- 
tangular waveguide with a width a = .9 inches and height 
b= . 4 inches . 

A typical discontinuity of each type was fabricated 
and its characteristics measured. A comparison of the 
evaluation of Equations (III-34) to (III-48) for the re- 
spective physical parameters. Chapter V, Table 2, shows 
good agreement with the results measured by the compari- 



son ref lectometer , as described in Chapter V. All curves 
given apply to the admittance Smith chart plot. 
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Equivalent Circuit 


(a) Capacitive Iris 



(b) L = .1 cm 


Figure 60. Capacitive Obstacle of Finite Thickness, Shunt Susceptance with a 
Fixed Thickness = .1 cm, Varying Iris Height d' in cm 
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(a) L = .075 cm 
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Figure 61. Capacitive Obstacle of Finite Thickness, Shunt Susceptance 

with a Fixed Thickness L, Varying Iris Height d' Centimeters 
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(a) Magnitude 
L = . 1 cm 


(b) Phase 

L = . 1 cm 


Figure 63. Magnitude and Phase of the Current Reflection Coefficient of a 

Capacitive Iris, Thickness of .1 cm Varying Height in Centimeters 
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(a) Magnitude (b) s Phase 

L = .15 cm L = .15 cm 

Figure 64. Magnitude and Phase of the Current Reflection Coefficient 
of a Capacitive Iris, Thickness of .15 cm Varying Height 
in Centimeters 
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(a) Asymmetrical Inductive Iris 


(b) L - .05 cm 


Figure 65. Inductive Iris of Finite Thickness, Shunt 

Suspectance with Fixed Thickness L. Varying Height 
d' in Centimeters ' 
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(a) L = .15 cm (b) L = .15 cm 
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Figure 66. Inductive Iris of Finite Thickness, Shunt Susceptance 
with Fixed Thickness L. Varying Height d'/2 in 
Centimeters 
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Figure 67. Magnitude and Phase of the Current Reflection Coefficient of a 

Inductive Iris, Thickness of .05 cm. Varying Height in Centimeters 
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x = • 20 
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x = .20 


Figure 72. Magnitude and Phase of the Current Reflection Coefficient of a Metal 
Inductive Post, Sidewall Distance x = .20 cm. Diameter Varies in 
Centimeters 
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(a) Magnitude 
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(b) Phase 
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Magnitude and Phase of the Current Reflection Coefficient of a Metal 
Inductive Post/ Sidewall Distance x = .50 cm. Diameter Varies in 
Centimeters 
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(a) Dielectric Post 


(b) Dielectric Constant = 4.07 
Diameter = .114cm 


Figure 75. Shunt Susceptance of a Dielectric Post with a Fixed 

Dielectric Constant, Fixed Diameter, Varying Sidewall Distance 
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Figure 76. Shunt Susceptance of a Dielectric Post with a Fixed 

Dielectric Constant, Fixed Diameter, Varying Sidewall Distance 
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(a) Dielectric Constant = .9 
Diameter = .328 cm 


(b) Dielectric Constant = 16.0 
Diameter = .114 cm 


Figure 78. Shunt Susceptance of a Dielectric Post with a Fixed 

Dielectric Constant, Fixed Diameter, Varying Sidewall Distance 
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(a) Dielectric Constant = 16.0 
Diameter =.228 cm 


(b) Dielectric Constant = 16.0 
Diameter = .328 cm 


Figure 79. Shunt Susceptance of a Dielectric Post with a Fixed 

Dielectric Constant, Fixed Diameter, Varying Sidewall Distance 
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(a) Magnitude 
d = . 114 cm 


(b) Phase 

d = .114 cm 


Figure 80. Magnitude and Phase of the Current Reflection Coefficient of a 

Dielectric Post with 3.85 Dielectric Constant, .114 cm Diameter, 
and Varied Sidewall Distances in Centimeters 
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(a) Magnitude 
d = .228 cm 


(b) Phase 
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Figure 81. Magnitude and Phase of the Current Reflection Coefficient of a 

Dielectric Post with 3.85 Dielectric Constant, .228 cm Diameter, 
and Varied Sidewall Distance in Centimeters 
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(a) Magnitude 
d = .343 


(b) Phase 
d = .343 


Figure 82. Magnitude and Phase of the Current Reflection Coefficient of a 

Dielectric Post with 3.85 Dielectric Constant, .343 cm Diameter, 
and Varied Sidewall Distances in Centimeters 
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(a) Magnitude 
d = .114 cm 


(b) Phase 

d = .114 cm 


Figure 83. Magnitude and Phase of the Current Reflection Coefficient of a 
Dielectric Post with 9.0 Dielectric Constant, .114 cm Diameter 
and Varied Sidewall Distances in Centimeters 
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(a) Magnitude (b) Phase 

d = .228 cam d = .228 cm 


Figure 84. Magnitude and Phase of the Current Reflection Coefficient of a 
Dielectric Post with 9.0 Dielectric Constant, .228 cm Diameter, 
and Varied Sidewall Distances in Centimeters 
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Figure 85. Magnitude and Phase of the Current Reflection Coefficient of a 
Dielectric Post with 9.0 Dielectric Constant, .343 cm Diameter, 
and Varied Sidewall Distances in Centimeters 
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(a) Magnitude 
d a .114 cm 


(b) Phase 

d = .114 cm 


Figure 86. Magnitude and Phase of the Current Reflection Coefficient of a 

Dielectric Post with 16.0 Dielectric Constant, .114 cm Diameter, 
and Varied Sidewall Distances in Centimeters 
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(a) Magnitude 
d = . 228 cm 


(b) Phase 
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Figure 87. Magnitude and Phase of the Current Reflection Coefficient of a 

Dielectric Post with 16.0 Dielectric Constant, .228 cm Diameter, 
and Varied Sidewall Distances in Centimeters 
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Figure 88. Magnitude and Phase of the Current Reflection Coefficient of a 

Dielectric Post with 16.0 Dielectric Constant, .343 cm Diameter, 
and Varied Sidewall Distances in Centimeters 
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COMPUTER PROGRAM LISTING 

Three computer programs are described on the follow- 
ing pages. The programs are written in Fortran V computer 
language for the Univac 1108 compiler. The first program 
described is referred to as Test 1. This program calculates 
the mean reflection coefficient generated as a function of 
distance. That is, it determines the electrical location 
of discontinuities in a waveguide. A functional block 
diagram of Test 1 and Test 2 are shown in Figures 89 and 
90 respectively. The second program, Test 2 calculates the 

complex reflection coefficient as a function of frequency at 

12 

a predetermined position in the waveguide. The above two 
programs utilize voltage measurements taken by the compari- 
son ref lectometer . 

The third program listed is the impedance matching 
program. This program uses data generated by Test 2, as the 
disturbance to be impedance matched, Chapter III. The match- 
ing program then calculates the physical dimensions and lo- 
cation of a predetermined matching element. The location is 
relative to the original disturbance. Generally, this pro- 
gram requires five seconds run time on the Univac 1108 to 
realize an impedance match. 

Comment statements have been provided throughout the 
program listing to partially explain the meaning of para- 
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meters used. 

Test 1 begins with a listing of constants to be used. 
Statement (9) and (10) give values to XMAX and XMIN, the 
range in distance along the waveguide relative to the ref- 
erence step. A table of sines and cosines are generated by 
statements (27) to (33). The use of this table reduces com- 
puter run time, since it requires less time to obtain this 
information from memory than to calculate the sine and co- 
sine when needed. An envelope shaping function,'*' when used, 
reduces side lobe levels at the expense of broadening the 
main peak of the locating curve. This part of the program 
has no direct application to this research, therefore, state- 
ment (38) prevents it from entering the following calcula- 
tions. The constant NO read in at statement (43) is the 
number of sets of data to be read by the computer. The con- 
stant NF is the number of data points measured. REFPT 
the magnitude of the reflection coefficient of the ref- 
erence step “used. 

Data is read in without a decimal point present. 
Therefore, statements (57) through (60) place the decimal 
point in the proper location. Statements (69) through (72) 
calculate Equation (III-31) . Statements (72) and (78) 
determine the wavenumber and the wavelength for each fre- 
quency corresponding to a voltage measurement by the compar- 
ison ref lectometer . Equation (III-31) is normalized to the 
average value of |r| (YRS(N)), for purposes of the computer 
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Figure 89. Functional Flow Chart of Test 1 Computer 
Program 
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calculation, by statements (80) to (90). The total mag- 
nitude of the reflection coefficient measured is plotted 
out by statements (91) and (120) using a plot subroutine 
not given here. The wavenumber range is adjusted by state- 
ments (121) to (159) in order to satisfy Equation (III-29). 

A single distance integration is performed by state- 

17 

ments (162) to (178). It has been found by D.L. Hollway 
that linear interpolation between values of YRS(N) cannot 
be used without loss of accuracy, therefore, integration is 
performed by considering each value of YRS(N) to apply over 
the internal from (WN(NN-l) + (WN(NN))/2 to (WN(NN) + 

(WN (NN+1) ) /2 . The range of integration WNRR extends from 
WNT (top) to WNB (bottom) . 

The discrete Fourier transform, or the fast Fourier 
transform, was not used in this program because the discrete 
Fourier transform requires that data samples be taken at 
equal increments of the argument. While samples are taken 
at equal intervals of frequency (50 MHz) , the transformation 
is from wavenumber to distance. The wavelength of the 
energy in the waveguide deviates from the free space wave- 
length due to the dispersion of energy within the wave- 
guide, therefore, the wavenumber intervals are not equal in 
increments . The wavenumber intervals corresponding to 
equal frequency intervals, increases as the frequency is 
increased. 

H(L,1) and H(L,2) are the real and imaginary parts 
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of the reflection coefficient generated at the distance 
Y (L) , statements (185) and (186). Depending on the slope 
of the magnitude of this reflection coefficient the dis- 
tance Y (L) is increased by . 5 or .25 by statements (197) 
to (202). If Y (L) is equal to or greater than XMAX in 
statement (204) the calculations are completed and the re- 
maining plots and printout are executed. 

Test 2 is very similar to Test 1, the difference 
being in the size in the internal WNRR. In Test 1 the en- 
tire spectrum measured was used for each distance calcula- 
tion. The result is a magnitude and phase of an average 
reflection coefficient generated at the distance of interest. 
In Test 2, however, a single distance or location is of 
interest and this is read into the program by XSET. For 
this particular distance XSET, the reflection coefficient 
is calculated over overlapping subintervals of 500 MHz. 

Each subinterval overlapping by 100 MHz, Figure 13. The 
value of the reflection coefficient calculated over any par- 
ticular subinterval is assigned to the center frequency of 
that sub interval. The result is that at approximately 100 
MHz intervals the reflection coefficient is known. Since 
each adjusted subinterval WNRR must satisfy Equation (IH- 
29), where the distance L n is the predetermined location, of 
the discontinuity, the resulting reflection coefficient is 
known at unequal subintervals of frequency. 

Calculations are actually made at ten locations along 
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Figure 90. Functional Flow Chart of Test 2 
Computer Program 
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the waveguide at intervals XINCR centered around the loca- 
tion, of the discontinuity, XSET. This is done to provide 
insight as to the effects of dispersion, on the reflection 
coefficient of the disturbance measured, as the energy tra- 
vels through the waveguide. 

In order to provide a smoother, more accurate curve, 
a linear interpolation is made between the 50 MHz points, 
resulting in data at 10 MHz intervals. This does not add 
any information about the total reflection coefficient, but 
it reduces the ripple in the resulting waveform, inherent 
in the algorithm used to perform the integration. 

The upper bound of the first subinterval is the upper 
bound of the measured spectrum, point NF. The lower bound 
of this subinterval is given by the statement (109) and is 
referred to as the NS point. The remaining calculations 
through statement (184) are duplicated from Test 1. The 
subinterval is shifted down the spectrum 100 MHz by the 
statement (216), and the program returns to statement 
Table (302) to determine the reflection coefficient from 
this new subinterval. When subscript NS is less than or 
equal to NSS., the subscript corresponding to the lowest data 
point taken, the program then exists from statement (110) to 
statement (220) . The program prints and plots results and 
ends . 

The impedance matching program calculates the phy- 
sical dimensions of a predetermined matching element and the 
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location of the element in the waveguide relative to the 
location of the disturbance as calculated by Test 1. The 
results of Test 2 on punched data cards are the input data 
for this matching program. 

The program begins with a list of constants. Comment 
statements have been provided in the program listing to iden- 
tify selected parameters. The matching element to be used is 
identified by the constant MACHEL, the number of data points 
by LE. Since the data points are not known at equal fre- 
quency increments, the frequency points are read in at state- 
ment number (60) . The real and imaginary part of the reflec- 
tion coefficient is read in corresponding to the highest 
frequency point first, by statements (62) and (63). The 
complex reflection coefficient of the test element is 
referred to as GAM(N), where N varies from 1 to LE. 

This program is equipped to begin calculations from 
any position along the waveguide, the starting position 
being denoted by ALENGH(l). The constant NOO read in at 
statement (108) denotes the number of times the results will 
be recalculated at corresponding starting points 
read in by statement (111) . This is a particularly interest- 
ing feature, since the set of physical parameters and the 
quality of the impedance match is a function of Xg/2 changes 
in the ALENGH ( 1 ) starting point. Chapter V. Generally, the 
best impedance match will be achieved by setting the starting 
ALENGH to zero, i.e. the electrical location of the test ele- 
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ment. 

Depending on the matching element selected, various 
appropriate starting points are programmed in statements 
(116) to (119) and described in Table 2. The three dimen- 
sional array is generated by statements (128) to (139) . A 
description of this dimensional space array is given in 
Chapter III. There are two cases possible, the matching ele 
ment may be located on the generator side of the measured 
mismatch or the matching element may be located on the load 
side of the measured mismatch. The results of Appendix C 
are used accordingly. 

Statement (149) branches the program to the appro- 
priate section of the program pertaining to the matching ele 
ment designated by the value of MACHEL. The reflection co- 
efficient for this chosen matching element is then cal- 
culated for the seven points of interest in dimension space, 
see Figure 19. The value of the complex reflection co- 
efficient is carried by the variable four dimensional array 
GGAN (I, J, KK, N) , where the first, second and fourth index 
correspond to the axis of the dimensional space array of 
Figure 19. The third index corresponding to frequency. 

When the array GGAM is completed, control is transferred to 
statement (375) . Statement (378) to (400) calculate the 
RMS of the mismatch, resulting from the measured discontin- 
uity and the matching element, for each point in dimension 
space. Statements (405) to (418) determine the smallest 


RMS value, while statements (419) to (432) determine the 
value of the I, J, K indices that correspond to this small- 
est value. Statement (436) tests for termination. If the 
indices corresponding to the least RMS are the 1=1, J=l, 

K=l, then a minimum value has been determined, if not, con- 
trol is shifted to statement (439). If a minimum value has 
been determined, a check is made to determine whether the 
increment size has been reduced to 1/8 of the original size, 
if not, control is shifted to label 350 where the present 
increment size is divided by two. If on the other hand the 
increment size is 1/8 of the original value, control is 
shifted to statement (447) or label 200, and the program be- 
gins to write output. If a minimum was not determined at 
statement (436) control if shifted to statement (439) or 
label 108 and the physical parameters of the dimension space 
(1, 1, 1) point are set equal to the physical dimensions of 
the point corresponding to the LRMS mismatch. Control is • 
then shifted to lable 250 and calculations are again per- 
formed for the seven points in dimension space centered 
around the new (1, 1, 1) point. 

Calculations to determine the physical parameters of 
the matching element are performed over the bandwidth of in- 
terest having the index points JBEG and JEND. When the match 
ing element size has been determined the results of the mis- 
match are calculated over the entire known spectrum origin- 
ally read in. Therefore, the characteristics of the chosen 
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matching element must again be calculated as determined by 
statement (468) , and the resulting mismatch again calculated. 

The output consists of printout of RMS mismatch over 
the band of interest, size and location of the matching ele- 
ment, and a listing of the reflection coefficient and the 
VSWR of the mismatch over the full band. A plot is made of 
the VSWR of the original disturbance and the matched dis- 
turbance over the full band. Also, data cards are punched, 
recording the real and imaginary parts of the resultant mis- 
match. Some typical results of this program are given in 
Chapter III. 
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Test 1 


Statement Label 
Number Number 


Statement 


*+ 2 * 

. 1 * 3 * 

<*<** 

_ 45 *_ 
4b* 

_4.7* 

48 * 

. 49 *. 
50 * 

....51* 

52 * 
53 * . 
54 * 

._. 55 *__ 

56 * 


DIMENSION I SUF ( 2000 0 ) ... ._ 

DIMENSION PT I TLE ( 12 ) »RTITLE<12> >MTlTLE<12) 

DIMENSION X ( 40 0 ) #VTITLE.<12> »HTITLE(12) 

DIMENSION P ( 200 > r 5 ( 1260) # OKS (100) i YREF ( 90 ) »YU(90> 

JIMENSIO.N-.r.RilQQ.j >WNaoi)).i.;vi^aQai-i-H.L4ll0.fJiXj.X1.40Ql-«FR(.lQQ?^r.SS.(10ilJ 

DIMENSION NREF ( 100) f NYU(IOO) 

...C PRESET . CONST ANTS- ANO - TABLES. OF-FUNC T I ONS * 

C XMAX IS THE MAXIMUM DISTANCE IN CMS FROM THE REF. REFLECTION 

XMAX-10Q • , 

XMIN=5. 

a 0 tl — 0 • 9 

PuRF-1 • 0 * ” ” 

-X PLPX PHASE ANGLE OF- THE .REFLECTION- COEFFICIENT-AGAINST-DISTAnCE* 

PLPX = 1. 

PLPX S-JL. 

NS- 1 

PI=3. 1415926536 .. 

C RPV f RADIANS PER VOLT. IS A CONSTANT OF THE RATIOMETER • IT IS THE 

_..C _ CQNDUC T I ON ANGLE AT FULL-SCALE -DI VIDED-BY— THE-- CORRESPONDING 

C VOLTAGE OUTPUT. 

— - .. RP V = 56. 6 *P I/. (6. .926* 180.-) 

C ATTEN CORRECTS FOR A WAVEGUIDE ATTENUATION OF 0.036 09/FOOT. 

ATTEN=0 . 036*2. 3026/ (10. *12. *2. 54) 

C TABLES OF SINES AND COSINES 

D06N- 1 * 500 

S(N)=SIN(N*?I/500.) 

S .( Nt5 0 U ) =rSXN ) 

6 CONTINUE 

J06N=1»26Q 

S(N*1000)=S<N) 

.6 CONTINUE 

C ENVELOPE SHAPING FUNCTION 

DO35N=Xf20Q ___ 

PANGLE=2, *PI*N/20Q. 

.. P ( N> =1 . -0 . 889*C0S (PANGLE) +0 . Oli2*COS (2. *PANGLE> 

P ( N ) - 1 . 

35. CONTINUE 

C 2 

C NO IS THE FIRST DATA CARD. 

HEAD ( 5r 601 ) NO 

601 FORMAT ( 1 3 ) 

READ (5.400) NF. 

400 FORMAT (13) 

READ l 5 f 12 ).R£FPI 

12 FORMAT ( F10 .5) 

REFPTS = REFPT*RErPT 

UOfaOOK-1 t NO 

READ (5r 20 0) <NR£F ( N ) # N=1 » NF.) 

RE AD ( 5 » 20 1 ) (NYU(N) »N=lrNF) 

2QQ. ..P.0RMAT-(.4 X r. I 4.I.XX f 1.4 iJJCrX4.iXX.».l-4 t.iXj I4U.X tJUU.lX «_1 tjxJJUX 4 <-lX xliUl X *.I U 

l»lX#I4*lX»I4»lX»I4tlX04#lX»l4»2X> 

Z01„ FORMAT. (4Xr 14 r IX r IA.UXlX 4' iXr.I.4_f 1X».I4 rlxa4».lXrI.4 r i>Cr 14 #XXi 14 >.1X r 14 

lrlX»I4.lX»I4.1XrI4#lX»I4.lX»I4.2X) 
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57* 00500N=1#NF 

58* _ YREFtN) =NREF(N)/1000.« _ 

59* YJ<N>=NYUlN)/1000. 

. 60* __ 500 CONTINUE ... 

61* w R I TE ( 6 # 502 ) 

.... 62* 502-f OR MAT ( 1H1 > 30X >.1R£FER£NC g-YOLT AGE-DATA.1a//J 

63* WRITE (6. 501) (YREF(N) *N=1 rNF) 

£>4* ..... . _ «RITE<6/503> 

65* 503 F0RMAT(////*30X»*V0LTAGE DATA WITH UNKNOWN IN LINE’’//) 

66* wRITE<t>* 501 ) ( YU<N) »N=1»NF) - .... 

t>7 * 501 FORMAT (15(3X»F5«3) ) 

h 6 * : B R ( i ) = 7 . 9 7 5 

69* i)0202N=l»NF 

70* ... CKS(N)= (REFPT**2*C0S(YREF(N) *RPV> ) /SlNt YREF ( N> *RPV>- 

71* Y.RS<N)-(CKS<N> *SIN(YU<N) *RPV) ) /COS < YU ( N ) *RPV ) 

72* YR(.N>=bORT( YRS(NJ) ) - - .... . 

73* FR(N+1)=FK(N) + «050 ' 

.... 7 9 * 202... CO. NT I NUE . 

76* nFP=NF+1 

76* U056.N=lrMFP 

77* wN(N>=50RT< < FR < N ) /29 . 9695 ) **2- (1 . / < 2. *2 . 54*WGW> ) **2 > 

. 7b* .. . «L (N)=l ./WN<N> . .: 

79* 56 CONTINUE 

— a a* yrw=o. 

til* YRSW=0. 

62* . _ u049N=1»NF .. 

63* YR*V = YRw*YR<N)*(WN(N + 1)-WN(N)) 

64* Y.R5W = YRSW*YRStN)*(WN(N + l)-WN(N) ) . - 

65* 49 CONTINUE 

— 66* YRSM = Yri5w/(W.N(NF.tJj^W.NiNSU 

67* YV*CYR,V/(^N(NF*1) - WN(NS) ) 

03*.. - D050N-1 1 nF - 

89* YRS(N)-YR5(N)"YRSM 

90* . 50 CONTINUE 

91* 1 F ( PLRF ) 16 * 16*51 

— 92*. 51-C0NTINUE . 

93* REAL LTR5Z. LENGTH 

94* LOGICAL 30Xf FlXLlMfUSELlH 

95* OATA HTITLE/ f FREQUENCY IN GHZ 

96* 2 - ‘ .. 

97* DATA VTITLl/'MAGNITUDE OF TOTAL REFLECTION COEFFICIENT 

9.6* 2 11 

99* CALL PLOTS( I9UFU ) >20000*2) 

100* D0401N = if NF 

101* x<n)=yr<n> 

102* 401 CONTINUE _ 

103* CALL PL0T( 10. *-3. »-3) 

.104* >!0DE=1 

105* LTRSZ-. 1 

105* UP = 2 ♦ 5 

107* UVER-0 . 

106* ... . HEIGHT-5. . ... 

109* LENGTHS. 

.110* _"_PE.RCiNT = #.9 

111* NJMIN T-b 

112* ML3N0=FR(I) 

113* HU6Ni)=M<<69> 


189 


114* I COUNT-89 

lib* ... dOX = .TRUEt ..... 

116* fixlim=. false. 

117* .. U5ELIM- • FALSE*... . 

118* 403 CALL PLOTX ( MODE* X 1 1 COUNT * UP * OVER * HEIGHT * LENGTH *80X . FI XLIM * USELIM » 

120* 16 CONTINJE 

121* X FUNCTION! OF DISTANCE*. ----- 

122* aNR=a\N|(NF)-WN<NS> 

123* £LDASH=0.5/WNR — 

124* IST = ELDASH*-1. 

_125 * i_=l 

126* YRS<NF+i>=YRS(NF> 

127* XF(XMlN-IST)86r86*84 

126* 84 r(L)=XMXN 

129* *0 TO aa — — _ 

130* 8b Y<L>=IST 

-131 * C — DISTANCE --LOOP. BEGINS 1 

132* 88 CONTINUE 

133* NN=NS . ... 

134* RM=(Y(L>M).04>/ELDASH 


136* XA=0. 

-137* X8 = 0. 

138* w.\lT=( WN(NF + l)+WN<NF) >/2. 

139* *'NRR=tfNR*M/RM - 

140* C WNR- TOTAL WAVE NUMBER RANGE. 

141* . C WNRR = A FUNCTION OF DISTANCE#^ THE RANGE REDUCED TO MAKE. THE ANGLE AN . INTEGER 

142* C *2*PI 

-1 4 3 * a N3=W N T - wN R R 

144* ATTER=1.*ATTEN*Y(L) 

145* .. 85 lF(WN(NN + l)+WNCNN)-2.*l»!NB)8lr9lf 95 

146* 81 NN=NN+1 

147* GO TO 85 .. . 

148* 95 IF <NN-l>301r301f89 

-14 9* 3ni_.NN = NN+i 

150* GO TO 91 

151* 89 IF (WN(Nn)+WN(NN-1)-2.*WN3>- 91*91*93. _ - . 

152* 93 ivN=NN-l 

153* . . GO TO 95 

154* 91 CONTINUE 

-135 *— * NA = I WN ( m ) ± W3I t NN.t.ll. .) JZ 

156* w N S 0 = W N ( N N * 1 ) 

157* wNST = WN(NN) - .. 

158* «NINN)=Wn3*WNA/2. 

159* aN(NN-1)=WN3-WNA/2. — .... 

160* C 6 

-161* _jCL£I&Gl£ D I S T A NC £_IN I££RAZXjQN 

162* U094N=NN*NF 

163* .. ALP=2.*(Y(L>-0.G4)*WN<N) - ....... _ 

164* C ALPHA=4,*PI*Y(L)*WN(N) # NALP IS THE MULTIPLE OF 2.*PI REMOVED. 

165* NALP=ALP _ 

166* AA=1000.*(ALP~NALP> 

^167* XF_- CAA) ...90 

166* 90 A A-A A *1 0 00 • 

169* . 92 NA=AA 

170* AANA=(AA-NA) 
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171* SALP=S<NA)*(S<NA+1)-S<NA> >*AANA 

172* . NASNA+250 - 

173* CALP=S(NA) MS(NA+i)-S<NA) ) *AANA 

174* IP = i. + UN(N>-WN(NN) )*199*/(WN<NF)-WN(NN>) . 

175* YRSPI = YRS<N>*P<IP>*<WN(N*U-tfN(N-in 

17b* XA = XA + YRSPI*CALP_ 

177* X3=XB+YRSPI*SALP 

17S* 94 CONTINUE - - — 

179* WN(N.N) =«NST 

180* hN<NN~l>=WMSO - - 

181* HL1=XA/(REFPT*WNRR*2. > 

182* HL2=X3/(REPPX*rtNRR*a.l 

163* HL50=HL1*HL1+HL2*HL2 

184* CORRN=ATTER/{l.-HLSOrREFPTS+HLSO*REFPTS>~ 

185* rt f L * 1 > -HLl *CORRN 

166* H(L*2>=HL2*C0RR.N . ... - 

187* H(Lf 3)=SORT<H(Lrl>**2.+H<L»2>**2.) 

.188* ANG = 57,29b7a*AT.ANCH(L-r^JL/HCLjJl)J 

189* IF <H(L*1>> 96 * 102 * 102 

190*. 96 lr (H<L*2>) 100*98*98 

191* 93 rl(L*4>=ANG + 180. 

192* _ GO. TO 104 . . 

193* 100 H(L*4)=ANG-180. 

194* GO-TO 104 

195* 1 02 H(L#4)=A:NG 

196* 104 CONTINUE -- 

197* NXL=Y(L) 

196* IF ( Y(L)-NXL-O.Gl) -U2»112*U6 - - 

199* 112 IF (H<L*3>-0.005) 114*114*116 

20Q* 114. Y (L+1>=Y(L) *Q»5 

201* GO TO 118 

20 2* 116 Y ( L+ 1 ) -Y ( L) *0 • 25 - 

203* 118 L=L + 1 

204* IF <XMAX-Y(L)> 122*122*88 - 

205* C 7 

20o* IL_DXSJ.ANCE=LOOR-EUDS - 

207* 122 CONTINUE 

208* LE-L-i . 

209* wRITE(6*109> 

210* 109 FORNlATIlHlr ' Y ( Ll • *6X * Y A ( N ) lr 6X * 'B(N).! #3X* * MAG* REFP-T • -COEFF « 1 #3X*. 

211* 1 ' PHASE ANG* * * /// ) 

-2 1 2 * 1Q6_W RITElbdQ 8 U Y < L ) * (H ( L> J) ..^s ijA ) iLZl » LE) ... 

213* 108 FORMAT(F1Q.2*3F10.5*F10.1> 

-214* 110. CONTINUE 

215* IF (PLRF) 411*411*412 

216* 412 CONTINUE - 

217* DATA RTITLE/* DISTANCE IN CENTIMETERS 

218 * 2. .♦/. 

219* DATA PT I TLE/ 1 MAGNITUDE OF AVERAGE REFLECTION COEFFICIENT 

220* 2 . _*/ 

221* U0404N=1*LE 

222* ... .. X(N)=H(N*3) 

223* 404 CONTINUE 

.224* _AR = X (.1 ) 

225* j01000N=1*LE 

226* _ IF(X(N>.ST.AR)AR=X<N) ... - 

227* 1000 CONTINUE 
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226 * 

VLBND- 0 . 0 


229* . 

V J3ND=AR* . i*AR 


230* 

CALL PLOT ( 10 • » -3. .-3) 


... 231* . 

MODE=l 


232* 

LTRS2=.l 


233* 

uP -3 . 


234* 

OVER-O . 


... 235* . _ 

LENGTH=5 . 


236* 

HEIGHT=5,0 



237* HL8ND=Y(l) 

233* HJ3ND=Y<LE) 

_239* NJMINT = 12 

240* ICOUNT=LE 

241* ._ .. 60X=.TRUE. . . * _ .. . . 

242* FIXLIM=. FALSE. 

243* USELIM=. FALSE. .... ._ .. . 

244* fixlim=.trje. 

-245* uSELIM=.TRUE. 1 

246* PERCMT- • 99 

247* CALL PLOTY (MODE » X . Y # I COUNT . UP * OVER . HEIGHT * LENGTH » BOX 

24 6* 1 #NUMlNT/PnTLE^TITLt/LTRSZrVL8NOrVU3NO#HL8NO*HU8ND 

249* 411 CONTINUE 

250 * IF (PLPX) 408»408»410 

-251 * 41CU.C 0 : 4 T I NUE — 

252 * uATA MTITLE/ 'PHASE ANGLE 

253 * 2 •/ 

254* n04U6N=lrLE 

256* a<N)=H(N>4) 

256* 406 CGNTINUE 

-257* MODE- l 

256* CALL PL0T(10.*-3.»-3) 

259* hEIGHT"b. 

260 * LENGTH=b. 

261 * l.TRSZ=.1 . 

262* uP”3 . 

-263* UVER-Q , 

264* HL3N3=XMIN 

265* HJBND=Y(LE) 

266* IC0JNT=L 

267* uox- • True • 

266* FlXLlMz.TRUE. 

-2u9* LiSELl M- . TRUE* 

270* VJ3ND=130. 

271* .. VL3ND=-180. ' 

272* PERCNT = . 99 

273*. . _ CALL PLOTY (MODE.Xr Y r I COUNT . UP r OVER # HEIGHT * LENGTH. 30 X 

27 4* 1 . MUMInT . MTITLE»RTITLE»LTRSZ# VLBND » VUBND » HLBND » HU3N0 

-275* .403. CONTINUE 

27b* 600 CONTINUE 

277* . CALL PLOT ( 0 • 0 r 0*0*999) 

273* END 


* FI XL I M # l 
fPERCNT) 


. FIXLIM . 
* PERCNT ) 


U SELIM 


END OF compilation: 


NO DIAGNOSTICS. 


Test 2 
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Statement Label 

Number Number Statement 

J* DIMENSION PTITLEU2) 

2* dimension YL<20>»B<iO»50»4).p<200>.ENN<6> 

_3* DI.MEN519N_FM.I0J!) 

4* DIMENSION rllOO) 

- .5* DIMENSION IBUF ( 20000 ) > NREFilOO.) .'NYU ( 1Q0 )j.YREF.I IQjO LtJfU.i IOQJ-lSJ 1260) 

6* 1* CHS (100) »FR(500> »YRS<500> * YR< 100) »WN«500> *WL<SOO> »H<100.4> #X<100) 

7* . 2.HTITLE(12)fVTITL£(12)»RTITLE<12) 

6* DIMENSION ARYR(500> 

-J9* CALL PLOTSdRU Ftn tanonn.gi 

10* MSM=0.9 

u* _. . ns-i : 

12* PI=3. 1415926536 

13* RPV = 56,&*PI/<6.92B*160.«J 

14* ATTEN=0.036*2.3026/<10.*12.*2,54) 

.15* D06N = 1 r 500 

16* S(N)=SIN(N*PI/500.) 

17* .. . S(N+500)=-S(N) 

16* 6 CONTINUE 

19*. .. D08N=1»260 

20* S<N+1000)=S(N> 

.21* ^-CONTINUE 

22* READ(5.2000)NO 

23* 2000 F0RMATII2) 

24* READ ( 5. 12) REFPT 

25* 12 FORMAT (F10 .5) 

26* C XINCR IS THE INCREMENT IN CM. FROM THE LOCATION OF THE MISMATCH. 

.27* READ (5r 601 )XINC& 

26* 601 FORMAT ( F 10 • 5 ) 

29* REFPTS=REFPT*REFP.T 

30* D02001 1 J=1 > NO 

..31* C — NF IS THE UPPER . FRE8UENCY..UMIX. 

32* READ ( 5. 9 ) NF 

-33* 9 FORMAT 113) 

34* C XSET IS THE LOCATION IN CM. OF THE MISMATCH FROM THE REFERENCE. 

.35* READ ( 5’ 602) XSET 

36* 602 F ORMAT ( Fl 0 . 5) 

. 37* READ 1 5> 200 ) (NREF ( N). f N = 1 .NF) 

36* READI5.201) (NYUtN) *N=l.NF) 

-39.* 200..F ORMAT 1.4 X i X4. IX. Tu.iy.Tu.i*. fit. iy.T4.ly.m.»y.T4.iy. fl uty. fa. ly.fq 

“0* 1»1X»14»1X»I4»1X.I4»1X»I4»1X*I4»2X) 

‘*1* 201 F0RMAT(4X»I4.1XiI4».lX_*I4»lX»14»lXa l UiXA.14^J.X».Iiua)U14#.lxai»axa4. 

92* 1.1X.I4.1X.I4.1X.I4'1X'I4.1X'I4.2X) 

.43* D0500N=1.NF . . . 

44* YR£FtN)=NREF(N)/1000. 

_45* YU <N ) =NYU<N )/100Q. 

46* 500 CONTINUE 

.47* .. FR(1)=7,975 

40* D055N=1»200 


49* 

50* 

51* 

52* 

53* 

54* 

55* 

_ 56* 

57* 

. 50* „ 
59* 

-JM* 

61* 

62* 

63* 


PANGLE=2.*PI*N/200 

P<N>rl*-0t£ 

P(N)=1. 

55 CONTINUE.. 


D02Q2N=1*NF 

US ( N ) = ( REFEI* *2*CQS ( Y RFF ( hi) » RP-V 1 ) /SI N ( YREF ( N) * R Ptf ) 

YRS(N)=(CKS(N)*SIN<YU<N)*RPV) )/COS(YU<N)*RPV) 

YR ( N) = SQRT ( YRS(N) ) 

202 CONTINUE 

. . NFP=NFU 

D01 002N = 1 r 461 

£R * N.tl l=FRlNl±.J)UQ 

WN(N)=SQRT ( <FR(N)/29.9695>**2-(l./<2.*2.54*WGW) )**2) 

WL(N)=1./WN(N> . ... 

1002 CONTINUE 


. 64* . jr i 

65* 00100QN=1*8B 

_66* DELY = YRS<Ntil^iaSlHl 

67* 0010011=1*5 

68* ARYR(J)=YRS(N) +DELY.*«2*(I*lJ__ 

69* J=J*1 

. 70* 1001 CONTINUE 

71* 1000 CONTINUE 

—72* ARYR< J).= YRS (fl9i 

73* NFF=J 

74* .... NF=NFF 

75* NFP=NFF+1 

76* D01004N=1*NFP 

77* YRS(N>=ARYR(N> 

_78* 100.4 CONTINUE 

79* YRSW=0. 

SO* D049N = NS * NF 

81* YRSW=YRSW+YRStN)*(Wt(N+l)-WL(N) ) 

82* 49 CONTINUE 

83* YRSM=YRSw/(*LtNFU>-WL<NS> ) 

_D05QN=NS*NF 

85* YRS(N)=YRS(N)-YRSVI 

86* 50 CONTINUE 

87* YRS(NF+I)=YRS(NF) 

88* _ YL(1)=X5ET-4.*XINCR 

89* NFF=NF 

-90* NSS=NS — 

91* CALL PLOT<0.0*3.0*-3> 

-92* D05000LL-1 * 10 

93* C PLOT AX 1 5 • 

94* _ CALL PL0TU0.Q*0.0^il 

95* CALL PL0T<0. 0*2. 203*3) 

—9 6 * _ C AL L_ .P LO T C 4...4 0 6 *_2 .203_t2J 

97* CALL PL0T(2. 203*0. 0*3) 

98* CALL PL0T<2. 203*4. 406*21 

99* CALL PLOT (0.0»0.0*3) 

100* YL(LL*1)=YL(LL)4XINCR 

101* *RITE(6*04) YL(LL) 

102* 84_.FQRMAT ( 1H0 » 1QX** DISTANCE* LEh+Z±l CM.*) 

103* ATTFR=i.*ATTEN*YL<LL) 

104* „. ._ XINT-10Q . 

105* lFCY(LL) .GT.21. ) XlNT=bO. 
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106* NF=NFF 

107* . .. L-l 

108* C BEGINNING OF NEW WAVENUMBER RANGE IDENTIFIED BY CL) 

_109*.„ __ 302 NS=NF-XINT 

110* IF CNSS+l-NS) 303001 #301 

__i 1 1 ♦ 30 3_C ON T I NUE 

112* NN-N5 

113* wNR = (WN(NF+l)+WN(NF)-WN.(Nai^ttN.tNS=i.UZ 2 * ; 

114* ELDASH=0.5/WNR 

115* RM=<YL<LL)-0.04)/.ELDASH 

116* M=RM 

—111* 1 E(M) 3 flliJSQij_afr 

118* 86 XA=0. ■ 

. 119* ... .. X3 = 0. ..... 

120* WNRR=WNR*M/RM 

121* WNT=tWN(NF + l) +WN(NF).)/2. 

122* C WNR=TOTAL WAVE NUMBER RANGE. 

-123* C_ WNRR-A, FUNCTION-OP— QISTAUCEi IS T HE RANGE REDUCED TO M AKE T HE ANGLE AN ■IN T EGER. 

124* C *2*PI 

125* * N3 = WNT* WNRR — — ‘ . 

126* 85 1F(WN(NN+1>+WN<NN)-2»*WNB)B1’91'95 

127* 81 NN=NN+1 . - ___ 

128* GO TO 85 

-129* 95 JF.(NN-1)3QL*-3Q1jA9 

130* 89 IF (WN(NN>+WN(NN-1>-2.*WNB> 9l#9li93 

131* ._ 93 NN=NN-1 __ _ 

132* GO TO 95 

133* ... 91 CONTINUE 

134* *NA=(WNtNN)+WN<NN+l) J/2.-WNB 

—135* wNSO=WNlNN-JU 

136* ivNST=WN(NN> 

137* _ wN(NN)=WNB+WNA/2. 

138* WN CNN-1 ) =WNB-WNA/2. 

139* ENN(1)=WNT 

140* ENN(2)=*NB 

—141* £NN..(3)=WNA 

142* C 6 

.. 143*. D094N=NN»NF .1 

144* ALP=2.*(YL(LL>-0.04)*WN(N) 

-145* NALP = ALP 

146* AA=1000.*( ALP-NALP) 

—147* 1FCAA)222>222>223 

148* 222 AA-AA+1 000# 

..149* 223 NA = AA . 

150* AANA=(AA-NA) 

- 1 51* SALP=5 ( N A ) 1 1 S C NA + ll-5 ( NAU.tAANA 

152* NA=NA+250 

—153* CALP=S(NA) + (SCNA*JJ_^SlNAlJ_*AANA 

154* 1P=1.+(WN(N)-WN(NN) ) *199 . / t WN <NF) -WN ( NN ) ) 

155* YR5PI = TRSCN)*P(lPJ^tfhL(Ni:ll^WNli^lJJ 

156* ENN ( 4 ) = X A 

157* ENN(5)=XB 

158* XA=XA+YRSPI*CALP 

—159* XB=XB + YRSPT.*SALP — - 

160* 94 CONTINUE 

—161* WN(NN)=WNST 

162* WNCNN-1)=WNS0 


163* 

164* 

HL1=XA/(REFPT*WNRR*2.) 
HL2-XB/ ( REFPT*WNRR*2. ) 

165* 

HLS0=HL1**2*HL2**2 

166* 

CORRN=ATTFR/(l.-HLSQ-REFPTS*HlSO*REFPTS) 

167* 

_168* 

C THE FOLLOWING NEG. SIGN CONVERTS FROM IMPEDANCE SMITH CHART TO 
C ADMITTANCE SMITH CHART. 

169* 

170* 

H(Lf 1)=-HL1*C0RRN 
H(Lf 2)=-HL?*CORRN 

171* 

172* 

B(LL>L.3)=SQRT(H«L*1)**2*H«L»2)**2> 
C PREPARATION for THE SMITH CHART PLOT. 

173* 

_174* 

ANG=ABS(ATAN(H(L»2)/H<L*1) )) 
IF ( H ( L # 2 ) >700.701.701 

175* 
176* 
177* 
178* 

701 IF ( H ( L * 1 ) )703»702*702 

... 702 ANGsANG 

SO TO 710 
703 AN5-PI-AMG 

179* 

J60* 

GO TO 710 

700 IF ( H ( L# 1 ) )70Uf705#7n*5 

181* 

704 ANG=ANG+PI 

. 182* 

. _ GO TO 710 ..... _ 


163* 705 ANG=2.*PI~ANG 

164* 710 CONTINUE 

165* 8(LL»L»4)=ANG*57.29578 

-166* FN ( L > =29 , 9695*£QRlXltfN (NF ) tV HRSUZ * ) * » Z±Ll^ALS 


187* 

168* 

189* 

190* 

191* 

192* 

X <L> =2.203+ (H<L» 11*2.203/. 22) 

Y<L>=2.203+<H<L.2)*2.203/.22) .... 

C PLOT CURVE WITH 9.10.11.12 GHZ MARKERS. 

IF(L-1 >950.950.951 

950 CALL PLOT(Xtl) * Y < 1 ) *3) 

951 CONTINUE _ _. _ 

193* 

AL=ABS(FN(L>-9.001) 

194* 

IF(AL-*05)507/507*512 

195* 

512 AL=A3S(FN(L)-10.001) 

196* 

IF< AL'.o^i 5nAf5nn. 513 

197* 

513 AL=ABS(FN(L)-11.001) 

..198* 

IF(AL-.05)509f 5n9f51U 

199* 

514 AL=ABS<FN<L>-12.001> 

200* . 

IFIAL-.05) 5ln#5ln» 50ft 

201* 

507 CALL SYMBOL<X<L)*Y<L>*.07#1»0#0»“1) x 

202* 

CALL SYMBOL(X(L) f Y (L) 9 ill 8H__9* 0_ GH£ #0.0# A) 

203* 

GO TO 506 

-204* 

5Q8 CALL SY^AOI (X(t ) • Y (L ) # * 07# } . 0 - fl # \ ) _ 

205* 

CALL SYMBOL<X<L> *TtL> • .1»9H 10.0 GHZ.0.0.9) 

206* 

GO TO 506 _ 

207* 

509 CALL SYMBOL(X(L) » Y(L) » . 07# 1 # 0 .0 # -1 > 

208* 

CALL SYMROL(X(L) # Y(L> # .1# 9H ll.Q GHZ*0.0.9) 

209* 

GO TO 506 

.210* 

_ 510 CALL S YMAQL (l(L) i Y{[) » a 07# 1 “1 ) 

211* 

CALL SYMBOL(XtL) #Y<L> * .1»9H 12.0 GHZ»0.0#9) 

212* 

506 CONTI NUF 

213* 

CALL PL0T(X<L>.Y(L>»2> 

214* 

C END OF SMITH CHART PLOT. __ 

215* 

LE=L 

216* 

NF-NF* 1 0 

217* 

L=L + 1 

216* 

C END OF REFLECTION. COEFFICIENT. AS. A FUNCTION OF FREQ. AT DIST.YL 

219* 

GO TO 302 
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220* 301 CONTINUE 

221* I F ( LL • NE • 5 ) GO TO 5000 __ 

222* WR I TE ( 1 » 1 071 ) LE 

. 223*_ . 1071 FORMAT ( 12) 

22*1* *RITEtl» 1070) (FN<I> » I=LE*lr-l) 

225* 1,070. FORMAT ( BF10 * 4 ) 

226* #RITE(l#1072MH(I#l)#IsirLE) 

227* WRITE C 1 r 1072) <HUr2) rI=lrLEl 

228* 1072 F ORMAT ( 8F1 0 • 6 ) 

229* 5000 CONTINUE .. 

230* wRITE<6>952) 

— 2.31* 952.-F0RMATlAHli_l.QX r * FR EQUENCY f r lQXiA-MAfiNlTUQS- OF REFLECTION COFFFICIEN 

232* IT * r /// ) 

. 233* D01010L=1*LE 

234* WRITE (6*953) (FNtL) # <B<LL#Lf 3) rLL=l*l0>> 

235* 953 F0RMAT(10XrF6.3r 10Xrl0F10.4) 

236* 1010 CONTINUE 

237* rtRITE ( 6 r.96Q ) 

238* 960 F0RMAT(1H1*10X»* FREQUENCY * # lOXr ’PHASE ANGLE**///) 

239* . D01020L = lrLE 

240* WRITE (6* 1021) (FN(L) » ( 9 ( LL# L* 4) r LL=1 r 10) ) 

. 241* ,1021 FORMAT (10X r F&« 3i.lOX# 10F-10*4J 

242* 1020 CONTINUE 

_243* 00229 J"1.#.LE 

244* I=LE+1-J 

245* ....... X(J)=B(5rIr3) 

246* 229 CONTINUE 

. 247* 00231 J = 1 f LE . ... 

248* I=LE*1-J 


..-.249* Y(J)=FN(I). ... : 

250* 231 CONTINUE 

.251* REAL LTRSZ t LENGTH 

252* LOGICAL BOX * F I XLI M * USELI M 

253*.. DATA HTITLE/ • FREQUENCY. INU6HZ. 

254* 2 */ 

—255 * 0 ATA__V.T 1TLE/ 1 MA SMI.TUQE nr RFFLFCTIilN rnFFPTr.TFNT 

256* 2 */ 

257* CALL PL0T<1Q.jl-3.jl=*>. 

258* M0DE=1 

259* LTRSZ=.l 

260* UP=4.5 

— 261* DVER = 0. 

262* HEIGHT=5. 

263* LENGTH=5. . 

264* HL3N0=Y(1) 

265* HU3N0=Y(LE) 

266* ICOUNT=LE 

_ 267* BOX = .TRU£. : 

268* USELIM=.FALSE. 

_ 269* FIXLIM=,FALSE* 

270* USELIM=.TRUE. 

... 271* FIXLlM=tTRU£. : 

272* VLBND=0. 

__273* VUBND= t2 

274* NUM I NT- 16 

...275* P£RCNT= * 9 

276* CALL PLOTY ( MODE# X*Yr I COUNT » UP# OVER r HEIGHT* LENGTH* BOX *FIXLIM* USELIM 
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277* 
_ 273* 

1 »NUMlNT * VTITLE»HTITLE»LTRSZ» VLBND»VUBND#HLBN0>HUBN0«PERCNT) 
D0230J=lrLE 

279* 

1=LE+1-J 

. 280* 

_.. X(J)=B<5f If4) 

28i* 

230 CONTINUE 

282* 

— _.CALL_fl.OJLilO.»-3.»-3> 

283* 

DATA HT I TLE/* FREQUENCY IN GHZ. 

.284* 

2 • / 

285* 

DATA RTITLE/’ANGIE IN DEGREES 

286* 

2 •/ ._ .. ... 

287* 

MODE-1 

_268* 

LTRS2=.i 


2B9* UP-4 > 5 

290* . OVER=0. 


291* " HEI5HT=5. ' 

-292* LENSTH=5._. 


293* 

294* 

HU8N0=Y<LE> 

HL9ND=Y(1) ..... ..... ... .. 

295* 

PERCNT=.9 

296* 

ICOUMT=LE 

297* 

298* 

299* 

_3QQ* 

BOX=.TRUE. 

USELIM=. FALSER - 

FIXLIM=. FALSE. 

VJ3UD=1B0» 

301* 

.302* 

VLBND=-180. 

NUMINT=16 _ 

303* 

304* 

305* 

306* 

CALL PLOTY ( NODE* Xf Y i I COUNT * UP » OVER* HEIGHT «LENGTH*BOX*FIXLIM*USELIM 

1 »NUMINT fRTITLEf HT.ITLE»LTHS2«VLBND«-VU3ND»HLBND»HUBN0»PERCNT1- ..... 

1006 D01007LL=1f10 

D010Q8L-1 »LF 

307* 

J=LE+1-L 

308* 

.... B(LL#Lrl)=8(LLf J*4) 

309* 

310* 

311* 

_312* 

1008 CONTINUE 

.. 1007 CONTINUE 

CALL PLOT (10. *»3* *-3) 

DATA PIITLE/'DISTANCE IN Cm- 

313* 

314* 

315* 

2 •/ 

LTRSZ=.l .... ..... _ ... 

UP=4.5 

316* 

317* 

_31 8* 

OVER=0 ..... _ ... 

HEIGHT=5. 

LE.MGTH=5. ....... 

319* 

320* 

HLBND=YL<1> 

. _ HU3ND=YL(10) 

321* 

322* 

PERCNT=.9 

...... ICOUNT=10 ........ 

323* 

80X=.TRUE. 

_324* 

__ ... USEL1M- • TRUE . 

325* 

_ 326* 

FIXLIM=.TRUE. • 

. ... VUSND=360. 

327* 

328* 

VLBND=0.0 

. ... . . NUMINT-6 .. ....... 

329* 

_330* 

□010111=1*4 
M00E=1 _ . 

331* 

332* 

IF(I.GT.1)M0DE=4 
F1=9,0tl-1 

333* 

LF = 5 
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334* IF(I.GT.1>LF=12 

335* I F < I * GT * 2 > LF = 16 

336* I F ( I . GT • 3 ) LF=23 

337* U01012LL=lrl0 

338* X<LL>=B<LL#LF»1) 

...339* 1012 CONTINUE 

340* CALL PLOTX(MODE*X> ICOUNT# UP » OVER* HEIGHT 'LENGTH #90X> FI XLIM#USELIM» 

341* . 1 NUMINT t RTITLE f PTLTLErLT-RSZi-VLBNDi VUBND rHLBNOfHUBNO* PERCNXJ 

342* 1011 CONTINUE 

.343* 600 CONTINUE 

344* 2001 CONTINUE 

34 5 * C ALL— PLOXIGaHj Q » 0 1 9.99 j 

346* ENO 


END OF COMPILATION; “ " NO DIAGNOSTICS* 



Matching Program 


Statement Label 
Number Number 


Statement 


1 * 


2* 

_ 3 * 

4 *“ 

5 * 

““ b*“ 

7 *__ 

8 * 

9 * 

10 *“ 

_ 11 *_ 

12 * 

_ 13 * 

14 * 

lb*_ 

16 * 

17 * 


DIMENSION D(100) 

HITMENS 1 0 OEXOW1 

DIMENSION FR(90) #WN(90) #WL(90) fFSWL(90) »DP(4) fWlDTH(4) ,ALENGH(4) 

' DIMENSION AA (90) V COT90TTSTT9D)TARMS (3> 3» 3T »SET< 90 ) »BETA ( 90 ) 

DIMENSION XREL(90) *YIMAG(90) 

DIMENSION lBUF(20 000),x?90)VYf90rfHTlTLE(12 ) >VtITLE( 12) 

DIMENSION ALPHA ( 3 > 90 ) , BETE ( 3 » 90 ) • BES0A ( 3 , 90 ) # BESOB ( 3 * 90 ) » BESl A ( 3* 9 

1 0 J * ’3Esi 8 < 3* 9n V 1 50PT3#’90 F7SI I'CZ » 90 J * 

COMPLEX gam ( 90) * AGAM(4,90) » GAMMA (90) • GAMMAH* BGAM » Yl # Z \ I t YT 

“"COMPLEX GG AM ( 3 1 3 V 90 *3)' 

COMPLEX Zl,ZJ 

"“COMPLEX X3AM ( 100) “ 

COMPLEX YII 


TTHTs PROGRAM DEtErM I NES~THE — mMEJJ5TtfN5 '0F~THE' - DE5TGNATED“ 

CALL PLOTS ( I 3UF (1)#200q0»2) 


“Matching element?* in order to “imped ancet~match the measured"" 

MISMATCH, 
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c_ list of constants. 

“ PI = 3. 1415926536' ~ 

C=2.996»(10**10) 

~C" waveguide width in inches=wgw 

wGy; = 0 , 9 
A - w G W* 2 . 5 4 

_C__ WAVEGUIDE HEIGHT IN _inches=wgh l 

B5wGH*2 . 54 

"C "E IS THE BASE OF' THE'NATUrAL' LOG". 

E=2. 7182818284 

5Q"f SOssqRTj 275 

C*******«»*************************************** **************** ****** 
c** *************** ************************* ***********•*****'•********** 

C** THE VALUE OF THE FIRST DATA CARD IS MACHEL, * 

"C** THE VALUE OF MACHEL' SELECTS THE' TYPE OF' MATCHING "ELEMENT; * 

C** THE VALUE OF MACHEL MAT RE 1 OR 2 OR 3 OR 4 WHERE * 

■'CVi'y'CORREsPONDs' TO'THE CAPACITIVE "IRISY ’ * 

C** 2 CORRESPONDS TO THE ASYMMETRICAL IND. IRIS. * 

"c** 3 corresponds to the solid inductive post; * 

C** 4 CORRESPONDS to the dielectric post, _ * 

C* ***•*,,,*,*•«***.*,**«** **,««**.*.**+«*****•*»**«*«•«**+***•********* 
C««*. *«»,<»**«.«,,«,.**.****,«****»*********,***«*«•**•******** **«•**«* 

~c data-card"order; : 

c first data card=machel 

"C" second data card=eprim 

c third data card=le(no of data points) 

~C~ FREQ, data cards • — 

C real and IMAGINARY parts of unmatched reflection coeff. 

"C BaND'LTMITS "To ' BE "MATCHED-OVER; 

C N00" THE number OF times THE DEVICE Is TO BE MATCHED 

~C " AlENGh the STARTING VALUE "FOR'LrNGTH'DOWN'THE"GUIDE; 

C NOTE" The number of alengh data cards shoulo BE = TO THE NO.'NOO 

READ(5»560) MACHEL ~ 

56 p FOrmAT(II) u l 

REaD(5i601)EPRIM 

"6O1 FORM*T(F10.3) " " ' 

READ 1 5» 1010) LE 

1010 FORMAT (12) "" 

WRITE(6.2050)MACHEL>LE 

~SV 5o _ FO R M A T i 1 H 1 , 7 ? 7 H 0 X i • M A'cH EL = ' , 1 1 > 10X » » NUM BER O F OA T A PdI N T S"= ' . 13 
READ<5#1011 ) (FR(I)»I=i,LE) 

1011 F0rmAT(8F10.4) ' 

REaD<5,1)(XREL(I)»I=1,lE> 

READ(Srl) (YIMAG(I) ?IS1,LE) 

1 FORMAT(8F10.4) 

WRi TECS* 2052)' 

2052 FORMAT!///, 10X, 'REAL PART OF THE MISMATCH',///) 

WRITE(6,205l) (XREL(I)7I5LE* 17«n 

__205l FOrmAT!10X,8F10.7) 

WRITE! 6,2053) 

2053 FORMAT!///, 10X, 'IMAGINARY PART OF THE MISMATCH',///) 

~WR I TE f 6 ,2051 V! YI MA G( I m =LE, 1,-11 

D08N-1 , LE 

GA M !N) =X R£LfNlTfO', ir*YIMA3fN) 

8 CONTINUE 
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75* 

”76* 

77* 

76 * 

79*_ 

SO* 

81* 

” 82 *" 

S3* 

84* 

85*^ 

86 * 

87* 

88* 

89* 

90* 

91* 

92* 

93* 

94* 

9b* 

' 96* 
97* 
98* 
99* 
'100* 
101 * 
102 * 
103* 
“104* 
105* 
“106* 
107* 
108* 
1 09* 

Ho * 


001020N=1*LE 
I=LE+1-N 
XG6M(N>=GAM<I) 
X02F CONTINUE 

OOl 02 1N=1 t LE 

^AMiMr=x^AvrmT 


1021 CONTINUE 

C~ FrBEG AND ~FREnD“ARE“ THE FrEOUENCY”BANO "LIMITS "IN' GHZ. ~0P THE“WATtHtN<5 

c bandwidth of interest, 

READt5,3 )F RBEGVFREND 

3 FORMAT ( 2F10 # 3) 

D04l = l';CE 


OIF=A35(FR0EG-FR(I) ) 

if(jif;le,o.o7)go to - 
CONTINUE 
JBEG-I 
D06l=l»LE 


‘Dip=A3$(FRtND-TTrrrn — 
IF(DIF.LE.0.07)GO TO 7 
CONTINUE 

jend=i 

IF{ JBEG.GEVJEND) JBEG=r 
IS frequency in ghz. 


6 

7 

C Fr 

T w is INSUIDE wAVENUMBEITTNrCM; 

C WL IS INGUIDE WAVELENGTH IN cm, 

C FsWL IS FREE SPACE WAVELENGTH^ IN TW.” 
0021=1, LE 

WN(I)=SQRT((FR(I)/29VC 
WL(I>=1./WN< J) 


FSwLT I ) =29 , 9695VFR (!) 

BETA(I)=2.*PI/WL(I) 

2 CONTINUE 5 

c initial values for the physical parameters of the matching element, 

REaO<5,4001)NOO ~ 

4001 FORMAT ( 12) 

* D030 0$NR = I , N§0 


111* READ(5,4000)ALENGH(1) 

112* ' “4'bOo FORMAT (FlO, 5) 1 

113* WRITE(6»5001)ALENGH(1) 

■ 114* ^001 FORMAT(///,10X, ♦STARTInS”VALUE OF~AL^NGH~ ^ F TO , 4) 

115* MM=0 

~IT 6 * 0 P ( 1T£732" 

117* IF(mACHEL,GE.3)0P(1)=,i6 

lie* “width ( 1 ) = • o5 

119* IF < MACHEL • GE, 3) WIDTH( 1)=, 47 

“120* c increment physic al“*parametERS^rounO“ThE“ ( ivm ) point xn — 

121* c dimension space, 

122* c this program preforms^anhi norememting “addi istment to zoom in; ” 

i ANroMT-i 

~124* IF(M«CHEL.LT.3)AMUM0R=;i 

125* IF(V|ACHEL.GE.3)AMUMBR=.05 

'126* ' 30 TO 250 "" 

127* 35o ANCR^TrANCRMT*2, 

~T26* 25o~OP“l 2 V=3P ( 1 j“ANUMBR7flTlCFfMT 

129* FA L 3D=0P(1) 

’130* FALBW=wiOTH(ll 

131* XF(0P(2) .LT..04,AND.MACHEL.LT.3>0P(2>=FALBD 
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132* IF(DP<2) ,LT. .005.AND,MaCHEL.GE,3>OP(2>=FALBD 

"133* 0P(3)=DP(1)+AMUMBR/ANCrMT 

134* WlDTH(2)=WIOTHtl)-,01/ANCRMT 

135* IF ( wI 0TH(2) ,LT,.03,AND.>1ACHEl.iLTr3)WlOTH(21-=FAtBW 

136* IF(wlDTH(2) ,LT. . 46. AND. MACHEL.GE. 3 ) wIDTH(2)=FAlBW 

137* WIaTH(3)=WlDTH(l)+.01/ANCRMT 

138* ALEnGh ( 2 ) =ALENGH 1 1 > -, 1/ANCRMT 

139.’ ALENGH ( 3) =ALENGH ( 1 ) + , 1/ANCRMT ~ 

140* C RETURN TO RECALCULATION, 

I4i* c Fixed mismatch is gam; — 

142 * c the fixed mismatch at a point alengh toward the generator is agam, 

TH* D0302M=173 

144* D0303n=J3EG» JEND 

145* AGAMtM,N)=GAM(N)*CEXP(fO>-lT*2';'*BETA(NWUENGH(Mn- 

146* IF { ALENGH ( M) ,LT,0.)AGAM(M*N)SGAM(N) 

147* ' '303 CONTINUE — ' ~~ 

148* 302 CONTINUE 

'149* “ GO TO ( 399 1 49 < 5i“5997699r>MACHEL 

150* C**************«******************************************************** 

"151* " " 399 CONTINUE " 

152* c matching element"Capacitive iris* 'Complex reflection coeff. 

153* C" As A FUNCTION OF DIMENSION SPACE AND FREQUENCYr'7 

154* D0300L=J3EG> JEND 

155* A9 = 50RT(1.-(3/WLTCT7**2T 

156* AA(l)=(1./A8)-1. 

157* "300 CONTINUE " "" " "* 

158* C V. I OTH ( I ) is EQUAL TO THE THICKNESS OF THE IRIS, 

159* ' C Dp{ I ) IS EQUAL TH TWICE THE IRIS HEIGHT* 

160* 00252J=1*3 

161* D(j)=3-0PCj> 

162* 002511 = 1.3 

163* IF(I .EQ.l.AMD.J.EQ,l,0R.I.EQ,l.'AND, J,EQ.2.0(r.T.EQ.17ANOTJ.EGrr3,0Rr 

164* 1I.E0*2.AND.J.EQ.1.0R.I,E9.3.AND.J.EQ.1)G0 TO 2003 

165* """ GO TO 251 

166* 2003 CONTINUE 

~167* COnD = wIDTH ( I J /DP ( JJ 

168* IFfCOND.GE.l, )DP( J)=FA l3D 

169* IF ( CO.NO ,Gt, 1,1 WIDTH (I) =FALBW 

170* GM;1.+(WIDTH(I)/(PI*DP(J) ) )*AL0G(4,*Pl*DP(d)/tE*WI0TH<I))> 

'171* C0( J)=C0S(PI*0P(J)/(2,*3) ) — 

172* SI ( j)=SIN(PI*DP(J)/(2,*3) ) 

T73* SE(‘J)=i./C0S(PI*DPr3J*sH/T5T*ffrJ 

174* _ 0030 1<K=J3EG. JEND 

175* 90nE=PI*0P(J)*WIDTH(I)/(2.*R*D(J)) 7 

176* 8TwO=AA(KK)*(SI ( J)**4)/tl.+AA(KK)*(C0(J)**4> ) 

177*' 3TRE=d./16.)*( (3/wL(KK))**2)*t'(TT-3.*(C0(jr**2TT»*2T*TST'njT^'«r) 

178* 3N0RM=(2,*9/wL(KK))*(Al0G(SE< J) ) -BOnE+RTWO+BTRE) 

“179* 3A = 3Norm+(B/0(J) r*TAN(pI*WIDTH(T)/WL<KKn 

180* 33=(S/D( J) ),1./SIN(2,*PI*WI0TH(I)/WL(KK) ) 

“181*' YI = (0,D*3a + 1, 

182* Zll=< (0.1)/B3)+1./TI 

183* YTr (0, i ) *BA + 1 ,/ZII 

184* 00260‘1=1»3 

“las* 1 b g a m = ( rTJiTT? rrr+TTi 

186* GGAMd, J,KK.N)=BGAM*CEXPt < 0 » 1 ) *2 . *BETA (KK) *ALEnGH(N) ) 

"187* I F ( ALENGH < N) , GE. 0 ; ) GGAM (I *U7KK", N)=BGAM 

188* 280 CONTINUE 
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189* C I CORRES, to THICKNESS OF THE IRls. 

190* C J CORrES' TO HEIGHT* OF~TRI ST 

191* __ C <K COrrEs TO FREQUENCY, 

192* " C SsAM is THE COMPUErREFLEcTrON-C^FFr^RERT^HITTANCr^^ 

193* c smith chart plot, 

~T94* 30~i CONTINUE 

19S* 25i CONTINUE 

196* ”' 252 CONTINUE ' ' ' “ 

197* GO TO 404 

198*'“ C* ******************** * * ***** **'***'***+ * * *>* *****r*'* ****** * * **** *********** 

199* 499 CONTINUE 

2oo* c” matching element _ TASYMMmrcAr^rroucTTvE~TRrs T “coMPLE)r*EFi-ECTTON 

201* C As A FUNCTION OF DIMENSION SPACE ANO FREQUENCY, 

202* ‘ * “ D0400 1 rl * 3 ' " : " 

203* D040lJ=l»3 

204* " " IP{I.EQ,1.AND,J.EQ;1,0R.T;EQ,I;AND,J.EQ,2,0R,I,EQ.1.AND,J,EQ.3,0R, 

205* 1I.EQ. 2. AND. J. EQ.l. OR, I. EG. 3. AND. J.EO.l) GO TO 2007 

206* ' ' 30 TO 4ol “ 

207* 2007 CONTINUE 

208* COnD=vJIDTH(I)/DP(J) * 

209* IF(C0nD,GE.1.)DP(J)=FAlBD 

2ili* ’ IF(COND.GE.l t )WIOTH(I) = FALBW “ 

211* CApDP=(DP< J)/SQTWO)*tl.+(WlDTH(I>/(Pl*DP(J) ) ) *AL0G(4.*PI*0P< J)/(E* 

212* 1 width cm r> 

213* _ CAp;)=( (4,/(3,*Pl) )*WlDTH(I)*DP(J)**3)**(l,/4, > 

214* D0402KK=J3EG, JEND “ 

215* C XA AND X3 ARE NORMALIZED IMPEDANCES. 

216*" " XA=(4.*A/WL<KK))*((A/(pI*CAPDP>>**2) 

217* X8=( A/{ 16,*WL(KK) ) )*( (PI*CAPD/A)**2) 

“21b* “zl"=(0,-i)*xa*17 

219* _ ZXls<l a /( (l,/( (0>1)*XA) )+l./ZI) ) 

220* ~ ' ZT=( 0»-l ) *X9+2II 

221* 3GAM=(1,-2t>/( l.+ZT) 

222* " D0403n=1>3 “ 

223* GGAM(I f J,KK»N)=BGAM*CE xP( (0»1)*2.*BETA(KK)*ALEN6H(n) ) 

224* " IF { aLEnGH <n) .GE’.OVTGGam CITJf KK7N)"=H5^W 

225* 403 CONTINUE 

226* " “<T I CORRES, TO THICKNESS”0F THE “IRIS, “ 

227* C J CORRES. TO HEIGHT OF IRIS, 

228* ' C” Kk CORRES, TO FREQUENCY.' 

229* C GGAM IS THE COMPLEX REFLECTION COEFF, AS PER ADMITTANCE 

“230* C - 5mITH“'CHARt“PlOTT 

231* 402 CONTINUE 

' 232* 40 1 CONTINUE 

233* 400 CONTINUE 

"234* GO TO 404 — 

235* C********** ************************************************************* 

236* 599 CONTINUE * 

237* __ c matching element 'solid inductive post* complex reflection coeff. 

236*” c As A FUNCTION OF" DIMENSION P 

239* C As A FUNCTION OF DIMENSION SPACE AND FREQUENCY. 

240* C DP ( I ) IS DEFINED TO BE THE DIAMETER OF THE POST. “ 

241* C WIDTH ( I ) IS DEFINED TO 9E THE DISTANCE OF THE CENTER LINE OF THE POST 

242* c from the side-wall; ~ 

243* 005401=1,3 

" 244* IF(wlDTH(I) ,GT.1.82)WlDTHfI)=FALBW 

245* D054lKK=J3EG, JEND 
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24fc* M=2 

'247* M=2 

246* XSH = 0,0 

"249+ X500-0 . 0 “ 7 

250* 543 XlrS^RTi l./( (N**2)-( (2 # *A/FSWL(KK) )**2) ) ) 

251+ X50=((SlN{N*Pl*Wl0TH(ir/An**Pr*t){Iit;7NT 

252* XS00=XS0+XS00 

253+ ~ IF(xSo.UE,.OOl) SOTO" 5<f5 ~ 

254* 544 N = N + 1 

”'255+ GO TO 543 . 

256* 545 XIrSORT(M**2/( <M**2)-< ( 2, +A/FSWL < KK ) ) +*2 ) ) ) 

2 3 7+ XSl =TTsIR( 2 , *M*PI +WIDTH < I J /A ) ) * (XI-1 • M 

256* XSH=X5I+XSI I 

259* A3xsI=ABS<XSI ) 

260* IF ( ABxSl .LE. .OODGO TO 547 

261* “ 546 M=V|+1 ~ " — — 

262* GO TO 545 

"263+ 547' CONTINUE 

264+ SOR( I ,KK)=XS00*2 t -2.*{ (SIN(PI*WlOTH ( I ) /A) ) **2> 

265* SII (IrKK)=xSlI + C t 5*(l./TANCPX*wlOTH(n/Aj }-SlM(2V*Pl*WI0TKTimTT 

266* 54 1 CONTINUE 

267* 54o CONTINUE . 

266* 005461=1,3 

269+ CT=1 • /TAN(pI * WIDTH! TT7 Al 

270* CSS=< l ,/5lN(Pl*WI0TH(I )/A) )♦*? 

271* 5NS=(5IN<P1*WIDTH(I)/A) )**2 

272* 00549J=1,3 

273* ' IF(I.E;0.1.AND.J,F:Q,I,OR,I.EQ.l.ANOiJ,EO,2 # OR;T.rQrrrANa^7rOv3,t)Fri 

274* H.E0.2. AND, J.EG.l.OR. I . EQ . 3 . AND . J . Eq . 1 ) GO TO 2008 

275* '"GO TO '54 9 

276* 2008 CONTINUE 

277* IFoP(J) , 3£ , # 22)DP(J)=pALBD 

278* IF (DP( J) ,L£, ,0)GO TO 200 

279* PD0A-(PI*Up(J)/A)**2 " ” ' 

_2Q0* PDo= < PI *DP ( J) /(A*2 t > )+*2 

261* 50pp=AL 0G( (4.*A/(PI*DP( jn*SrN(Fr*WrOTRTrT7AT)T 

2Q2* D0550<K=J3EG, JENO 

263+ ~'C" XA AND X3 ARE NORMALISED" IMPEDANCES* 

264+ XB=(A/ W L(KK) )*POOA*5NS 

265* ‘ XA = xB/ 2.+A/(2.*WL<KK) } *CSS* < SOPP+SOP ( I , KK) - < (Pl+DPT J) /T2. +70 * **2) + 

266* 1 ( ( ( SOPP + SQP ( I # KK ) ) +CT-SII ( I *KK ) )**2) ) 

287* ■ 21= ( 0 ,- 1 ) *xB+l t 

266+ ZIl = (l,/( d # /( (OtU+XA) >+l./ZI> ) 

269* 2T= ( 0 , -1 ) ♦ X3+ZI I 

290+ BGAM=(1,-Zt)/<ZT+,1) 

29l* 0055l< = l , 3 ” 

292* GGAMd, J,KK#K)=3GAM*CEXP( t 0 • 1 ) *2. +9ETA ( KK ) *ALENGH< K) ) 

293* rF(ALENGH(Kr,GE,OVO)T3GAMTT7DfK<VK)=B5AM 

294* 55i CONTINUE 

'295+ “ ” C I CORRES TO DISTANCE OF PoST FROM SIDE" WALL”, 

296* C J CORRES TO DIAMETER OF PoST * 

"297* C Kk corres to frequence; 

296* C K CORRES TO DISTANCE BETWEEN MISMATCHES* 

* 299*“ 556‘ CONTINUE 

300* 54g CONTINUE 

“301* 54q CONTINUE " 

302+ GO TO 404 



205 


303* c*** ********** ********************************************************** 

“304* 699 CONTINUE 

305 * c matching element ’dielectric post* complex reflection coeff. 

"306* " ~C A s A fu n cti on oftiimensi on~space~and~frequency» - given the dielectrics 

307* SQEPRm=SQRT<EPRIM) 

309* C WIDTH is DEFINED TO 9E THE DISTANCE FRoM ThPsIDEwALL^TO THE CENTER- 

310* c "Line of the post; — 

311* C 3ESLA IS THE 9ESSEL FUNCTION WITH ALPHA AS ITS ARGUMENT. 

312* C BESLB IS THE BESSEL FUNCTION WITH BETE AS ITS ARGUMENT. ~ 

313* 00602j=l»3 

“314* IF<dP(J) . G E 773 4 2TDP PJ rEF ALB D 

315* D0fcn3N=JBEG. JEND 

316* ■”* ALPHA( JrN)=Pl^DP(J)/FSwL(N) ” “ 

317* __ BEtE<J,N)=ALPHA(J>N)*SoEPRM 

318* 3ESOA( J > N)=3SSL(ALPHA(j#N)"#IT 

319* BES0B< J»N)=9SSL(RETE( J,N) ♦ 1) 

320*" BESIM Jf N)=9SSL( ALPHA CjVN)T3'T 

321* _ _ BF.SI3C J>N) =B SSL (BETE ( J,N) >3) 

322* t' “J CORRES TO DIAMETERS N CoRRESTD'TREZ^SPACF WAVELENGTH; “** 

323* _ _ 603 CONTINUE 

324* '* 602 CONTINUE “ 

325 * c Calculate so the partial sum. 

326* C I CORR£S'TO DISTANCE'TROM SIDEWALL. 

327* C L CORRES TO rREE-SPACE WAVELENGTH. 

'326* 'D06U 4 I = lr3 ■" - 

329* IF( w lDTH (U .LT. .46) WIDTH (I )=WI0TH(1) 

330*“' DObOSL-JBEGr JEND *~ “ 

331* N=2 

“332* XS60=070 

333* 606 XI=SQRT<1./<(N**2)-((2.*A/FSWL(L>)**2>)> 

334* XS 0 =MSIN(N*PI*WlDTH(I)/An** 2 )*CXI-l./NJ 

335* XSOO=XSO+XSOO 

336* IF(xSO,LE..OOl)GO~TO“6o8 

337* 607 Ns n 4-1 

336* 50 TO “6(1 6' * 

339* _ _ 600 CONTINUE 

340* “ SOp ( I » L ) = XS00*2'.-Z7*”( ( sIN(PI*WIDTH( J) /A)T**2) 

341* ____ 605 CONTINUE 

'342* 6 O 4 CONTINUE ~ ~ “ 

343 * c calculate the shunt and series impedance for the post model. 

“344* "006091 = 173 

345* 0061 0 Jr 1 » 3 

346* IF ( I . EO. 17 AND. JYEQ71 .OrTITEQTT. AND7UVE0V2.“0R .l7E07T7^ND“.tJVETT73,T5R'; 

347* _ __ 1I.EQ.2.AND # J,EQ.1.0R,I.E0.3.AND.O,E0.1)GO TO 2009 

348* “ “* 30 T° 610 

349* 200g CONTINUE 

35o* c r correS“to 'Di stance~from the - sidewall; ' 

35i* _ _C J CORRES TO DIAMETER OF DIELECTRIC. 

352*' ' " SOpp = ALOG( (4,*A/(PI*DP{ J) )*SIN(PI*WIDTH(I)/A) ) ) 

353* 0061 lMr J3EG» JEND 

'354* C M’ CORRES TO FREQ. 

355* XB=(2,*a/wl(M) )*( (PI*D p( J)/A)**2)*(SlN(PI*WlDTH( I)/A)**2)/( ( (ALPHA 

"356* TCD7M) ♦♦2)*tiESIR(J» WJ7BESIA( J^M) ) 4^(1 ;/( ALPHA (J»M)*BESOA(J»M)* 

357* 2BESI3( J#M)-BETE( J»M)*BESOB( J»M)*9ESlA( J,M) ) )-2, ) 

358** * ’ XA = ,5*X34< A/(2.*WL(M) ) ) * ( 1 . / ( SIN ( PI *WinTH ( I) /A ) **2) )* ( (BES08 1 J*M) / 

359*_ __ _ IBEsqA { J » M ) )*( l./(BETE(jfM )*BESOA(J,M) *RESlB(JfM) -ALPHA (JfM)* 
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360* 

361* 

362* 

" 363* 
_ 364 *_ 
365* 
366* 

" 367* 
366* 

~ 369*~ 
370* 
~zin~ 
372* 
373* 
374* 
375* 
376* 
~377*“ 
376* 
379* 
380* 
381* 
382* 
383* 
384* 
385* 
396* 
'387* 
388* 

~389*“ 

390* 

391* 

392* 

393* 

394* 

~ 395*” 
396* 
397* 
398* 
399* 
400* 


2B£509( J#M)*RESI A( JrM) ) )-<SOPP+SOP<I #M) > ALPHA < ) ) 

Yl = l./<(0a>*X3+i;) 

Yll=l./i (0,-l)*XA)+YI 

YT=l,/( ( 0 » 1 ) *XB4l */YII ) 

BGam=<YT-1.)/(1*+YT) 

30bl2Mil»3 

GGAM< I , J»M,N)=BGAM*CEXP( < 0 # 1 ) *2 . *BEtA ( M) *ALENGH ( N ) > 

IF( ALEMGH(N) ,GE.0* )GGAM(1» J'»MtN)=8GAM- 

612 CONTINUE 

X'“ GGAM IS THE-'COMPLEX'RErCECTlOfrCOEFP'A^ PER~AWlTTANCr SMXTH-XHARt^ 

61 1 CONTINUE 

SIFtONTlNUE 

609 CONTINUE 

“ GO TO 4 04' 

C****** ****** *********************************************************** 

"" 404 CONTINUE 

c Calculate the vismatch and the rms of the mismatch, then determine 
t the smallest value anetexitor - return TO~TNCREWENT"Ttf oimv space, 

DO20I I = l» 3 

' ■ D02G2j = l»3 

D02n3<=l,3 

IF( I *EQ, 1 . AND.U.EG, 1 ,AnD,K*EG, 1. OR,I,EQ,IiANO, J^EQ, 1,AND,K.EQ* 2VOR 
1,1 .£0,1 . ANQ. J.EQ.l.AN0.<.£Q,3.OR*I.£G.l.AND,J.EQ.2.ANO.K,EQ,l,OR,I 
2 ; Eq . 1 ; A ND . J . EG • 3 .“AND , <7 EG .1 .DR ~T V E 3 . 27h ND .LT.EQ .17*NO.KT EG^TTOr , T,E“ 
3Q.3.AN0. J.EO.l. AND,K,E q,1)GO TO 2001 

GO T° 203 1 

200) CONTINUE 

PARTSv.rO.O 
0020 4 N=J8 £g» JEND 

G A vi M A H - ( A G A M ( X » N )+GG AIM ( I # J , N , K ' ) -2 .* AG AM ( K7NT* G G AM TT»s 
1AM(<»N)*G3AM(I, J»NrK) ) 

SNGLVL=(CA3S{GAMMAH) )**2 
PARTSvrPARTSM+SNGLVL 
204 CONTINUE 

RS:S°RT<PArTSM) 


arms < I • J» K ) S ( r; /SORT {UEND-JREG JT*Rr 
C The ROOT yrA.N SO. OF THE MISMATCH OVER THE BAND OF INTEREST HAS BEEN 

~c Calculated, — — 

203 continue 

20 2 CONTINUE 

20l CONTINUE 


4oi* c Determine the" 

402* c center point in dimension space, 

403* ‘ C Arms is the root mean SQUARE MISMATCH OVER theband of interest, 

404* C DETERMINATION OF the LEAST OF THE ARMS ENSEMBLE* 


405* 

406* 

ATEST=ARMS(1,1 

002 0 5 1 r 1 » 3 

407* 
408* 
”409* 
410* 
“ 411* 
412* 

D0206J=1»3 
D 020 7 K =1 / 3 
IF ( I *EO« 1 * AND*) 
1*I,EG*1*AND. J.l 
2 ,Eg, 1. AND. J.EQ 
39. 3. And, J.EO. l 

413* 

GO TO 207 

414* 

2004 CONTINUE 

415* 

IF( ARMS(I VU»K)‘ 

416* 

207 CONTINUE 
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417* 206 CONTINUE 

41b* ” ' ' 205 CONTINUE 

419* D0l00l=l»3 

'420* OOl0U = lr3 

_42 1 * DOio2K = l#3 

42 2 * IF( X • E~Q , 1 *aNd*U»E 0 * 1 «AmD,K*EQ* i«ORtI*E ff» I. • "A NO • J •EQ«1*AND*X# E® * 2 • OR 

423* l.I.EO, l.AND* J.EQ.l, AND . K ,EQ . 3* OR • I .E&. 1 . AND. J.EQ. 2. AND. K . EQ, 1 . OR. I 

424*' 2.Eo, 1 , ANO/J,E3 # 3VAND7K.EQTI70R;i .EO;2VAND7J.EQa ,AND.K,£Q.1.0R/I.E 

425* _ _ _ 3Q.3, A.NP, J« EQ. 1« AND , K. EQ* 1)60 TO 2005 

426* ' ' 30 TO 102 ” ' ' ” 

427* 2005 CONTINUE 

"428* IF ( A3S ( A TEST- ARMS Ft VJ # K ) - ntT7T;E*Wr®0“TO“ItJ3 

429* 10? CONTINUE 

430* ' 1 0 1 CONTINUE 

431* 10o CONTINUE 

432* 103 CONTINUE 

433* MM=MV+i 

'"434* C TrST'"FO R TERMINATIONT 

435* IF(\lH.EQ.l)GO TO 108 

4 36* IF(I.NE,1.0R. U.NE;1*0R;K7NE71TS0 TO" 108 

437* IF(ANCRMT,GE.Q. )30 TO 200 

436* GO TO 350 

439* 10a continue 

440* WIOTH(l)=wiOTHTI) “ 

441* _ _ OP(i>=DP(J) 

442*' ALeN^H( 1 )-ALENGH 0C) 

443* ITi=lTIME(IT2*IT3) 

444* ““ ’ IF ( I Tl # LT • l , 5E*5) GO” TO 200 

445* GO TO 250 

446* Z ReSUL f ANT" MI SMATCH“PROCEE dURE 15 C»*RIED OUT, 

447* 200 CONTINUE 

44B* WRITE ( 6# 370) I > J*K 

449* 37o FOrmAtuH!,//# 10 X# ’OPTIMUM CHOICE HAS INDICES ( • » U# 9 t * » lit • • * r II* 

450* !♦)♦) “ — " T " 

451* WRI T= ( 6 # BOO ) 

‘452* 8 O 0 FORMA t ( // # 1 OX #’* COORDINATE rrVJTX7 f »T0X7 r ROOT” MEAN SQUARE OVER MATCH- 

453* 1E0 a fl ND » ,//) 

454* DOfiO^I 1=1 *3 ' ~ — 

455* _ 00603JJ=1#3 

456*'““ D060 4 KK = 1 * 3 

457* IF(lI # EQ,1, ANO. JJ,E0,l.AND # KK*EQ*1.0RtII # E0*l«AND,JJ*E0,l.AND,KK t E 

“456* 10 7 27 0 R 7 1 1 V E Q I7A N D7JJ7F07T7A Ntr,X X 7E0T3 . 0 R7TI7EQ , 1 . A Nt>7J J . E Q • A N D ; 

459* 2KK,EQ.l,0fMI ,EQ. l.AND, J J. EG. 3, AND f . EQ . 1 . OR . 1 1 . EG . 2 , AND . JJ . E<5 . 1 . 

460* 3AND.KK.E0, i,OR,II .EQV3. AND."UUVEQ*1.AN0 # KK,EQ # 1)G0 TO 20 06 

461* GO T° 604 

462*"““ 2006 CONTINUE * “ 

463* WRITE (6* 801 ) 1 1 # JJ > KK # ARMS ( 1 1 # Jj# KK ) 

464* BUi FORMAT i 21X7* ( ~' f 7 1 lT*Ti f J IT*# * VI I * *" ) **7 1 0X7 FI 0 • 7 ) 

465* _ 8 O 4 CONTINUE 

466*' ' " 8 O 3 CONTINUE """' ‘ "" “ 

467* 602 CONTINUE 

468* GO TO (396,498# 598,698) RACHEL 


4 69* c*** ****♦*♦♦♦♦*****♦****<{**♦*****♦♦ ************************************* 

47 C* 39$ CONTINUE 

471* _ D(J)=6-DP<J) 

472* D0395KK = 1»LE 

473* AB=50RT( 1.-(9/WL(KK> )**2) 
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474 * AA<KK)r<i./AB)-l. 

“ 475* GM = 1. + (WIDTH( X)/(PI*DP{ J) ) )+Al.oG(4V*PI*t)Pl“J)/(t*WlOTKnm 

476* C0(J)=C0S(PI*0P<J)/<2.*B>> 

477* 5I<J>=SlN(pl*0P<J}/f2.“*BJJ 

478* SE( J)=1./C05(PI*0P< J)*G^/<?.*B) ) 

“479* 10mE = PI*3P( Jj*WlDTH{rr/(27*A*0TJ)T 

430* ‘ BTwO=AACKK)*(SI(J)**4)/(1.+AA(kK)*<CO{J)**4)) 

491* 3tre=(i./ 16.>*< (3/wL(KK))**2T*(a.-3.*(Cocj)*#2m*2msrnn»»Tn 

402* BN 0 R m = ( 2 , *3/WL ( KK ) ) * ( AlOG ( SE ( J) > -BOnE+BTWO+BTRE) 

" 483* 8Ar3N0R'*t+ (B/D ( J)) *TAN( pi* WIDTH ( I ) /WL (KK7T 

404* _ BB=<a/ 0 ( J) )*l./5IN(2«*pZ*WI0TH(X)/WLCKK)) 

485* YI= ( 0, i ) *3 a+1 , 

486* 2 ZlsC fo,l)/B3Ul,/Yl 

487* YT=(0,l)*3A+l./ZII 

480* 3Gam=(YT-1,)/C1 # *YT) 

489* SGAMt l, J^KKrK)=3GAM#CtxPM0#l)*2.*9rTA(KK)*ALEMGH(K)7 

490* IF{ALrNGH(K) , SE , 0 . ) GG Am ( I # J # KK f K > =B gAM 

491* AGAMtK,KK)=GAM{KK>*CExp((0r-ir*2i*3ETM'XKT*^UEW5mKT7 

492* IF{ALE.^GH(K) .LT^O.) AGAm(KpKK)=GAM(KK> 

493* 3A^mA(kK) = ( AGA^(K*KK) 4-GGAM(I# J,KK#K)-2.*A5AMCK>KK:)*GGAM(T7-JrKK7X7 ) 

494* l/(l .*AGAM(K#KK) *GGAM(i , J,KK»K) ) 

495* “ 395 CONTINUE ' “ 

496* DHr }P ( J ) / 2 , 

■"497*' WRlT^(6»226l ' 

490* 226 FO^MAjdHl,///, lOXf *CApAClTIVE IRIS IS THE MATCHING ELEMENT*) 

499* ’*”*""* WR1TE{6»242) ANCRMT ~ : * 

500* 242 FORMAj(// 1 jOXf • NUMBER OF INCREMENT = *,F4,0> 

501* WRITE(6#27o)FR( JBEG) VFR(JENr)) - 

502* 27o FORMAK//H0X# ’MATCHED BAND IS FROM *,F 7 t 3 ,» GHZ* TO *>F 7 . 3 , f GHZ* 

“503* I 1 ) “ 

504* WRITE ( 6 » 227) WIDTH (I ) , Dm # ALENGH ( K ) 

505* 227 FORVAT(//, 1 0Xr*THICKNEsS= • »F6#4f lx* ♦ IN”CM» # 10X7*HEIGHT= * *F6T4 

50b* _ lr* IN CM» » lOXr ‘DISTANCE TOWARD THE G^N. FROM MISMAtCHz *#F8*4) 

507* ““ ' GO TO 999 

506* c************* ******* *************************************************** 

'“509*“ “499 CONTINUE ' 

510* CApDP=(DP( j)/SOTWO)*( 1 ,+(WIDTH(I)/(PI*DP(J) > ) *AL0G(4,*PI*DP( J) /<E* 

511* " ' IWIdtH(i)))) “ 

512* CApD=( ( 4 * / ( 3# *PI ) )*WIDTH(I)*DP(J)**3)**d./4.) 

513* DO 496 KKXI rtE 

514* C Xa and are NORMALIZED IMPEDANCES# 

515* XA=(4 # *A/wLUK))*C(A/(pI*CAPDPn**2T 

516* _ X3=(A/(16.*WL<KK) ))*( (p I *C APD/A > **2) 

517* ZI=(0#-l)*x6*l, “““ 

518* ZIl = ( 1 ,/ ( < 1 ,/C (0» 1)*XA) )*1,/ZI) ) 

519* ZT=(o,-n*xa*zn ~ 

520* BGAM=d,-ZT)/(l.+ZT) 

“521 * GGAMCIV J,KKr<)=9GAM*CExPrtOTrT*2T*BETATKK)*ALENGH(Kn 

522* IF ( ALENGH (K) ,GE,0. JGGAmCI »J rKK,K)=BGAM 

523* ASam(K,KK)=5AM(KK)*CEXP< (0# -I U2. *BeTA^XKT*ALENGH(kT1 

524* IF( ALtnGH (K) ,LT.O.)AGAm(K*KK)=GAM(KK) 

525* GAMMA(KK) = (AGAM(K#KK)+GGAM(I,J f KK,K)-Z7*"AGAMTIC7KX}iG^AMtlTavK^Kn~— 

526* 1/<1,*AGAM(K»KK)*GGAM(I, J#KKrK) ) 

"527* 496“C0'ntInuE 

528* DH=QP(J)/2. 

"529* WRITE ( 6» 405) 

530* 405 format<ihi» ///» iox# * inductive asymmetrical iris is the matching el 
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531* 1 EMEmT») 

“ 532* WRITE ( 6 » 406) ANCRMT 

533* 406 FORMAl(//r 10X» ‘NUMBER OF INCREMENT = *,F4.0> 

534* WR 1 7- ( 6 # 470 ) FR (^JBEG)TFftTJENTT) 

535* 47n FORMAT(//»iOX» ‘MATCHED BAND IS FROM f *F7.3*« GHZ, TO * *F7.3» * GHZ. 

536* I T T 

_ 537* WRITE<6»407)WZDTH(I>#DH*ALENGH<K> 

538* 407 FOR \lAT(//» l OX# * THICKNESSES - * VF 6 , 4 riyi ‘CM.iV10XV»HElGHT“5 -T7r674#' f 

539* _ 1CM.“10X» ‘DISTANCE TOWARD THE GEN. FROM MISMATCH = ‘#F8.4) 

540* ’ " GO TO 999 

541 * c* ****** ************* *************** ********cOMMON*OUTPUT***** ********** 

542* 59e continue 

543 * _ c matching element ‘solid inductive posTt complex reflection coeff, 

544* c As A FUNCTION OF DIMENSION P 

546* _ C A S A FUNCTION OF DIMENSION SPACE AND FREQUENCY. 

546*" “ C OpU) IS DEFINED TO 3E THE DIAMETER OF THE POST. 

_ 547* C WIDTH(I) IS DEFINED TO BE THE DISTANCE OF THE CENTER LINE OF THE POST 

546* £“'TrOM'th5“‘SIDE-WALL'. “ 

549* _ __ D0525KK = 1#LE 

550* * Ns 2 

55l* M=2 

"" 552* ' XSI I = 0 , 0 " ~ “ 

_ 553* XS00=0,0 

‘554 * 513 XI=S q RT f lV/( ( N**’2T-TT 2 “*A7FS¥C ( KK n **2 IT) 

555* XSO=t (SIN(N*PI*WIDTH(I)/A) ) **2 ) * ( XI~1 ,/N) 

556* XS00 = XS0fXS00 

557* IF(X50.LE. ,OOl)GO TO 5l5 

558* ' ‘ 514 N = N*1 

559* ' GO TO 513 

560 * 57 5 X I : S OR T ( M * * 2 7TT M**‘2’T~T (‘Z7^A7F3WL t KTyT**2TTj 

561* XSI=( (SIN(2.*M*PI*WI0TH(I)/A) )*(Xl-i. ) ) 

“ 562* “ - - XS1I=XSI*XSII 

563* A9xSl=ABS(XSI ) 

564*' " - ‘ IF(A3XSI.LE, ,001)30 TO'51’7 

__565* 516 M=vul 

56b* GO "TO “Si 5 

567*_ 517 CONTINUE 

568* ' SOp( I #KK ) =xS00*2.-2T*tTSIN(PI*Wl0THTlT7ATT**2l 

569* SI I (I #KK)=X5II + ( .5*(1./TAN(PI*WIDTH(I)/A))-SIN(2,*PI*WIDTH(I)/A) ) 

57C * *“ "525 CONTINUE " " 

571* CT= 1 ./TAN(pI*WlDTH(I ) /A) 

572* CSs^IiT/SlN^Pl^WlDTHriy/rJTT*^ 

573* SNS=(SIN(PI*WIDTHU)/A) )**2 

574* " PD0A = (PI*0p(J)/A)**2“ " 

575* PD0=(?I*DP(J)/(A*2.) )**2 

57b* - ' SOPP = AloG< ( 4,*A/(PI*DP( J) ) *SIN( Pl*WlOTH < IT/A) ) ) 

577* D0526KK = 1 »LE 

57b* Z XA And *3 ARE NORMALIZED' IMPEDANCES* 

_ 579* X3=(A/ w l<KK) )*PDOA*SNS 

580* XA = x 0 / 2 .+A/( 2 ,*WL(KK) ) *CSS* ( SOPP+SOP ( I , KK7 - { (Pl*DP ( J)7 ( 2, *A) )**2)* 

_ 581* 1 ( ( ( SGpp + SOp ( I > KK ) )*CT-sII C I f KK) ) **2) ) 

58?* ZX = (Oi.l)«xS*l. 

583* Zll=<l,/< (l./( (0»1)*XA))+1./ZI) ) 

“504* ZT=(0r-*l) *X3+ZII — 

565* 3Gam=( i ,-ZT>/<ZT**l> 

586* ‘ * GGAM(I t J,<K#K)=BGAM*CExPt to » 1 ) * 2 . *BpTA (KK) *ALENGHt« ) " 

567* IF ( ALENGH<K) .GE. 0 . 0 ) GGaM( I » J r KK »K)=0GAM 
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58B* C I C0RrE$ TO DISTANCE OF POST FROM SIDE WALL, 

589* C 'J CORRES TO DiAmETER OF PoSH 

590* C K* CORRtS TO FREQUENCY, 

591* C K CORRES TO DISTANCE' BETweEN~MI5«ATCHES7 

592* AGAM(KrKK):GAM(«)*CExp( (0»-l)*2**8ETA{<K)*ALENGH(K) ) 

593* IF(ALENGH(K),LT,0.0)ASAM{KVKK)i'SAM(K*n 

594 ♦ GAMMAUf<) = {AGAMUrKK) + GGAM(I# J,KK#K)-2 # *AGAM(K,KK)*GGAM(Ir J#KK,K) ) 

"'595* i/(l.-AGAM(K#KK)*GGAM(I"J,KK»K)J 

595* 526 CONTINUE 

~ 597* WR I TE ( 6 » 519) 

596* 519 F0RMATtlHl,///#10X#'SOLlD INDUCTIVE POST IS THE MATCHING ELEMENT*) 

“599* WR1 TET67 5 2 0 VANCRMT 

600* 52q FORmAT(//»1OX»*NUM0ER qF INCREMENT = *,F4,0) 

“ 601* WRITE (6lr 521 ) FR ( J3EG) #FR( JEND) 

602* 521 FORMAT (//# 10X* * MATCHED BAND IS FROM *»F6,3*« GHZtS* TO *#F6*3 G 

603* 1H2. ») 

604* WRITE(6»522)0P( J) ^WlOTHd > fALENGH(K) 

‘605* 5Z2"F0rviAT <7/ * lOX V r DI AMETE r"=~ t »F 679T DISTANCE "sr 

606* _ l#P6. t +»’ CM, ♦ » 10X» ’DISTANCE TOWARD ThE GEN, FROM MISMATCH = trFe.4) 

607* GO TO 999 ' ' 

606* C*********** ************************************************************ 

609 * ■' ‘ 69e continue 7 

610* DOb20N=l#LE 

611* ALp^A(j,N)rPi7DP[Ln/FSwL(N) 

612* BETE U,n)= ALPHA ( J#N)*SqEPRM 

' 613* BES0A( j»N)=B5SL(ALPHA{j»N)715 — - — — - - - 

614* BES03( J*N)=-3SSL(RETE( J,N) r 1) 

615* BtSI A ( J» N) zBSSL ( ALPHA ( J»N)73) 

616* BESI3< J»N)=BSSL(BETE( J,N) » 3) 

617* 52q~ CONTINUE' " ' 


616* 0062 1L=1 $ LE 

619* Ns 2 ' 

620* XSOO^O.O 

621* 622 XI=S^RT ( 1 , / ( <N**2) -1 ( 2, *A/FSWL (L ) ) **2) ) ) 

622* XSO= C (SIN(N*PI*WlOTH(I)/A) ) **2 ) * ( XI -1 . /N ) 

~£23* XSOO=XSO + XSOCT 

624* IF<x50.LE,,001)GO TO 624 

625* 623 N=N*1 ‘ ~ 

626 * 50 T ° 622 

627* 624 CONTINUE 

626* SOp( I »L)=XS00*2,-2 ,*( (sIN<PI*WIDTH(I)/A) )**2) 

629* 'CONTINUE “ 

630* _ S0PP = AL0G( (4,*A/(PI*DP( JJ ) *SIN(Pr*WlOTH< I ) /A) ) ) 

631* D0625Mzl/LE 

632* XB=(2, + a/WL<M) ) *( ( PI *D p ( J ) /A ) **2 ) * ( S I N { PI *Wl DTH ( I ) /A ) **2 ) / ( ( ( ALPHA 

633* ‘ 1( Jr^J ♦♦2)*BESIB( J»M)/BESIA(J»M) )*(l,/(ALPHA(JrM)*BESOA(JfM)* 

634* 2BESI3( J»M)*qETE( J#M)*be:S 08( JfM)*0ESlA(j,M) ) )-2, ) 

635* XA = .5*xb+(a/( 2,*WL(M) ) )7 ( 177 ( SIN { PI *WI0TH CIT/A3 _ r*2) T*TTBeSTmrJ7MTX 

636* 19ESOA ( J,M) )♦ ( l ,/( BETE (j » M ) *BES0A ( Jr M > *BESI0 < J * M ) -AlPHA ( J # M) * 

637* ~ 23ES03f J#M)*3ESIA{J,M)) j-(50PP+S0PfI»M))*(ALPHA(U,H)/4777 

638* Yl = l*/( C 0 r 1 ) *XB+1 • ) 

639* YH = l,/( (Or-l)*XA)+YI 

640* YT=l./( (U» 1 ) *X3>1,/YII) 

641 * BGa'M= ( Yt-1 . ) / ( 1 , +YT) 

642* GGAM< I r JrM,K)=8SAM*CExP< (Or 1>*2 ,*BEtA<M)*ALENGH<K) ) 

“ 643* IF ( ALENGH (K ) ,5E.O, ) GGAm ( ITJ iMiK )-BGAM 

644* AGAM<K,M)=GAM<M)*CEXP< ( 0 * -1 > *2, *BETA ( M) *ALENGH<X ) > 
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645* IF(ALENGH<K).LT.O.)AGAmOOM)=GAM<M) 

' 646* ~~ fi'AMMA (M>5{AGAM(K»Mr*GGAM(l7JVM;K>-r;‘*AGAMmMI*GGAMU7J'»MrK> >/( l.-“ 

647* 1AGaM<K#M)*GGAM(I.J.M.K) ) 

' 648* " 625 CONTINUE 

649* WR I t £(6«626) 'uifllCthiC kOsi Is El ) 

651* WRITE(6>627) ANCRMT 

652* 627 FORMAT {//# 1 OX r* NUMBER oF“ INCREMENTS" ’ iF4VO) : 

653* _ _ WRITE(6>628)FR{ JREG) .Fr(JENQ) 

654* 62Q FORMAT(//»lOX» ’MATCHED BAND IS FROM ♦ .F7. 3, * GHZ, TO "* .F7,3» * GHZ. 

65b* 1 • ) 

65b* WRITMG'p 629) WIDTH mTDpTjrt*CE>JGHTin 

657* 629 FORmAt(//» lflX» ’DISTANCE FROM SlOEwALL = • #F6.4, IX# ’CM, ’ , 5X. *DlAM£T 

658* ~ 1ER = ’ » F 6 • 4 » * CM 5X ,♦ DI STANCE TOWARD THE GEN, FROM MISMATCH = *# 

659* 2F8,4> 

660* C ************************************************************ *********** 

661* 999 CONTINUE 

' 662* WRITE ft* 22*8) 

663* 22fl FORMAT(//»iOXr ’RESULTANT MISMATCH* ) 

664* WRITE(6>229) 

66b* 22b FORM a T ( / / * 10X* ’FREQUENCY IN GHZ. * ’ 10* ’ • MAG. REFL. COEFF , • » lQX , * V$w 

666* 1R’ rlOxi «MA3.REFL. COEFF. DISURB.’ »//T"“ 


6b7* 00240N=1»LE 

660* “AMAS = C A3S(GAvMACN)*7 ~ 

669* _ VSWR=(1.+AMAG)/U.-AMAG) 

670* " ** AM»aG = CA3S(GAM{N) ) “ “ 

671* _ _ VSwRR=( 1.+AMMAG)/(1,-AmMAG) 

672*"” X(n)=VSWRR ~~ 

673* WRITE ( 6» 230) FR(N) ’ AM AG , VSWR » AMm AG 

'674* ■“ 23o FORMAT ( 10X»F 8.3 # 1 55f V F l o75VI2X7Pa • 47T0 x • FI 6 . 5 ) ~ 

675* _ __ 24o CONTINUE 

676*' WRlTC(6#23l) * “ '** " ’ 

677* 23 i F0 rmAT( 1H1,///,45X > ’MATCHING ELEMENT* ’ 3flX# ’UNMATCHED ELEMENT*) 

670*”' " WRl TC ( 6 ♦ 232 ) 

679* 23 2 FORMAT(//rlOX» ’FREQUENCY* ’10X» ’REAL REFl. COEFF, * » lOX* * I MAG REFL. 

660* r COEFF . ’ t iox* *REAE~REFiTr~tOEFF7’ >TDyV’ IMA^7"REFL, COEFF, * t//) 

681* D024lN=l»LE 

682* WRITE (6»233)FR(N) i GGAM ( IYJ»N»K) >GAM(NI 

683* 233 FORMAT ( iox»Ffl. 3’ 14X» Flo. 6’ 16X»F10.6» 15X’F10 .6* 15X#Fl0.6) 

684* 24i CONTINUE 

685* c Punch data cards. 


666* 

687* 


"D'Oi"001N=l.LE 
REl<N)=reAl<GAM(N) ) 



688* 

689* 

1001 

YImAG(N)=AIMAG(GAM(MTV 

CONTINUE 



690*“ 

691* 


DOio03N=1>LE 

REL(N)=REAl<GAMMA(N) ) 



692* 
6 J3* 

1003 

YImaG (N) = AI mag (GAMMA (Nfy 
CONTINUE 




694* w R I T E ( 1 • 1 0 0 4 ) ( R E LTN r# K)=um ) " 

695* _ WRITE( 1» 1004 ) (YIMAG(N) , N=LE ’ 1 # -1 >__ 

696*''"" 1004 FORM* T ( flFlo .6) ~ 

697* C Plot COMPARISON 8ETWEEN unmatched and matched disturbance. 

696* *' D'0'300on=1V"lE 

b99*_ __ Y<N)=FR(N) 

700** ~ 3000 CONTINUE 


7ol* 


R3 = xU) 
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702* DO3001N=1»lE 

703* IF { X C N) .GT # RB)R3iXCN) : 

704* 300l CONTINUE 

' 705* VLa»\|0=l # " " * 

706* VU6nD=RB 

“7(1 7* “CALL PLOT (id # ♦ -3iT^3) 1 

706* «EAL LTR52, LENGTH 

709 * logical BOX# FI XLIM# USETlTW 

710* DATA HTITLE/'FrEQUENCY in GHZ. 

'711 ♦ Z TT ; 

712* DATA VTITLE/*V0LTAGE STANDING WAVE RATIO 

713* 2 *7 

714* D03002NA=1 > 2 

715* IF(NA.EQ.1)G0' ‘"'TO" 3003 " ' 

716* D03004N=1#LF 

717* AR = CABS(GAM^A(N)')" ’ 

716* X (N) =( 1. +AR)/( 1,-AR) 

7 19* 3004 CONTINUE 

720* 3003 CONTINUE 

721* M0DE = 1 

722* IF (NA # E0«2)^0DE = 4 

723* LTRS?=.l “ “ " ~ ‘ ‘ " ‘ r " 

724* UPr4 • 5 

“725* OVer = OTO' 

726* HEigHT=4,5 

727* " ’ LENGTH = 4.5 

728* HLbnD=Y(1) 

729* HUBNDz y { LE ) 

730* ICo'jNTrLE 

“731* “BOx=.TRUE." 

732* USeLI’^.TRuE, 

733* FlxLlM=,TRUE. ' " * 

734* NU.M1 NT = 20 

735* CALL PLOT Y ( MODE * X# Y# icount» up# over* height# lengthvbox.fi xlimv usflim 

736* 1 #NU S1 INT# VTITLE»HTITLE,LTRSZ* VLBND# vUBND#HLBNO#HUBND#PERCNT) 

“ 737* 3002 CONTINUE “ ' ‘ 

738* CALL PLOTtO.O# 0.0#999) 

* 739* END 

END OF COMPILATION! - NO" DIAGNOSTICS 
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APPENDIX C 

DERIVATION OF THE IMPEDANCE MATCHING EQUATION 

The objective of this derivation is to establish the 
key equation which has been used in the impedance matching 
technique. Due to the comparison ref lectometers ' unique 
ability to locate discontinuities in an otherwise matched 
waveguide and to determine the reflection coefficient as a 
function of frequency for each disturbance, it will only be 
necessary to deal with two disturbances at a time, the ori- 
ginal mismatch and the matching element. 

The reflection coefficient measured and calculated 
in this thesis apply to the admittance Smith chart. The 
following derivation also applies to the admittance Smith 
chart. 

The following equations will be required: 

Phase propagation constant 

B = 2 7T / X (C-l) 

y 

Current reflection coefficient 

r io = l r i o l ex P [j 0 o ] 


(C-2) 
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Current reflection coefficient generated by the disturbance 
having the shunt admittance Y, where the characteristic ad- 
mittance is Y q 

F I = (Y-Y 0 )/(Y+Y q ) (C-3) 

Generalized current reflection coefficient 

r l = rio ex Pt"3 2BS 3 = | r Io jexp[ j0 o ]exp[-j2BS3 (C-4) 

Given a transmission line: 



Figure 91. Transmission Line with Disturbance 

The quantity Y a is a single shunt admittance in an otherwise 
uniform transmission line. The total admittance at point A 
is that formed by Y a and the characteristic admittance of the 
line, and is given by: 
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Y a = Y a +Y o 

This admittance defines the current reflection co- 
efficient Tj a at A by Equation (C-3) . (C-6) 

r la = Y a /< Y a + 2 Y b> 

The reflection coefficient at point B is given in 
terms of Tj a and L (the distance from A to B) as: 

r iab = r Ia expt-j2BL] (C-7) 

Similarly, if the only disturbance in the transmission 
line were a shunt admittance at B, the reflection co- 
efficient at point B is given by: 


r ib - V<V 2 V 


(C-8) 


Consider now, both Y a and Y^ in the transmission 
line as in Figure 89. 


Generator 

Side 



Figure 92 . Transmission Line Containing a Discontinuity 
at Point A and Point B 
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The total admittance at point A looking away from the gener- 
ator is given by Equation (C-5) . On the load side of Yj., at 
point B the admittance is given by : 

Y/Y© = (y a /y 0 +j tan 3L) / [1+j (y a /y Q ) tan 3L] . (C-9A) 

Defining normalized admittance y = Y/Y Q and y a = Y a /Y Q , 
Equation (C-9) can be expressed as: 

y = (y' +j tan BL)/(l+j y' tan £L) . (C-9B) 

d a. 

The total admittance on the generator side of point B is 
given by Equation (C-10) . 

(y t = Yb + (ya +1+ 9 tan 3L)/[l+j (y a +l)tan BL) . (C-10) 

Simplifying the above we have: 

y . = y a +yb+i+j (yb+yayb+Dtan bl 

‘ l+j (y a +l) tan El 

Substituting y^ into Equation (C-3) , the reflection co- 
efficient as pertaining to the admittance Smith chart is 
given by: 
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r = 2y a exp (- j 2 6L) +y a y b -y a Y b exp (-j2BL) +2y fc 
4+2y a +2y b +y a y b -y a y b exp (-j2$L) 


(C-ll) 


From Equation (C-6) and (C-8) , y and y can be expressed in 

a b 

terms of T and r_. . 

Ia lb 


y a " 

- 2r Ia /' r iB-D 

(C-12) 

*b = 


(C-13) 


Substituting (C-12) and (C-13) into (C-ll) and simplifying: 


F It= 


F Ia eXp ^ 2BL) +r ib -2r ia r Ib exp 2 6L) 


l-r ia r ib exp (-j26L) 


*C-14) 


Equation (C-14) expressed the total current reflec- 
tion coefficient at point B as a result of y^ and y^. This 
total current reflection coefficient is in terms of the re- 
flection coefficient T and r , which are defined by 

la lb 

Equations (C-6) and (C-8) respectively, r F , and F 

It Ia lb 

are reflection coefficient referenced to the admittance 
Smith chart. 


Expressing r in the form of Equation (C-14) is 
particularly compatable with the comparison ref lectometer . 
For example, for a positive distance L, r corresponds to 
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the mismatch measured by the comparison ref hectometer and 

F corresDonds to the matching element a distance L toward 
la 

the generator from the mismatch. The' total mismatch result- 
ing from this match is given by r ^ . The objective of this 

research is to reduce r_. to the lowest possible value over 

It 

the bandwidth of interest. 



APPENDIX D 


COMPARISON REFLECTOMETER OPERATING PROCEDURE 

The setup procedure is referenced to the reflec- 
tometer setup of Figure 23. In the following paragraphs 
it will be assumed that the operator has recently assembled 
the comparison ref lectometer and a complete calibration of 
the system is required. It will also be assumed that 
the individual components of the system are in good con- 
dition and meet their individual specifications. 

In order to reduce 60 Hz pov/erline ripple from inter- 
fering with the 1 KHz modulated r.f., it was found necessary 
to float the comparison ref lectometer above ground. This was 
done by plugging all of the units into the powerline using 
two-to-three prong adaptors, leaving the pig tails ungroun- 
ded. 

When energizing the comparison ref lectometer it is 
essential to turn on the coherent synchronizer first and then 
the BWO. This prevents possible damage to the BWO helix. 
Consequently, when shutting down the units the BWO should be 
shut down first followed by the coherent synchronizer. 
Coherent Synchronizer 

The LFE Model 244 should be set up using the follow- 
ing procedure when used in the comparison ref lectometer set- 
up . 

(1) Switch the front panel switch to multiply. 


Switch the reference mode to tunable. 
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(2) Switch the reference mode to tunable. 

(3) Switch the frequency tune to fixed. 

(4) Turn signal level approximately 1/4 revolution 
from full on. 

(5) Switch frequency fine tune to the 8-16 GHz 


range . 

(6) Switch the front panel switch to 50 MHz marker 
position. 

(7) Switch mode switch to 16-32 MHz/volt oscillator 
sensitivity . on back panel. 

(8) Connect 25 MHz tap TP1 to the input of a high 
impedance scope. 

(9) Connect power line plug to powerline stabiliz- 
ing transformer. 

(10) Switch on power to unit. 

(11) Turn phase lock error voltage to off. 

(12) Allow approximately one hour to warm up. 

This unit is now ready for the final operating procedure. 
Type 210 Sweep Oscillator 

The AIL BWO should be set up using the following 
procedure when used in the comparison ref lectometer setup. 

(1) Insert the 8.0-12.4 GHz plugin. 

(2) Push sweep select button to slow. 

(3) Turn rate inner knob clockwise, sweeper will 
stop after a single sweep. 

(4) Outer rate knob is used for fine rate adjust. 



(5) 


Push off button under AM. 


(6) Depress blank button. 

(7) Depress XTAL button. 

(8) Attach output from crystal detector on direc- 
tional coupler C2 to input of oscillator plugin, Figure 23. 

(9) Connect 200 micro farad external capacitor to 
range ext. BNC on back panel. 

(10) Plug in power line connector. 

(11) Depress line button and wait for the r.f. 
light to go on. 

(12) Switch r.f. knob on plugin to on position. 

(13) Adjust level control to obtain -10 dBm 

(.1 raw) power level at the output of directional coupler 
C3, Figure 23. 

The unit is now set up for final operating procedure. 

Audio Oscillator 

Set audio oscillator to 1 KHz. This can be done by 
connecting the output of the audio oscillator to the input 
of the Hewlett-Packard 5211A electronic counter and adjust- 
ing the oscillator for a 1 KHz counter reading. 

Turn the gain down on the audio oscillator and 
connect its output to the bias connector of the waveguide 
modulator. Connect a high impedance oscilloscope ' s input 
to the crystal detector on C3, Figure 23. Increase the 
audio oscillator gain until a square wave is observed on 
the scope. Care must be taken not to exceed the maximum 
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rating on the bias level of the modulator. 

Ratio Meter 

The Hewlett-Packard Model 416B ratio meter should 
be set up using the following procedure when used in the 
comparison ref lectometer setup. The reference step should 
be terminated in a matched load. 

(1) Turn the error voltage gain up until a phase 
lock is registered on the level meter of the coherent synch- 
ronizer. This will stabilize the frequency, and therefore, 
also the amplitude of the reflected wave. 

(2) Connect the incident and reflected inputs to 
detectors on Cg and respectively. 

(3) Plug in ratio meter power line using adapter. 

(4) Switch detector to XTAL. 

(5) Switch excess incident attenuation to 10 db. 

(6) Switch reference attenuator to -30 db. 

(7) Switch on power to the unit. 

(8) Allow unit to warm up for one hour. 

(9) Adjust reference adjust for a 100% scale 
reading. 

(10) Measure voltage at output using the Dymec volt- 
age measuring system. 

(11) Adjust output voltage, using the control on the 
back panel, to 6.928 volts. This voltage is Emax. 

(12) Readjust reference adjust to center of scale. 
This unit is now set up for the final operating procedure. 
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Voltage Measurement System 

The Dymec model DY-552A voltage measuring and record- 
ing system is comprised of the model DY-2210 voltage-to- 
frequency converter, the H.P. model 5211A/B electronic 
counter, the model DY-2540 scanner/coupler, and the motor- 
ized tape punch. This system should be set up according to 
the following procedure when used as part of the compari- 
son ref lectometer setup. 

(1) Switch on power to voltage-to-f requency conver- 
ter. 

(2) Switch range switch to 10 volt. 

(3) Connect output of ratio meter to input of volt- 
age- to- frequency converter. 

(4) Connect output of voltage-to-f requency converter 
to input of electronic counter. 

(5) Make proper zero adjustment and calibrate volt- 
age-to-f requency converter. 

(6) Switch on power to electronic counter. 

(7) Turn sample rate to the hold position. 

(8) Turn sensitivity half way. 

(9) Switch function switch to frequency/second posi- 
tion. 

(10) Switch on power to scanner/coupler. 

(11) Switch start switch to automatic. 

(12) Switch decay switch on back of panel to fast. 

(13) Switch back panel switch to tape. 
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(14) Plug in push button to remote reset according 
to the diagram on back panel of electronic counter. 

(15) Switch on tape punch. 

This unit is now set up and ready for the system final oper- 
ating procedure. 

Coherent Synchronizer Spectrum Test 

If a reliable phase lock is to be made at 50 MHz 
intervals, the coherent synchronizer must be in good oper- 
ating condition. It is important that the spectrum be 
relatively flat and that all phase lock points be achievable 
This is checked by tapping the 25 MHz signal from the 25 MHz 
amplifier and displaying it on a high impedance oscilloscope 
A representation of the varactor spectrum can be then ob- 
served by sweeping the leveled RF input of the BWO and 
observing the IF signal on an oscilloscope. A small wire 
was replaced in the coaxial varactor mount, shorting the 
center post to the grounded outer conductor. The distri- 
butive inductance and capacitance of this wire serves to im- 
pedance match the varactor to the varactor mount. 

The coherent synchronizers RF input was initially 
connected to the directional coupler by a type N adapter 
and a type N cable to the type N connector of the coherent 
synchronizer. It was experimentally determined that the 
frequency response of the cable interferred with the opera- 
tion of the varactor. A solid type N coaxial line re- 
placed the cable, and the varactor spectrum was then 
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adjusted so as to provide a flat harmonic spectrum from 
8.0 GHz to 12.4 GHz. These adjustments allow phase lock 
points from 7.975 GHz to 12.425 GHz. 

Adjust the lock range to zero on the front panel of 

the coherent synchronizer. Depress the fast button on the 

i 

AIL BWO sweeper and turn the sweep knob to the clockwise 
position for repetitive sweep. Use the sweep output of 
the swept oscillator to trigger the oscilloscope and depress 
the blanking button on the sweep oscillator plugin. 

Adjust the oscilloscopes' time sweep until one full 
sweep of the BWO fits within the scope. Figure 93. What is 
actually seen is the level of the 25 MHz signal as the rf 
to the input of the coherent synchronizer is swept. There- 
fore, the characteristics of the bandpass of the I.F. amp- 
lifier is reproduced every time the BWO frequency is within 




Figure 93. Display Showing the I.F. Band Pass as the RF 
Input is Swept * 
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25 MHz of a varactor harmonic, Figure 25. The level con- 
trol on the front panel of the coherent synchronizer should 
be adjusted to achieve as much gain as possible and yet to 
keep the noise level down. This can easily be done, since 
there is a point where increasing the gain increases the 
noise level sharply. The varactor should be adjusted such 
that the amplitude level of the peaks of the response of 
Figure 90 are relatively constant. Since this adjustment 
for the model 244 is a tedious one of trial and error, it 
is suggested that as much care as possible be taken with this 
instrument. 

System Operating Procedure . With all individual in- 
struments operating according to the previously discussed 
procedures, the system is now ready to be set up for measure- 
ments. 

(1) Depress the slow button on the sweep oscillator. 

(2) Turn the rate inner knob clockwise, to single 
sweep and stop position. 

(3) Turn outer knob half way. 

(4) Attach the test element to the comparison re- 
flectometer as in Figure 23. 

(5) Depress sweep button on sweep oscillator. 

(6) While oscillator is in the process of a slow 
sweep, monitor the meter on the ratio meter and adjust re- 
ference control to ensure that the needle remains in the 
range of the meter scale. Steps 5 and 6 may need repeating 



so that in a single sweep, without changing controls, the 
needle swings symmetrically on the scale. 

(7) Reduce the sweep rate such that it takes approx- 
imately 15 minutes to sweep the band one time. This time 
can be reduced as the operator gains experience with the 
equipment. 

(8) Set the oscilloscope to internal trigger and set 
the sweep rate to maximum. Adjust the scope gain such that 
the 25 MHz wave will half fill the screen when a phase lock 
has occur ed. This can be done by turning the lock range to 
maximum. When the scope adjustment has been made return 
the lock range to zero. 

(9) Adjustment of the lock range determines the 
phase lock duration. During the phase lock duration a mea- 
surement is taken with a normally closed push button. On 
the particular model 224 used in this research the control 
was set approximately 1/64 inch above the zero volt line. 
This setting is best made after working with the control 
while making a sweep. 

(10) Depress the front panel start button. This 
starts the slow sweep. As the frequency increases the wave- 
form will increase on the oscilloscope and the level on the 
level meter will increase. When a phase lock is made the 
phase of the sine wave on the oscilloscope will change 
phase by 180° and 'will remain very stable in amplitude and 
phase. The level meter will jump as the phase is locked 



and remain at that level. The phase meter on the coher- 
ent synchronizer which is an indication of the error volt- 
age on the helix of the BWO, will begin to move, as the 
internal sweep attempts to change the BWO frequency. At 
this time a measurement is taken by depressing the hand 
held push button, Figure 23. When the error voltage is at 
its maximum setting and can no longer hold the BWO, the 
phase lock is broken. The phase meter centers, the level 
drops to zero and the sine wave disappears from the scope. 
The BWO jumps to a higher frequency and continues to in- 
crease slowly until again phase locked automatically. This 
is done 89 times and gives one set of measurements. 

The second set of measurements is taken with the re- 
ference step terminated in the matched load, Figure 23. The 
procedure is the same as that described in the above para- 
graphs. It is important not to change any adjustment on 
the ratiometer or measuring equipment during the two sets 
of measurements, since the two will be compared during the 
calculations (III-31). 

When measurements are completed, the information on 
the tapes is punched on data cards by a tape to card con- 
verter. These data cards are then used as data in the re- 
flectometer computer programs, Appendix B. 
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